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ABSTRACT 


This report presents the results of experimental and 
numerical investigations of the flow field in the head-end star 
grain slots of the Space Shuttle Solid Rocket Motor . This work 
provided the basis for the development of an improved solid 
rocket motor ignition transient code which is also described in 
this report . The correlation between the experimental and 
numerical results is excellent and provides a firm basis for the 
development of a fully three-dimensional solid rocket motor 
ignition transient computer code. 
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NOMENCLATURE 


b 

c x ,c 2 

c p 

e 

e gp 

K 

k 

K 

M 

Nu 


P 

Pr 

r 

R e 
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t 

T 

u, v, w 
V 


x,y, z 
a 

y 

e 

o k , c e 

<*p 


subscripts 


= star slot width 

= empirical constants appearing in turbulence model 
= specific heat 
= gas internal energy 

_ internal energy of the gaseous burning propellant 
= convective heat transfer coefficient 
- kinetic energy of turbulent velocity fluctuations 
= thermal conductivity 
= Mach number 
= Nusselt number 
= pressure 
= Prandtl number 
= propellant burn rate 
= Reynolds number 

= source term(s) in governing conservation equations 
= time 

= temperature 
= velocity components 
= (u 2 + v 2 ) w 
= spatial coordinates 
= thermal diffusivity 
= ratio of specific heats 

= rate of dissipation of turbulent kinetic energy 
= constants appearing in turbulence model 
= constant appearing in propellant burn rate equation 


c = mass 

e = energy 

i = initial value 

l = laminar 

P = propellant 

ref,R = reference condition 

star tip = propellant grain star tip 


t 

wa 

xm, ym 


turbulent 
adiabatic wall 

x, y components of linear momentum 


superscript 

* 


dimensionless quantity 
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I . INTRODUCTION 


This report describes the results obtained from experimental 
and numerical analyses conducted for the purpose of developing an 
improved solid rocket motor ignition transient code. The 
specific objective of this work was to improve ability to predict 
the influence of the star grain on the ignition transient; in 
particular, the calculation of the flame spreading rate on the 
propellant surfaces inside the star slot,. This report is divided 
into three main subject areas: 1) experimental analysis, 2) 
numerical analysis and 3) computer code modification, development 
and documentation. 

Section II presents the results of a series of experiments 
conducted at NASA's George C. Marshall Space Flight Center 
(MSFC) . These experiments utilized a model designed and 
constructed in the Aerospace Engineering Department at Auburn 
University. This model is a one- tenths scale simulation of the 
Space Shuttle Solid Rocket Motor ( SRM) head-end section. The 
model was tested in the special test section of the 14 xl4 
trisonic wind tunnel at MSFC. The tests were cold-flow, using 
air, simulations of the internal flow through the igniter nozzle 
and the head-end section of the SRM. The air was supplied 
through a special high pressure line connected to the wind 
tunnel. There was no external flow a r ound the model. The tests 
were designed to provide both qualitative and quantitative data 
on the interaction between the igniter plume and the star slots 
and the flow field within the star slots. Qualitative 
measurements were made using oil smears, Schlieren photography 
and by seeding the flow field with aluminum particles which were 
illuminated with a laser system and recorded on video tape . 
Quantitative measurements were made of the pressure distribution 
within a slot and of the heat transfer rates to the wall of a 
slot. The correlation of the data between the various 
experiments performed is excellent. 

Section III provides the theoretical basis for the 
Computational Fluid Dynamic (CFD) model which was developed to 
analyze the flow field in a star slot and the results obtained 
from the CFD analysis. The CFD model was verified using the 
experimental data obtained from the cold flow tests described 
above. The primary objective of the CFD model was to provide 
data regarding the spread of the combustion process throughout 
the star slots. Based on these calculations a burning surface 
area versus time model is obtained for use in the ignition 
transient performance model of an SRM. The results of the 
analysis are compared to existing motor data and are shown to be 
in good agreement . 

Sections IV- VII describe the modifications, interface, _ 
program description, and instructions for an improved ignition 
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transient computer code. The modifications and interface are to 
an existing one - dimensional ignition transient computer code. 

The code was chosen because it is well documented and provides an 
efficient means of implementing the modifications to the burning 
surface versus time calculations as determined from the analyses 
described above. The ability to account for the flame spreading 
in the star slot of the head-end section represents a significant 
step forward in the ability to accurately model the early portion 
an SRM's ignition transient.. The last two sections ^ describe the 
development of a pre — and post-processor for an ignition 
transient code. The code presented here uses a modified one- 
dimensional ignition transient code modified to account for the 
flame spreading in the star slots to solve for the ignition 
transient performance of an SRM. However, it would be relatively 
straightforward to modify the current version of the code to use 
a more sophisticated ignition transient model as a solver. The 
computer code is currently implemented on a SUN workstation using 
X-Windows. Instructions necessary for using the program in an X- 
Windows environment are included in Section VII. A sample input 
file and the code for the pre- and post— processor are included in 
the Appendices . 

The work described in this report provides an excellent base 
for the development of a fully three-dimensional ignition 
transient performance prediction code using CFD techniques. The 
experimental database which has been generated from this work 
provides significant new insight into flow field phenomena 
occurring in the star slots of SRM's which have head-end star 
grains and head-end igniters. 


1-2 



II. EXPERIMENTAL ANALYSIS 


One of the primary limitations of existing ignition transient 
prediction computer codes is an overestimation of the ignition 
delay and departure of the predicted pressure-time history curve 
from measured data in the early part of the transient. One of the 
reasons for this is the lack of data on the flow field m star 
grain sections of solid rocket motors. The experiments described 
in this section focused on the flow field in the star gram section 
of the Space Shuttle solid rocket motor (RSRM) . The purpose of the 
experiments was to obtain a credible data base for the f low field 
patterns and heat transfer rates within the star slots. A one 
tenth scale model of the Space Shuttle RSRM head-end star grain was 
designed and constructed for the experimental study. ?he 
availability of experimental flow field data during the ignition 
transient in solid rocket motors is very scarce. Conover 
conducted cold flow tests using a one-tenth scale model of the 
Space Shuttle solid rocket motor's head-end star grain section. 
Conover's tests used both a single port igniter, such as found on 
the RSRM and a three port igniter. The series of tests included 
Schlieren photographs of the igniter mounted in a plenum, oil smear 
data, pressure data and heat transfer coefficient measuremen s. 
Fifteen pressure ports were fitted in one side plate of a slot, 
while fifteen calorimeters to measure heat transfer coefficients 
were located in another slot. The data were taken at pressure 
levels from 100 to 1500 psi in 100 psi increments. The actual 
Space Shuttle igniter operates at approximately 2000 psi. Limits 
on the test facility prevented taking data at higher igniter 
chamber pressures. The static pressure data obtained provide some 
qualitative trends, but there was considerable scatter m the data 
when the igniter chamber pressure exceeded 1100 psi. This was 
probably due to the fact, as Conover states, that "above this 
pressure the side plates used to form the star grain slots were 
deflected to produce a one-sixteenth- inch gap where the plates come 
together at the star points," thus causing some leakage. There 
were also some apparent inconsistencies in the temperature data due 
to model warm-up during the test. However, the oil smear data 
provide a good indication of the recirculating flow pattern inside 
a slot. Even though Conover's experimental data provide useful 
information on the flow field inside the star slot, they should be 
considered preliminary in nature, and a starting point for a more 
in-depth investigation. 


The present work describes a series of tests directed toward 
the collection of qualitative and quantitative data documenting the 
main characteristics of the flow i n the head-end section of a solid 
rocket motor. In particular, the following objectives were to be 
accomplished with the test program. 
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1) Obtain information on the igniter plume structure and 
shape and its interaction with the star grain geometry. 

2 ) Determine the region of the igniter 

the side walls of the slots of the RSRM star grain 

section model. 

3) Determine the flow field characteristics of the subsonic, 
recirculating flow within the slot. 

4) Measure the heat transfer coefficients at several 
locations inside the slot. 

The above objectives were to be accomplished using the 
following techniques: 


1) Flow visualization. 

2) Oil smears 

3) Static pressure measurements 

4) Heat transfer coefficient measurements 

5) Velocity measurements 


Th ~ test article used in this investigation was a one-tenth 
scale h cold flow model based on the geometry of the Space Shuttle 
head-end section. The test article had four slots, as opposed 
to^eleven slots in the actual Space Shuttle motor. A single port 
igniter model and two four-port igniter models were used m the 

teStS m he scale factor of 1:10 was derived from an analysis which 
matches the Reynolds number between the model flow and the full 
scale flow in the star grain section. Besides geometrical 
similitude! the primary scaling Parameter to 

cimilaritv between the cold flow model and the real, full scale 
moSr is the Reynolds number. Compressibility effects are 
• mnortant onlv in the igniter plume region, since the flow inside 
Ts?ar slot is essentially subsonic. However, the igniter mass 
?low ratf is an imp” rtant parameter, since it is thought to be 
responsible for entrainment of the flow, which determines the 
^Srculftion pattern. The value of the Reynolds number determines 
££ l ”tu?e Sf the viscous effects. The viscous effects in turn are 
the mature or tne vi* convective heat transfer coefficient and 

therefore to the amount of heat transfer from the hot gases to the 
solid oropellant. An exact match of Reynolds number between the 
ii? scale flow field is not necessary for similarity. 

Instead generally good agreement between Reynolds numbers is a 
sufficient* condition^ for the general studies of these flow fields, 
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as has been extensively documented in the 
literature . 2 ' 5 


Because of the complicated nature of the flow in the head end 
section, it is difficult to define an overall Reynolds number for 
the entire slot. However, it seems appropriate to consider a 
representative Reynolds number, which can be defined at any point 
in the star slot, as: 


J? = 


e 


pVL 

U 
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where p, V and p are the local density, velocity and viscosity, 
respectively, and L is a geometric reference length (for instance, 
the P distance from a given point to the motor centerline) . The 
product 0J can be viewed as the local mass flux. It is assume 
that this local mass flux is proportional to the overall mass flux 
that enters a star slot, which in turn is proportional to the mass 
flow rate that enters the star slot divided by the open area of the 

slot, so that: 


til 


pV = 


slot^act 

bl 
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where b and 1 are the slot width and slot length, respectively, as 
shown in Figure II-l. 

The mass flow rate that actually enters a given slot depends 
on several factors, such as igniter mass flow rate, ih. g , geometric 
dimensions of the slot and motor, and igniter plume shape. The 
functional dependency between these parameters is usually not 
known, but can be generally expressed as: 


tii 


tii 


dimensional effects are negligible and the mass flow available to 
the slot ih , is the portion of igniter mass flow encompassed 
in the circular “sector facing the same slot (see Fig. II-2) , or 
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simplify the analysis, it can be assumed that three- 


0 b 

^siot, avail c ^^211 * iff 2nr 
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The factor b/(2rr) can be considered as the fraction of the port 
perimeter occupied by a slot. 

Using Eqs . (II-2 ) , (II-3) and (II-4) , the Reynolds number can 

be expressed as: 
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R * 
6 


— flplumeshape ) 
2nr l\i 


Note that the Reynolds number, and thus the flow pattern inside the 
slot, is independent of the number of slots. 


For proper scaling, let R el = R e2 , where the subscript 1 refers 
to the actual motor and 2 refers to the model. Then the following 


assumptions are made : 


1) All dimensions are exactly scaled. 


2 ) 


The igniter plume shape is scaled for 1 and 2. 


Note that the plume shape is determined mainly by the igniter 
nozzle area ratio, igniter stagnation pressure, P q , and back 
oressure. When all dimensions are scaled, the area ratio is the 
same for the real srm and the scaled model. The back pressure is 
also the same, and equal to the ambient value. This latter 
condition is true at least until the point in time when a first 
ignition occurs. This suggests that the igniter total pressure 
ideally should be the same, that is 2000 psi. 


Under the assumptions above, f (plume shape) drops out when 
equating the Reynolds numbers, and the following relationship is 

obtained: 


hhh = h 

r i\h~ A igl ^2 

Then, defining the scaling factor as S, so that 


h h . fi 

" h k r i 
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Equation (II-6) is reduced to: 

S = ^ II-8 

* igl ^2 

At this point, an appropriate viscosity- temperature relationship 
must be chosen. Following Caveny 5 , an empirical curve fit for 
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typical combustion gases is such that viscosity, y, is proportional 
to the gas temperature, T, as 


so that Eq. (II-8) becomes: 


S = 




i stf 


iti 


igl 



0.6 
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qince Ih T and T are known values, it appears from Eq.. (11-10) 
that the 1 ' scale factor, S, can be arbitrarily chosen simply by 
varying the igniter model mass flow rate, However, yZ 
igniter throat area is dependent on the scale factor to be 
determined. Furthermore, the igniter chamber pressure should be 
close to 2000 psi because of the assumption that the plume shape 
does not change from the real motor to the scaled model. This is 
more evident if Eq. (11-10) is rewritten in terms of stagnation 
pressures and temperatures, as: 


S = 


01 


^2 T 02 ^ ^2 ) 0.6 


*0 2 \| *2 T 01 T 1 N Y 2 2 


( 


Yi +1 


Yl± 

2(y 1 -D 


w 1 

2(Y 2 -1> 
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y 2 +1 


Note that the term A\/h\ is, itself, S 2 . 

For the cold-flow tests, the values of R 2 and y 2 for air are 
chosen for the model, and typical values of the temperatures are 
T = 3 10°K, T = 298°K . Since P„ = P 02 = 2000 psi, the pressure 

term does not Contribute to the value of S. Representative values 
of the variables for the flow in the head-end section of the actual 
SRM are taken from Ref. 5. 

Substituting in Eq. (II-ll), a value equal to 0.098 is 
obtained for S, which is very close to the chosen scale factor 
1:10. Fortuitously, a one-tenth scale also defines a model size 
which is close to the largest size that would fit in the dimensions 
of the test section of the 14xl4-inch tunnel employed for the 

tests . 

The number of slots was determined based on different 
criteria. According to the previous analysis, the Reynolds number 
is independent of the number of slots. Therefore, one can ideally 
choose any convenient number. Even though the actual SRM head-end 
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section has eleven slots, only four slots are used in the scaled 
model. This represents a tradeoff between the desirabi y 
visualization in one (transparent) star slot and the requirements 
for all other test instrumentation. 

The entire star grain section model, as well as the three 
ianiter models, were fabricated from aluminum, except for th 
transparent slot which consists of two plexiglas plates. The total 
length of the model is 19.72 inches; the largest diameter is 16 
inches. Figures II-3 and II-4 show a schematic representation of 
the entire model and of the igniters, respectively. 

Each star slot is formed from two plates separated by a bottom 
spacer of suitable thickness. An insert at the head-end simulates 
the actual grain surface. The circular port 13 fo ^ ed fr °^ 
contoured pieces connected to the outer surface of the slot plates 
Two plates are located at the upstream and downstream end to 
the slots, and provide the attachment points for all the parts. 
The igniter is connected to the inner surface of the head-end 
plate 5 The single port igniter has a insert into which the nozzle 
Sfbeencut V/single port igniter has a throat drameter of 
0 6025 inches and a conical shape with a half-angle of 27.2 
degrees. Sr area ratio is 1.428. The four-port igniters are 
simply four straight holes drilled in the igniter casing, each with 
a 0 3012 diameter. The four-port igniters are .oriented so that 
each of the four jets centerlines are directed into a star slot. 
Figure II-5 shows the whole star gram section model. The three 
igniters used are shown in Fig. II-6. 

The model is fully instrumented for measuring the parameters 
of interest as defined above. 

Static pressure measurements are taken inside a single star 
slot and along two of the four contoured sectors 
circular port of the model. Three pressure ports are located along 
each of the two contoured sectors . Twenty-eight static pressure 
ports are provided in one wall of a star slot. _ The measurements 
are taken at three different slot depths and consist of eight, ten 
and eight pressure taps, respectively, as shown m Fig. II -7. Tfte 
three depths are equally spaced along the height of the slot. 

In the plate forming the wall adjacent to the one containing 
the pressure taps, twenty-eight calorimeters are installed 
obtain heat transfer coefficient data. The calorimeters are placed 
at the same geometric locations as the pressure ports (Fig. II 9) • 
Each calorimeter is mounted in a plug, flush with the inner surface 
of the slot plate. A detail of the calorimeter installation is 
aiven in Fig II-8. In order to get accurate measurements of the 
heat transfer coefficients, the calorimeters were pre-heated before 
each measurement was taken. 

A second star is used to obtain oil smear data. A silicone- 


II-7 





Fig. II-6 Igniter models 
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Fig. II-7 Location of pressure taps and calorimeters 



Fig. II-8 Detail of calorimeter installation 
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based oil was used to apply a matrix of oil drops to the finished 
guxfcice of one of the plates in the slot. After each run, well 
defined marks or smears indicated the local direction of the flow 
and gave a good overall picture of the flow field. 

The transparent slot previously mentioned was used for real 
time flow visualization. A laser sheet was projected from the 
aft-end of the model and illuminated most of the transparent 
slot Aluminum particles mixed with pure alcohol were injected 
into * the slot and the aluminum particles, were illuminated by the 
laser sheet. The movement of the particles was clearly visible 
in the transparent slot and provided an excellent qualitative 
measurement of the behavior of the flow field in the slot. The 
flow visualization obtained using the aluminum particles gave a 
more detailed picture of the flow patterns in the slot than was 
possible with the oil smears. In addition, the real time nature 
of measurements provided a means for studying the dynamic 
characteristics of the flow field. Video tape recordings of 
these experiments were made to document the measurements . 

The fourth slot was initially intended for making hot-wire 
aneometry measurements. However, because of difficulties 
associated with flow blockage in the small slot and difficulties 
in making measurements in or near the high speed (high subsonic 
or low supersonic) plume, this measurement was abandoned for the 
present investigation. 

Test Facility 

The cold-flow tests used the special test section in the 
NASA Marshall Space Flight Center X4xl4 -inch trisonic wind 
tunnel . The tunnel operated as an intermittent blow-down wind 
tunnel from storage pressure to atmospheric exhaust. The full 
Mach number capability was not needed for the test program which 
was carried out . Instead only a high pressure internal flow 
through the special test section was required. The high pressure 
air passed through the hollow centerbody of the tunnel, into a 
pipe connected to the head-end plate of the model. It was then 
exhausted into the special test section through the igniter 
models. There was no external flow around the model. A venturi 
was installed upstream of the model to determine the mass flow 
rate through the igniter. Additional information regarding the 
NASA/ MS FC trisonic wind tunel is given in Ref. 6. 


Test Plan 

The test program included two main series of tests. 

Table II-l shows the test matrix which was used for each series 
of tests. In the first series each igniter model was placed in 
the test section without the star grain portion of the model. 

Air flow at pressures of 100, 500, 1000, 1500 and 1800 psi passed 
through each of the igniters used in the experiments and mass 
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flow rates corresponding to each pressure were recorded. A 
Schlieren system and video tape recorder were used during each 
run to examine the plume shape at various pressures. This was 
done to establish a reference for the plume geometry which could 
be compared with the plume geometry observed with the star gram 
in place . 

The second series of tests used the entire head-end section 
model along with each of the three igniters. Oil smear, flow 
visualization, static pressure and heat transfer coefficient 
measurements were made at each condition shown in the test matrix 
of Table Il-l. It should be noted that the test condition at 
1800 psi approximates the design condition of 2000 psi at which 
simlitude between the one-tenth scale model and the actual Space 
Shuttle RSRM flow field is achieved. The value of 1800 psi was 
used because it represents the upper limit on the facility at the 
mass flow rates necessary for the tests. 


Table II-l. Test matrix 


Igniter 

Angle 

(deg) 

Chamber 

Pressure 

(psi) 

100 

500 

1000 

1500 

1800 

0 


1 

2 

3 

4 

5 

22.5 


6 

7 

8 

9 

10 

45 


11 

12 

13 

14 

15 


Oil Smear Test Results 

Oil smears were taken for each of the fifteen test 
conditions shown in Table II-l. The oil smears were generated 
from a pattern of oil droplets placed on one side of a slot. The 
spacing between the droplets was approximately one-half inch. 
Figs. II-9 through 11-23 show photographs of the oil smears 
generated for each of the fifteen test conditions. The oil 
smears provide considerable detail regarding the direction of the 
primary flow in the slots, the region of the igniter plume 
impingement, and the recirculation patterns which occur in the 
slot Although qualitative in nature, this data showed good 
agreement with the data from the CFD analyses presented in refs. 
7 _ 9 , which are summarized in Section III of this report. 
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Results From Static Pressure M easurements 

Static pressure measurements were made at 27 locations . on 
the surface of one of the slots. Fig. H-7 shows the location 
and numbering scheme used for the static pressure ports . Static 
pressure distributions obtained from these measurements are shown 
in Figs. 11-24 through 11-38. This data confirms the _ 
Quantitative data obatained from the oil smear tests with regard 
to the location of the main, flow paths, recirculation regions, 
stagnation points and "dead vXegipns ® thin the slot The data 
obtained agrees well with theCFD results which will be discussed 
in Section III, where a comparison will be given between the 
experimetal data and the CFD results . 

Results From Heat Trans fer Measurements 

Calorimeters were placed in a slot face adjacent to the slot 
face where the static pressure measurements were made. The 
calorimeters were located at points corresponding to the 27 
locations shown in Fig. II-7. Because of the lack of a 
sufficient number of calorimeters to measure data at all 2 / 
locations simultan eous ly, two runs were made using 15 
calorimeters in each run. T hre e calorimeters were not moved 
between runs to insure that consisten t data were being obtained 
between the two runs. The r esults for the measured heat transfer 
coefficients are shown in Figs . 11-39 through 11-53 . 
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Fig. 11-10 Igniter 2 (100 psi) 




Fig. 11-11 Igniter 3 (100 psi) 
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Fig r* 11-14 Igniter 3 (500 psi) 
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Fig. 11-16 Igniter 2 (1000 psi) 



Fig. 11-17 Igniter 3 (1000 psi) 
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Fig. 11-20 Igniter 3 (1500 psi) 


11-17 




Fig. 11-21 Igniter 1 (1800 psi) 



Fig. 11-22 Igniter 2 (1800 psi) 



Fig. 11-23 Igniter 3 (1800 psi) 
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Fig. 11-24 Igniter 1 (100 psi) 

Measured Pressure, 22 deg Igniter 
Planter = 1 00 psi 



Fig. 11-25 Igniter 2 (100 psi) 


Measured Pressure, 45 deg Igniter 
Plgnter =100 psi 



Fig. 11-26 Igniter 3 (100 psi) 
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Measured Pressure, 22 deg Igniter 


Plgniur = 500 psl 



Fig. 11-28 Igniter 2 (500 psi) 


Measured Pressure, 45 deg Igniter 
Plgnlttr = 100 pal 



Fig. 11-29 Igniter 3 (500 psi) 
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Measured Pressure, Single Port Igniter 



Fig. 11-30 Igniter 1 (1000 psi) 


Measured Pressure, 22 deg Igniter 
Plgnltw ■ 1 000 pSl 



Fig. 11-31 Igniter 2 (1000 psi) 


Measured Pressure, 45 deg Igniter 
Plgnltar * 1 000 psi 



Fig. 11-32 Igniter 3 (1000 psi) 
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Measured Pressure, Single Port Igniter 
. Plfli»H*r e 1 500 pel a 



Fig. 11-33 Igniter 1 (1500 psi) 


Measured Pressure, 22 deg Igniter 


D. ^ 1 5AD nil 



Fig. 11-34 Igniter 2 (1500 psi) 


Measured Pressure, 45 deg Igniter 



Fig. 11-35 Igniter 3 (1500 psi) 
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Fig. 11-37 Igniter 2 (1800 psi) 



Fig. 11-38 Igniter 3 (1800 psi) 
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Uuiurid HMt Tran afar, Sing l* Port Ignllar 
PorVi*r = 66 atm (100 pal) 



Fig. 11-41 Igniter 3 (100 psi) 
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Fig. 11-42 Igniter 1 (500 psi) 
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Fig. 11-43 Igniter 2 (500 psi) 


NMtwtd H*at TWuhr, W 

%rfc-.*0Cp»l 



Fig. 11-44 Igniter 3 (500 psi) 
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Mooturod Hoot Tronofor, Single Pori Ignllor 


Plfrrflar ■ 00 Stm (1 000 pit) 



Fig. 11-45 Igniter 1 (1000 psi) 


NiHind He* Trencfer, JU* loiter 
IV**. 1000 pel 



Fig. 11-46 Igniter 2 (1000 psi) 


Measured Heat Transfer, 43° Igniter 
Planner = M atm (1000 pel) 



Fig. 11-47 Igniter 3 (1000 psi) 
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amnd HmI Tr»n*l»», Unfl 1 * r °* >«">*•» 
Pi,^ ■ 1M0 p*l 



Fig. 11-48 Igniter 1 (1500 psi) 


HM«nd HMt TtWMtW, Mi* ientUf 
rVMr.lSOOpal 



Fig. 11-49 Igniter 2 (1500 psi) 

HtttHnd HmI Trantfer, «• IflnH* 

(VMr.HOOpal 
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Measured Heat Transfer, 8lngle Port Iflnllar 
Plant* • 122A aim (1800 pal) 



Fig. 11-51 Igniter 1 (1800 psi) 





Fig. 11-52 Igniter 2 (1800 psi) 


Measured Heat Transfer, 45 # lgnHer 
Fiona*. 12fc4 atm (1800 pel) 



Fig. 11-53 Igniter 3 (1800 psi) 
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III. THEORETICAL /NUMERICAL ANALYSIS 


Introduction^ ion transient of an SRM emp ioying a pyrogen igniter 

can be defined as the time interval from the onset of the igniter 
flow to the time a quasi-steady flow develops. The starting 
transient is traditionally divided into three phases: the induction 
interval, or ignition lag; flame spreading; and chamber filling. 
The induction interval begins when the igniter flow is initiated 
and ends when a point on the propellant surface reaches a critical 
ignition temperature, and a flame first appears. The flame 
spreading phase follows, ending when the entire propellant surface 
is ignited. Following this is the chamber filling phase, during 
which rapid chamber pressurization occurs due to the energy and 
mass addition from the burning propellant. A peak pressure may 
occur followed by a pressure decrease towards an equilibrium 
value! attained when mass production by the propellant equals the 
mass outflow from the motor nozzle. Numerous studies have been 
directed at the analysis of SRM ignition transient phenomena. As 
discussed by Peretz. et al. 1 , these analyses can generally be 
categorized into three groups: (1) lumped chamber 

models 2 ' 7 ; (2) one dimensional, quasi-steady flow, P x) mod si* s^, 
and (3) temporal and one-dimensional flow field, p (x, t) models . 

The simplest analyses fall into the first group; a uniform chamber 
pressure is assumed and an equation for dP chmber /dt is derived an 
integrated to obtain a pressure-time trace. The flame spreading 
speed is assumed to be a known constant. In P(x) type models, flow 
property distributions are considered at each instant of time and 
one-dimensional steady state conservation equations are solved 
along the motor axis. As in the P(t) type models, the flame 
spreading speed is not part of the solution but rather must be 
input. In P(x, t) type models both spatial and temporal property 
variations are considered. A series of control volume increments 
are assumed along the motor axis and a set of time dependent one 
dimensional conservation equations are solved. The flame spreading 
speed can be obtained as part of the solution if convective heat 
transfer to the propellant grain is taken into account A widely 
used model of this type is that developed by Caveny and Kuo. 

Jasper Lal, et al 12 , have recently developed a one-dimensional 
model which takes canted pyrogen igniters into account by 
modifying the heat transfer analysis to include direct igniter 
plume impingement on the solid propellant surface . Even more 
recently, Johnston 13 has presented a numerical procedure for the 
analysis of internal flows in a solid rocket motor wherein an 
unsteady, axisymmetric solution of the Euler equations is combined 
with simple convective and radiative fluid heat transfer models and 
an unsteady one-dimensional heat conduction solution for the 
propellant grain. Flow in the star grain slots is not directly 
calculated. Instead, burn rate constants are adjusted to account 
for the variable area in the star grain region. Results are 
nresented for a Titan 5% segment SRM, a Titan 7 segment SRM, and 
the Space Shuttle SRM. Of these three motors, the Space Shuttle 
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cpM has the more pronounced axial grooves in the star grain, and 

611 in general suSfas 1 ' tfoTe* usifff 

with “he exception o£ the time period which 

S?e|ly "prlmlrUy 

factors : ar 0- U ) the th | l0 s W tar i g 1 eometry “ifthe hJaTetd segment °f 

tpproximS?ed by variations in port area and burning perimeter of 

nra in- and (3) the igniter flow field is not taxen into 
account The present analysis seeks to address these issues. 

Conservation Equations 

Tn this investigation, the Space Shuttle solid rocket motor 
(qrm> is taken as the reference motor design. It is characterized 
<fflaroe length-to-diameter ratio and by a. small port-to-nozzle 

throat area ratio. m0t ?he h^ad-e^d? sta” °shaped 

rS^oIfof the solid propellant grain co^ains p^TlI- ^7 * Cr ° SS 

SeCti ?L°fflow e f h ie a i d d" 6 ?o re^na'lyzld 8 i°s extremely complex; it is 
unsteady, multi-dimensional, turbulent, and compressible . J^th* er, 
the flow field is divided into a supersonic core region, defined by 

inside th^stfr^ unsteady compressible 

cylindrical fordfnate ^ystemfor the port region from the motor 
centerline to the star grain tips, and a rectangular cartesian 
coordinate system for the region of the flow inside the star si . 

Although the flow- field is three dimensional in the head end 
section of the motor, the star grain cross sectional shape of the 

propellant segment “plies “o "restrict the domain to 

S y ^ing7e Heitor, as shown in Fig. III-2. The two-dimensional 
Navier-Stokes equations are obtained by averaging the full^ 
three-dimensional ^er-Sto^s ndi |Sf ff the pla'ne of sy^etry 

ff Se star sfot. Avefaginf fs carried out along the azimuthal 
^inai-e A in the port region, and in the z-direction inside 
thS star slot.' The calculated flow field variables thus represent 
an average value in the domain considered (area shaded in Fig. Ill 

2) * The governing equations can most conveniently be solved in 

existl° W 

S€htk%itn s a i e^ef f StK 3fS&SU 

11% rtfmofr centerline to the bottom of the slot taken as the 
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reference length. It should be noted that the igniter mass flow 
is itself a function of time; this will be discussed in more detail 

iater^he dimension i es s governing equations then take the following 
form: 

continuity: 


dp* + dp*u*^apV* = s III-l 

at* dx * dy* 


x-momentum: 

ap*u*^P*u* 2 , a P *uv: 
at* dx* dy* 


ap*^ i dp* r ^ §u*__ 2 av*j 
dx* Re dx* 3 dx* 3 dy* 

III-2 


Re dy* 


r du*^dv^ ] + 

ay* ax* 



4 d^u* , dy* •) 

3 3x* 2 ay* 2 3 dx*dy* 


+ 


y -momentum: 

dp*v* + dp*u*v* . dp*v* a 

at* ax* ay* 


ap*^ 1 dp* r 4 dv* _ 2 au* j 
lZe ay* 3 dy* 3 dx* 
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Re dx* dy* 




4 ffv* , dV* ^ 1 d ?u* j 

3 ay* 2 ax* 2 3 ax*ay* 






energy: 

3p*e ^ dp*u*e* ^ dp*v*e* _ 

at* dx* dy* 


— — r 1 

Y~1 Px Re M 2 


dx* 


(K*^) + 

dx 


dy 


... V dy * 11 * Ar* av‘ Se 


ax* ay* 


au-.j 2 .^) 2 


ax* 


ay* 
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i ( aul + irl) 2 ] 

2 ay* ax* 


+ 
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equations of state : 
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r = p 

yM 2 


e* = 


T* 


y (y-1) M 2 
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The conservation equation source terms are given by. 


S c = ~ 


p y 
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p*u*v*_2 l dp* v* 1 I 1 * { 1 ^ v ' + fo— ) 111-8 

y * 3 Re ax* y* Re y* 3 dx* dy' 


p*v* a _ 2 1 dp’v* ^* 1 P* ( d v '* _ v '*) m-9 

s y* y » 3 Re dy * y* 3 y* dy* y* 




p*e*v* . 1 l JC* 3ri _p.Zl-.J_n* ( Zl) 2 
y * y-l Pr Be *f a y’ 3y* y* * e y* 
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3 pe H ax* dy* y* 


for the circular port region of the computational domain, 
slot itself, these terms are given by: 


S c = 2P P * 


r* 

Jb* 


For the 

III-ll 
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where 


_ 2 1 3p* dvT_ + 2 jH-1) 

3 Re dx * dz* b* * e dz* 


a = 2 1 3u* 2 

>■ 3 Re c(y* 3z* i>* * e 3** „ 


• Re * 1 az‘ 2 3 z- dz • 


3 Re p 3x* 5*y* dz* 
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III-13 


III-14 


A ~ 2 P p 2>» e « > 

, , i i 

Y~1 Pr Re Af 2 Jb* dz* r 


(grain burning) III-15 
(no grain burning ) 111 


Turbulence Model 

The low Reynolds number form of the two-equation, k-e model 
developed by Jones and Launder 1415 is used in the present analysis. 
This model employs two variables, the kinetic energy of turbulent 
velocity fluctuations , k, and the rate of dissipation of kinetic 
energy e, from which the turbulent viscosity, \i t , is calculated. 
The turbulent thermal conductivity, K t , is obtained by assuming a 
constant value of the turbulent Prandtl number, a t -0.91. The 
turbulent kinetic energy and energy dissipation rate are determined 
from the solution of the following differential equations: 


k-eguation : 

a P ^ap*u^ ap*v2fc = _1_{ JL[ — ] +Sri (n,* + -^) -irr ] } 

" + ~dv* Re'dx* 1 *' ax* dy* dy* 


Re 




Re' 


dk \ 
dx* 




+ Sz 


JJJ-17 
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e -equation: 


ap*e ^ ap*u*e^ap*v*e 
at* dx* dy * 


_!_{_!_[ (p -*-] +-^-[ (p,* + — 

He ax* ' «e ax* dy * ' ° 


+ 1 (7 6 P +C — P C p»jj- 2 + ( j^Hl) 

~Re Cl lc Pl Cl k P 2 Re 29 k Re p* ax* 2 dy* 2 


9 2 


ll ax-‘ dy- 3 ax-y- 


• a 


where 


1 3 1 3x* dy* dx* dy * dy* dx* 


and 




The source terms are given by 

o - p*v**. 1 afc.i 2 

* y* Re y* 1 °* dy* Re 3 * y* dx* dy* 


a p‘v*.fc 1 ± u •Xl(duL + ^L) 

3 * - o £ ) dy* Re 3 c 'k 9t y^ dx* dy* 


in the circular port region of the domain, and by 


S k 


= P* 



■r'P' 


S ] 
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III-22 


III-23 


III-24 
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l . /4 r /dtf*\ a 5nr* du* _ dvr* + ( d«l) ' + ( _?5C) '} 

- m*' , 3 t< al 7) az-ay- 1 ‘az- 1 az- 


dv 


3 k az* 


for the portion of the domain within the. slot. 

Heat Transfer Relations 

The goal of the present analysis is to examine the interaction 
between the igniter plume, the developing flow field within th 
h!ad-Jnd star grain slots, and the rate of flame spread over the 
arain surface This, in turn, may provide insight into the 
lopiSprTateness of a particular grain design or a particular 
igniter design for a given SRM. Many previous analyses utilize a 
convection heat transfer model, while others, sueh asJoh J ston ' 
t~ -i 1 n 7 p pl combined convection-radiation model. Tne prese 
analysis utilizes a simple convection-only model developed by Kays 
and Leung”, and used previously to correlate heat transfer within 
the O-ring gap of the Space Shuttle nozzle- to-case joint , given 

by 




0.152 Re 0 - 9 Pi 


0.833 [2. 25 ln(0 . 114 J?e 0,9 )+13.2 Pr-5.8] 
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where the Reynolds number is based on the hydraulic diameter of a 
single star grain slot. Fluid properties are based on a reference 
temperature" denoted as 19 


T_ = T + 0.5 (T wlJ -D + 0.22(3Va-T) 
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and the adiabatic wall temperature is calculated from 


T„ a = T + *JPl 


V 2 

2c_ 
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The propellant burning rate is assumed to be of the form 

x = r ref ( -£—) ”exp [ o p ( T p -T ref ) ] III-29 

Pref 

The constants appearing in Eq. (27) are defined in Table III-l. 
Erosive burning is assumed to be negligible during the very early 
portion of the ignition transient. 
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Table III-l. 


Constants used in burning rate law 


Constant 

Value 

Units 

*-*ref 

0.01078 

m sec' 1 

Pref 

6898.2 

KPa 

n 

0.35 

- 

<*p 

0.002 

K' 1 

^ref 

300 

K 


In order to determine when a given element of the solid propellant 
reaches a critical ignition temperature, one must know the surface 
temperature of the solid grain. This in turn depends on the 
amount of heat transferred to the grain from the hot gases 

The grain is considered to be a semi-infinite slab whose 
temperature is initially uniform. Heat tranfer to the slab is 
assumed to be one-dimensional. Thus 


dT P 3% III-30 

at p dz 2 


with the boundary conditions 


T p (t,co) = T pi 


ar p (t/Q) 

dz 


- ^{T-T p {t,0 )] 

K p 
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and the initial condition 


T p (0,Z ) = T pi 
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The assumption of one dimensional conduction heat transfer implies 
that ignition of adjacent grain surface elements is attributed to 
direct heat transfer from the hot gas only. This appears to be a 
reasonable assumption because of the low thermal conductivity of 
the solid propellant. Assumed physical properties of the 
propellant g?ain material are given in Table IH-2. 
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Table III-2 . Solid propellant properties 


Property 

Value 

Units 

Pp 

1758 

Kg m~ 3 

K P 

0.4605 

W irT 1 K' 1 

(C p )p 

1256 

J Kg' 1 K' 1 

Tp, critical 

850 

K 


The initial propellant temperature is assumed to be T pi = 298°K. 

The (slot) wall shear stress must be approximated in orderto 
provide closure for the governing equations. A velocity profile 
across the slot width (in the z-direction) is assumed for the 
velocity components u and v. Pai's 20 polynomial form of the 
velocity profile for a flat duct is employed. For example, the u 
component of velocity is 


u 

^uz 


= 1 


p-g ( 

P-1 



g-1 

p-1 


( 



2P 
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where il is the velocity at the symmetry plane of the duct (slot) , 
p = 11 .O 6 , and q = 16 20 . Although Eq. (III-33) applies specifically 
to fully developed flow, it is used here to account for wall shear 
stress effects even though the flow is time variant. It should be 
noted that frictional effects at the slot walls are accounted for 
only in the unignited portion of the propellant; they are 
neglected for the ignited portion because of the strong transverse 
velocity component due to mass injection from the burning 

propellant. , . . , . 

Finally, the transverse velocity component w(z) is assumed to 

be linear, with w = 0 at the slot plane of symmetry; w at the slot 
wa ll can be determined from the known burning rate. Hence, 


w. 


wall 


= M. 


dw 

dz 


- - 2 w 
b ^ 
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Numerical Technique 

The numerical solution of the equations of motion is obtained 
utilizing the explicit, time-dependent, predictor-corrector finite 
difference method developed by MacCormack. This method has been 
widely used for the numerical solution of a variety of fluid 
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dynamics problems, including those containing mixed 
subsonic-supersonic flow regions, as in the present investigation. 
MacCormack's predictor corrector technique is well described in the 
existing literature. 22 - 23 In the present investigation, an explicit, 
fourth order numerical dissipation. scheme has been introduced into 
the set of equations to damp numerical oscillations induced by the 
severe gradients in the flow field associated with the developing 
igniter flow. The fourth-order damping scheme introduced by Holst 
and modified by Berman 25 and Kuruvila 26 is employed. This scheme 
involves certain free adjustable parameters, C x and C y , usually 
referred to as damping or dissipation coefficients. It is 
recommended in the literature 22 - 23 that C x , C y be such that 

A uniform rectangular grid is employed in the portion of the 
domain that represents the head-end section, which includes the 
star slot segment plus a distance equal to the port diameter. For 

region a 92x63 grid has been used, as shown in Figure III 3 . 
A grid of unequal spacing in the x-direction is used in the 
extended portion of the domain from the head-end segment to the 
motor exit. The scheme of Cebeci and Smith, 27 in which grid spacing 
is increased by a fixed percentage from an initial point, is 
utilized. 

Initial and Boundary Conditions 

In time marching problems, the initial conditions should in no 
way affect the steady state results. In theory, any initial 
conditions can be chosen. Of course, for an arbitrary set of 
initial conditions, the transient has no physical meaning and only 
steady state condition is a meaningful representation of the 

flow field. , . _ 

Since in this study the starting transient is of primary 
importance, the initial conditions must reflect the actual values 
of the flow field variables at time t = 0. Therefore, both 

components of the velocity have been set equal to zero, u = v = 0, 
and ambient values have been chosen for pressure and temperature 
everywhere in the flow field. Density and internal energy are 
obtained from the equation of state. The turbulent kinetic energy, 
k, and its dissipation rate, £, are also zero initially. However, 
imposing this value would lead to a singularity in the k and £ 
equations; thus, a very small value has been assigned to these two 

variables. . _ .... , 

Boundary conditions must also be specified for all of the 
dependent variables, u, v, p, e, k, £, along with corresponding 
values of p and T. Of primary importance is the specification of 
the time variant conditions from the developing igniter flow at the 
inflow boundary of the calculation domain. The conditions 
specified are those from the single-port igniter used in the Space 
Shuttle SRM. The igniter mass flow vs. time trace is shown in 
Figure III-4, using data from Ref. 28. From the known igniter mass 
flow rate vs. time trace and the geometry of the igniter nozzle, 
the values of the flow variables u, p, p, T at the igniter nozzle 
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exit are calculated using one-dimensional nozzle flow theory; the 
velocity is assumed to be zero, v=0, and the internal energy 
is calculated from the known temperature. Note that the numerical 
method employed is a shock capturing technique, so that the shock 
waves associated with the igniter plume are embedded in the 
solution. Inflow boundary conditions for k and e must be also 
specified. Since no experimentally measured profiles are 

available, the turbulent kinetic energy is assumed to be equal to 
a small percentage of the f reestream kinetic energy , and the energy 
dissipation rate is taken from an empirical relation used, for 
instance, in Ref. 29. Along most solid walls of the computational 
domain the "no slip" condition is enforced, that is, u — v — 0, 
a Q_go due to the low Reynolds number form of the two-equation 
turbulence model, k = 0 and e = 0. For these walls the temperature 
boundary condition is specified by assuming an adiabatic wall. The 
exception is the cylindrical port portion of the domain downstream 
of the star slot section. For this region, the wall temperature is 
determined in the same manner as in the upstream star slot, and a 
"blowing wall" (v * 0) condition is imposed— to a ccount for mass 
addition due to burning. The top boundary of the computational 
domain represents the centerline of the motor, so that a symmetry 
boundary condition is enforced here. 

The outflow boundary at the aft end of the motor represents 
the final location where values of the dependent variables must be 
specified. The approach adopted is to place a diaphragm at the 
motor exit, immediately preceding a fictious "nozzle". When the 
pressure differential across the diaphragm reaches some nominal 
value, say 1 atm, the diaphragm bursts. The nozzle is assumed to 
fill instantaneously. Continuous checks are made on the ammount of 
flow approaching the exit and the amount of flow the "nozzle" could 
pass given the existing upstream conditions. A simple one 
dimensional time -dependent mass conservation calculation is then 
utilized to provide boundary conditions for the next time step. It 
should be noted that, for the first 200 msec or so following 
igniter firing, the downstream boundary condition is not as 
important as it would be later in the ignition transient. 


Sample Results 


Cold Flow „ „ . m 

After verification of the overall flow model and computational 

technique by comparison with known solutions or measurements for 
problems such as supersonic laminar flow over a flat plate 30 and 
developing turbulent flow in a pipe 31 , attention was turned to the 
igniter plume as a measure of the model's ability to predict more 
complex flow fields. . Experimental data obtained both by Conover 
and in the series of tests described herein were utilized. 
Schlieren photographs of the exhaust plume (s) were taken for each 
igniter at various values of upstream total pressure. In 
particular, the Space Shuttle SRM igniter is a single port, conical 
nozzle configuration with an area ratio of 1.428 and a half angle 
of 27.2 degrees, with a design exit Mach number of 1.79. Figure 
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HI-5 compares the calculated plume Mach number contours for an 
igniter total pressure of 102 atm (1500 psi) and a total 
temperature of 316 °K with the Schlieren photograph of the plume at 
the same conditions. A qualitive comparison indicates that the 
results are quite satisfactory, although the shock wave appears to 

be smeared over several grid points. 

Next, calculations were performed for the cold flow _ (no 
burning) flow field within the star slot of the head-end section. 
The computed results can be compared with the oil smear data 
obtained in the present cold flow tests. Both the case of a sing e 
port igniter and a multi-port (canted) igniter are considered The 
multi-port igniter exhaust jet is aligned with the slot, at an 
angle of 45 degrees with respect to the motor centerline. Typical 
results are shown in Figures 111-6:111-9. Figures III-6 and III-7 
compare results for an igniter pressure of 6 8 atm (100 psi , while 
Figures III -8 and III-9 compare results for 34 atm (500 psi) . The 
igniter (not shown) is located near the upper right corner of the 
slot, exhausting from right to left. 


Attention was next turned to the problem of heat transfer and 
the calculated progression of the burning surface of the propellant 
grain within the star slot. As time progresses, the propellant 
surface temperature rises due to heat transfer from the hot gas 
flowing over it. Propellant ignition is assumed to occur when the 
temperature of the surface reaches a critical ignition temperature 
(850° K) . Obviously, the resulting flame spread is dependent on 
the heat transfer correlation assumed (eg., Eq. III-26) . Figure 
HI-10 illustrates a typical predicted burn sequence for the single 
port igniter presently used on the Space Shuttle SRM. First 
ignition of the propellant surface occurs in the vicinity of the 
shock located in the igniter plume, as might be expected. 
Interestingly, this burning progression pattern can be anticipated 
from calculated cold flow heat transfer contours, as can be seen in 
Figure III-ll, which illustrates calculated Nusselt number contours 
for cold flow over a range of igniter pressures from 6.8 atm (100 
psi) to 102 atm (1500 psi) . Again, the igniter exhausts from right 


Figure III-12 compares the calculated head-end pressure-time 
trace with values obtained from measurements taken from actual 
motor firings 31 , as well as with values calculated from a typical 
P(x,t) model 10 . The comparison is made for the time interval 
0<t<120 milliseconds, over which the first propellant grain 
ignition occurs and flame spreading begins on the slot surface. It 
can be seen that the P(x,t) model significantly underpredicts the 
head end pressure during this interval, while the present model 
matches the measured pressure trace with reasonable accuracy. 
During this period, the flame spread is fairly slow, as would be 
expected. A flame first appears on the propellant surface at 25 
msec after the igniter firing, and at 120 msec approximately 20 
percent of the star slot grain is burning. The rapid pressure rise 
shown by the P(x,t) model at approximately 115 msec corresponds to 
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Figure III-5: Igniter Flowfield Comparison (cold-flow) 
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Oil Smear 



Calculated 

Figure III-6: Single Port Igniter Cold-Flow, P ignite 


b 6.8 atm 
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Calculated 

Figure III-7* 45 Degree Igniter Cold-Flow, P ifln i t . r “ 6.8 atm 
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Oil Smear 



Calculated 

Figure III-8: Single Port Igniter Cold-Flow, P lBnlt . r - 34 atm 
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Calculated (Cold Flow) Heat Transfer Contours 
Igniter Pressure - 6.8 atm (100 psl) 



Calculated (Cold Row) Heat Transfer Contours 
Igniter Pressure 3^ atm (500 psl) 



Calculated (Cold Row) Heat Transfer Contours 
Igniter Pressure « 102 atm (1500 psl) 



Figure hi- 11s Calculated Nu Contours (Cold-Flow) Single Port Igniter 
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Figure 111-12: (Early) Ignition Transient Head-End Pressure Rise 
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the ignition of the CP portion of the grain in that model. First 
ianition of the CP portion is predicted to be at approximately 80 
msec in the present model. It should be noted that the present 
model tends to under-predict the pressure rise in the head-end for 
t > 120 msec. This suggests that the predicted flame spread for 
time t > 12o' msec is too low. This may be due to one or more of 
several effects, including under-prediction of the convection heat 
transfer coefficient, the presence of significant radiative heat 
transfer, or the significance of 3-D conduction effects within the 
propellant grain. 
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IV. CAVENY PROGRAM MODIFICATIONS 


To incorporate the results obtained from the experimental 
and analytical work modifications were made to the original so 1 
rocket motor ignition code. These modifications proceeded alo g 
two independent paths: changing the code so that it could work 
with the interface program, and changing the code to improve its 
accuracy when used with star grain segments. 

To significantly increase the accuracy of the prediction 
codes used by the Caveny program, two of its fundamental 
aqsumDtions had to be discarded and replaced with other 
computational models. First, the original code assumed {*** g* e 
Drooellant grain ignited at a given temperature, that the entire 
seoment was not burning until it reached this temperature, an 
that after reaching that temperature the whole segment would burn 
evenly t Second? rather Chan calculating the burning perimeters of 
the segments based on input motor geometries, the Cavenyprogram 
accepted an equivalent effective perimeter of the segment at 
various stations along the grain. This effective perimeter was 
modeled as a circular perforated grain, and the gas dynamics 
equations assumed one dimensional flow xn a cylinder. These 
assumptions led to inaccuracies for star grams. 

To correct the inadequacies of the first assumption, a new 
model for the burning area was introduced. A system was 
introduced in which the fraction of the segment burning was a 
function of simulation time rather than local temperature. A 
table of percentage burning surface as a function of time c 
be specified for any of the segments. This is obtained from the 
CTD Sllel Inscribed in Section III. A variable NBPTAB (number of 
burning perimeter tables) was added to the NAME namelist The 
variables NBPENT, BPTIME, and BPFRAC were added to the INPU 
common block, respectively representing up to thirty burning 
perimeter entries in a given table and up to twenty values of 
time and burning surface fraction per station. These tables 
contain discrete values of burning surface fraction as a function 
o! “me. A linear interpolation function INTERP we. also added to 
the Caveny code to calculate intermediate values. In addition, a 
new model for turning on burning of the grain had to be 
introduced to the code as well. As previously stated, in the 
previous model, the grain segments were assumed to begin burning 
when they reached the ignition temperature. 

A subroutine to calculate burning areas as functions of 
distance burned was also added. The original version of the 
Caveny program held the burning perimeter constant during the 
course of a run. This assumption was considered a valid 
aDDroximation since the simulated burning times and burn 
distances would be small in comparison to the total burn time and 

overall size of the motor. 
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To correct the shortcomings of this second assumption, a new 
model for calculating the burning perimeter was introduced. In 
this model, a number of different grain geometries can be 
spUifiSl in the input case file. A variable NBGTAB (number of 
burning geometry tables) was added to the NAME namelist. The 
variables NBGENT and BGVAL were added to the INPUT common block, 
respectively representing the number of burning geometry 
variables in a given table and up to ten variables that have 
different meanings depending on the geometry model selected, 
subroutine to integrate the burning rate,, BIXTALC and a variable 
representing the total burn distance at an axial station, BDIST, 
were also added to the Caveny program. 

The Caveny program reads a single input file that specifies 
all parameters for a particular case . This file consists of at 
least one namelist, followed by a number of tables, followed by a 
number of optional namelists. The program reads the first 
namelist and from it modifies the default values set in the 
program. ’ This namelist includes variables describing the motor 
geometry, simulation parameters such as time step and end time, 
options for printing results, etc. Following this initial 
namelist are the tables for igniter mass flow and motor geometry. 
After these tables are a number of optional NAME namelists, each 
of which can be used to modify any of the simulation parameters 
listed in NAME. These subsequent namelists are read when the 
Caveny program has reached the end of a run. The new set of 
variables is generally restricted to output parameters such as 
time increment, print interval, and updated end time. 

The tables associated with the new burning fraction function 
are specified by a single integer value for NBPENT ( I ) (read using 
a 1110 format) between 1 and 20 representing the number of 
burning fraction versus time pairs for that location. These pairs 
of values for BPTIME ( I, J) and BPFRAC(I,J) (read using a 2F10.4 
format) immediately follow, one per line. An example set of 
burning fraction tables is given in Table IV. 1. 


In each of the following example listings, a vertical bar 
and all text following it should be interpreted as comments and 
should not be inserted into the input files. 
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Table IV. 1: Sample Burning 


Perimeter Fraction Table 


&NAME 

NBPTAB = 4 

SEND 

5 

0.0000 0.0000 

0.2500 0.1000 

0.5000 0.2000 

0.7500 0.4000 

1.0000 0.8000 

4 


0 <= n <= 30 tables listed below 

5 entries in Table #1 
0% burning at time = 0.00 sec 


80% burning at time - 1.00 sec 
4 entries in Table #2 


These tables must be placed in the input file immediately 
following the initial NAME namelist, but before any of the other 
tables Much as igniter mass flow, grain geometry, etc.) that can 
be specified The number of tables provided must correspond to 
the Ser specified in NBPTAB, and the number of entr res per 
table must correspond to the value provided for NBPENT(J) . 

The tables associated with the new burning geometry function ar 
specified in a single integer value for NBGENT(I) (read using a 
i?in fJraat) between 1 and 10 representing the number of 
iariab?STin the tXle for that location. These variables are 
interpreted differently by the various motor )j S i? h 
be sDecified. In all cases, the first variable BGVAL(J,1) 
indicates the type of geometry specified at the station. The 
n^e? of variables that follow and their meanings are determined 
bv this first value. Table IV. 2 shows the valid values and their 
meanings . Tables IV. 3 through IV. 7 detail the meanings of the 
remaining variables in the table for each type of available 
geometry. When possible, variable definitions are linked to 
correspond ing quantities between various grain configurations. 


Table IV. 2: Caveny Code Input Tables 


BGVAL (1,1) 

Grain type 

1 . 0 

Circular perforated 

2 . 0 

Slotted 

3.0 

Star 

4.0 

Wagonwheel 

5.0 

Dogbone 
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Table IV. 3: Circular Perforated Grain Variables 


Variable 

Units 

Variable Meaning 

BGVAL ( I , 1) 

none 

1.0 = Circular perforated grain 

BGVAL(I, 2) 

inches 

Interior diameter (Dl) 

BGVAL ( I / 3 ) 

inches 

Exterior diameter (D2) 


Table IV. 4: Slotted Grain Variables 


Variable 

Units 

Variable Meaning 

BGVAL (1,1) 

none 

2.0 = Slotted grain 

BGVAL (1,2) 

inches 

Interior diameter (Dl) 

BGVAL (1,3) 

inches 

Exterior diameter (D2) 

BGVAL (1,4) 

none 

Number of star points (N) 

BGVAL (1,5) 

inches 

Slot width 


Table IV. 5: Star Grain Variables 


Variable 

Units 

Variable Meaning 

BGVAL (1,1) 

none 


BGVAL (1,2) 

inches 

Interior diameter (Dl) 

BGVAL (1,3) 

inches 

Exterior diameter (D2) 

BGVAL (1,4) 

none 

Number of star points (N) 

BGVAL (1,5) 

degrees 

Star half-angle 
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Table IV. 6: Wagonwheel Grain Variables 


Variable 

Units 

Variable Meaninq 

RGVAL ( 1 . 1 ) 

none 

4.0 = Wagonwheel grain 

BGVAL (1,2) 

inches 

Interior diameter (Dl) 

BGVAL (1,3) 

inches 

Exterior diameter (D2) 

RGVAL (1,4) 

none 

Number of star points (N) 

D O V rVJJ ( a / t / 

BGVAL (1,5) 

degrees 

Star half-angle 


Table IV. 7: Dogbone Grain Variables 


Variable 

Units 

Variable Meaning 

BGVAL (1,1) 

none 

5.0= Dogbone grain 

BGVAL (1,2) 

inches 

Interior diameter (Dl) 

BGVAL (1,3) 

inches 

Exterior diameter (D2) 

BGVAL (1,4) 

none 

Number of star points (N) 

BGVAL (I, 5) 

degrees 

Star half-angle 


The routine to calculate the burning perimeters as functions 
of time were based on the work presented in Reference 2. A 
sample set of geometry description tables is given in Table IV. 8. 


&NAME 


Table IV. 8: 


NBGTAB = 2 


Sample Burning Geometry Table 


&END 

1.0 

54.345 

23.456 

2.0 

54.345 

123.456 

6.0 

15.0 


Table #1: CP Grain (3 vars) 

Inner diameter (inches) 
Outer diameter (inches) 
Table #2: Star Grain (5 vars) 

Inner diameter (inches) 
Outer diameter (inches) 
Number of star points (-) 
Star half-angle (degree) 
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The original Caveny program produced two types of output 
files- a 132-column file reflecting key simulation variables as 
functions of time and (where appropriate) space, and a plot file 
that was designed specifically to work with the Princeton 
University plotter . 

Rather than changing the format of the default output file 
of the Caveny program, new routines to generate a separate plot 
file were added instead. Of. all the simulation variables 
calculated and printed in the original program, only the 
variables listed in Table IV. 9 were considered to be of interest 

for plotting. 


Table IV. 9: Variables Saved in caveny. plot File 


Time and Space Dependent 
Variables 

Time De 

pendent Variables 

Name 

Meaning 

Name 

Meaning 

PSIA 

Pressure 

PESTAP 

Exit Pressure 

T 

Temperature 

BMTOT 

Mass Burned 

M 

Mach Number 

M 

Mach Number (Nozzle) 

BR 

Burn Rate 

FLBF 

Thrust 

TAUB 

Burn Distance 

- none - 

(unused) 

XMSA 

Mass Flow Rate 

- none - 

(unused) 

TPS 

Surface Temperature 

- none - 

(unused) 

REP 

Reynolds Number 

- none - 

(unused) 


To make the Caveny program more maintainable and modifiable, 
it was restructured to take advantage of the program development 
features offered by the Unix platform. In its original form, the 
program was approximately 3500 lines of FORTRAN code in a single 
file All common blocks and namelists had to be repeated in each 
subroutine that used them. These common blocks and namelists 
often had incompatible lengths and variable names substituted 
between different subroutines. 

To make the code more readable and more easily edited, the 
source was broken up into separate files that contained only one 
subroutine each. The filenames had a basename that was the same 
as the subroutine they contained and a .f suf f: ix. In addition, 
all common blocks and namelists were also placed in individual 
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files. These files were given names corresponding to the name of 
the common block or namelist and a .com or .nam suffix 
respectively. These files were then referenced m the code 
throuqh INCLUDE statements. This organization insured that , 
identical copies of the common blocks and namelists were being 
used by all subroutines. 

This new organization offered another advantage in program 
development. With each subroutine placed m only one file, only 
the files that had been changed since the previous build had to 
be recompiled, rather than compiling all modules after any 
change. This change was further aided through use of the Unix 
makeutility. An appropriate makefile that was suitable for the 
whole project was created. This program organization made 
modifying the program much easier. A complete listing of this 
makefile is provided in Appendix B. 
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V. CAVENY PROGRAM INTERFACE 


The interface application Is designed to make running the 
Caveny program easier for potential users, both novice and 
expert. The three natural divisions of using the program were: 
building the input data files, executing the program, and 
collecting and interpreting the results it generated. The 
principal parts of the interface program, are designed to 
correspond to these same three areas. 

An input file editor which is included is more than a simple 
text processor. Ideally, it should prevent the user from creating 
an input file that the program could not run. In practice, this 
would be nearly impossible, but it could ensure that all 
variables fall within certain constraints. From within the 
interface application, the Caveny program runs in the background. - 
The interface program automatically reads in the data generated 
from the most recent execution. It then be able to produces a 
number of plots simultaneously. These plots are dynamic, in that 
the user configures the variables displayed, the range of values 
over which they are shown, and all aspects of the plots' 
appearance . 

The parts of the interface program that generate input files 
for the Caveny program are the most crucial in making the program 
easy to use. Configuring input files proved to be the most 
da unting and error-prone aspect of running the Caveny program. As 
such, attention was paid to providing the user with a number of 
features to make editing these files more convenient. 

The input file editor is designed to provide two critical 
functions: on-line documentation of all simulation variables and 
error- checking of the values provided for those variables. The 
documentation for the variables includes (but is not necessarily 
limited to) definitions and descriptive text, units in which the 
variables should be input, and default values for each. Limited 
error-checking is provided by comparing the input values against 
known valid minimum and maximum values. While this technique does 
not insure that a case is consistent and will run properly, it 
does provide a safeguard against obvious errors. 

A simple text file format has been chosen to support this 
part of the program. A file containing all data about the 
variables that could be modified in the Caveny program input case 
files would have significant advantages over compiling this 
information into the interface program. Descriptive text and 
default values could be changed by a user without recompiling the 
code, and information tailored to different audiences could be 
maintained. More importantly, if new variables are added to the 
Caveny program in the future, these variables can be easily 
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incorporated into the input file generation part of the interface 
program. 


The output capability of the interface program is provided 
to reduce the amount of data generated by the Caveny program and 
to present the data in a more easily recognizable and 
understandable fashion than simple printed data. 


The ability to generate plots of the data presents a 
significant problem. First, the data is not organized in a 
fashion that would be well-suited to work with an existing 
plotting package. Although the output file can be changed to 
contain only numbers, there are no existing programs that could 
accomodate the volume of data produced by the Caveny program 
One of the main considerations has to he reducing the amount 
data provided to the user. The most efficient way to do this is 
to create a number of plots that allows. the ^er to ^ of 
dat-a at one time. These plots are dynamic, and the variables tney 
display and the ranges over which the data are J 1 ® p l^® d b ® a “ h ^ e 
changed. More then one plot can be maintained at a time by the 

program. 

An integrated environment from which the user does not need 
to exit offers the most significant advantages m ease of use. 

The goal is to provide a shell program that encourages use by 
Sosl not lamiliar with the Caveny program or the Unrx operating 
system. To this end, a number of features that were not critica 
to using the program were incorporated into the interface. A 
simple text editor based on the xedit sample application has been 
designed and included. 
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VI. XPLOT PROGRAM DESCRIPTION 


The specifications set forth in Chapter V led to the 
development of a stand-alone application named xplot, a graphical 
user interface for the Caveny program. 

A number of issues had to be resolved at the outset, since 
they would have significant impact on program development. The 
computers available in the College of Engineering at the time the 
oroiect began were the IBM mainframe, Sun workstations, and a 
number of IBM PC compatibles. The mainframe offered no graphig? 
capability. The PCs did not have the memory, storage,. or 
processing capability necessary to run even the original Caveny 
proaram. The Sun workstations offered the power and graphic 
output capability needed. Since the operating system used on 
these machines was Unix, that would be the operating system for 
which xplot would be written. 

The X Window System from the X Consortium at MIT was . chosen 
as the graphics package for the interface program. Competitive 
packages available for Unix workstations that were considered 
include- Pixrect (Sun Microsystems' device-dependent screen 
drawing mechanism), PHIGS (a Programmable, Hierarchical 
Interactive Graphics System - designed for creating, 
manipulating, and rendering objects in both two- and 
three-dimensional space) , and GL (a multi-vendor general-purpose 
drawing environment from Silicon Graphics) . X is more appropna e 
for the specific application for a number of important reasons. X 
is widely available on most Unix platforms that support graphics 
output. The calls that are made to X routines are 

platform- independent ; as such, programs written on a machine that 
supports X should only require recompilation to work on a 
different machine. Revisions of X have been available at no 
charge for both internal use and worldwide distribution. Most 
important, however, X is the only system available that provides 
any support for developing full-scale, interactive user 
applications; most of the other systems only offer libraries of 
drawing calls. 

The choice of Unix for the operating system and X for the 
graphics system made C the most natural programming language to 
use in writing xplot. At the project's inception, C compilers 
were provided for use with the Sun workstations. More important, 
however, was the need for a language offering both data 
structures and dynamic memory allocation to use in conjunction 
with subroutine calls to X. X routines require data to be passed 
and returned using data structures, and they often require the 
calling application to both allocate and free memory. Since 
FORTRAN does not support either capability, C was the obvious 

choice . 
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From the outset of the program, xplot and the Caveny program 
were designed to work together but remain two separate processes. 
The Unix environment supports true preemptive multitasking, so it 
was possible to create an interface shell that executes the . 
underlying solid rocket motor prediction code. From the user s 
standpoint, there is no need to exit from xplot. All the details 
of building the input file, running the Caveny program, and 
loading the resulting data were made transparent. 

There are several reasons that support making the interface 
nrortram a seoarate process. There is no reason why the user 
should have to wait for the Caveny program to finish executing 
hpfore usinq other parts of xplot. Had the programs been 
integrated, there would have been no way for the user to continue 
building input datasets are creating plots from previous runs . 

The user would have had to wait for the simulation to finish 
before continuing the analysis. It would also have added to the 
complexity of the two programs if C and FORTRAN routines had been 
used together. The way in which arguments are passed to 
subroutines differs between the two languages, and special 
libraries that support both languages must also be used. Finally, 
theCaveny program could be left unmodified if two processes were 
used; any Y other implementation would have required making changes 
to the existing code. 

The Unix operating system provides true preemptive 
multitasking for multiple processes on the same machine. Any 
number of executables can be running concurrently, and 
time-sharing between them is handled by the operating sy . 
Associated with each instance of an executable is the environment 
under which it is run. Collectively, the executable code, its 
environmental variables, and the memory it has reserved for its 
internal data are referred to as a process. 

Under the Unix operating system, there is only one way to 
start a new process. The fork system command creates a duplicate 
of the process that is currently executing, including all of its 
environmental variables and data. The fork command is unique in 
that it has a single entry point and two return points, one in 
the calling process and one in the process it creates. These two 
processes are identical, but can distinguish themselves from one 
another through the value returned from fork. The parent process 
is returned the process ID of the child process if the call 
succeeds or -1 if the call fails, and the child process is 
returned a value of zero . 

The fork command only creates a duplicate of the parent 
process- it does not run the desired executable. For this to 
occur , one of the exec family of system routines must be used to 
change the executable associated with the child P ro ^ s s- Note 
that all environment variables remain the same across a call to 
one of the exec functions. 
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The fork and exec functions are almost always used with one 
another. Typically, all a child process will do upon returning 
from a fork is transform itself into the appropriate executable 
using one of the exec calls. Assuming that the child process is 
performing some task for the parent process, the parent will 
often block until the child has finished executing. Under Unix, 
the wait command provides this facility. 

The fork-exec-wait combination is used in conjunction with 
xolot and the Caveny code. When the user chooses to execute the 
Caveny code after building an input case, the interface program 
forks to create a clone process. The child process then _ •. ^ • 

immediately transforms itself to the executable name stipulated 
in a corresponding field in the execution dialog. This executable 
is named caveny by default, xplot does not wait until the child 
process has finished executing to restore control to the user. 
Note, however, that the current dataset, input files, and plots 
will not necessarily be updated unless specifically requested y 
the user. Although similar results could have been achieved using 
the Unix system command, system is, in general, less reliable and 
requires much more overhead than a fork-exec-wait combination 
does . 

No direct communication takes place between xplot and the 
Caveny program. Instead, xplot is used only as a pre an 
post-processing application. Input files can be generated using 
xplot 's built-in editor. When the user selects to run the Caveny 
program, a file name must be specified in the execution dialog, 
xplot renames this file to the name of the input file expected by 
the Caveny code, caveny. in. The Caveny code is then run, and it 
generates two output files named caveny. out and caveny. plot 
Either file can be viewed using the file editor dialog, and the 
plot file can then be loaded as the current dataset and used to 
create new plots . The two programs communicate through these 
three files only. 

The organization of data generated by the Caveny program was 
reflected in xplot, both in the features it offered to the user 
and in the way it stored and manipulated the data internally, n 
general, the Caveny program produces two types of data: variables 
that are functions of both distance along the motor and time 
elapsed in the simulation, and variables that are functions of 
time alone. The model chosen for the interface had to be flexible 
enough to allow for both of these types of data. Instead of 
modifying the output file of the Caveny program, write statements 
were simply added at appropriate points m the code to produce a 
plot file with the desired format (see 2.5). 

Data were stored in a three-dimensional floating point array 
named plotData. The X-dimension was used as discrete locations 
along the grain, so this dimension was limited to 30 values The 
Y-dimension was used for the different variables output to the 
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clot file. Since the maximum number of variables that could be 
output was limited to 8 that were functions of both space _ and 
time and 8 that were functions of time alone, this dimension was 
limited to 16 values. The Z-dimension was used as an index into 
the time at which the variables had been written. Since this 
number was dependent on a number of factors U-e ; , the total, run 
duration step size, and print interval of the simulation) , the Z 
coordinate had to have the greatest dimension. Ultimately, this 
value was limited to 1024 points, a value that was generally 
appropriate for runs of less than one second. Figures VI . 1 a 
VI . 2 graphically show this data organization. 


Distanc* 


Tim* = 0.00 00 — oonda 
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Figure VI. 1: Two-Dimensional Data Organization 
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Figure VI. 2: One -Dimensional Data Organization 

In Figure VI. 1, the two-dimensional sets of data (used for 
quantities such as pressure, temperature, and Mach number) have 
values (Y axis) for all pairs of space (X) and time (Z -The 
increasing values of time are shown as planes of (X, Y) pairs, 
which is consistent with the way the data is written to the file. 
In Figure VI-2, one -dimensional^ data; (such as exit pressure, mass 
burned, and thrust) is shown. Here, data is a function of time 
(Z) alone. 

This storage technique made it possible to change the way in 
which two-dimensional data was presented in plots . The curves 
could be drawn at various points in time, with the variable shown 
as a function of distance along the grain. However, the data 
"cube" could also be rotated, so the data could be drawn at 
various "slices" of distance, with each variable shown as a 
function of simulation time. Example plots that reflect this 
capability are given in section VII of this report. 

The X Window System is a set of library subroutines that 
provide graphics and interface capabilities on Unix platforms. X 
is based on the client-server model; client applications make 
requests to an X server to perform certain actions, and the X 
server notifies the client about particular events that are of 
interest. X works equally well on a single computer or across a 
network. Each host computer can have an X server running on it, 
so it is possible for applications to make requests on both local 
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and remote hosts. X provides a platform-independent set of 
subroutines that client applications can use. All 
machine -dependent drawing and event scheduling are handled 
through the X server, an implementation tailored specifically 
the hardware on which it is running. 

x actually provides three different levels of possible 
interaction with client applications. At the lowest level is e 
£ f?£arv (Xlib) a library that provides extensive routines for 
drawing graphics' primitives (lines, rectangles, polygons, arcs, 
text etc ) and basic windowing environment routines. Xlib alone 
is g^eraily used only for basic drawing. and as the foundation 

fs an Sample of one ouch mechanism. It provides the framework 

f£ea^tnd n L^ xt 

is C an < ^event-driven tl system? 1 in^that n the rl 'client S application creates 

all polling and dispatching of events. The third leve 
interaction is a widget set, which greatly affects the 
functional i ty of the user interface. Short of writing custom 
widcrets thethings that can be done in a program are limited by 
the functionality of the widgets employed. This three- lere 
organization of interface code is shown m Figure VI . 3 . 



Figure VI. 3: X Toolkit Application Design Model 
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The client application can interact with any or all levels 
of interface provided: it can change the value, appearance, or 
behavior of one or more widgets through the widget set, control 
the overall application through Xt, issue drawing requests in a 
window through Xlib, or perform any operating-system specific 
calls concurrently. 

The xplot interface program uses all three levels of X. The 
widget set chosen was the Athena widgets provided as examples 
Y^lth the distribution from the X Consortium. Although Athena 
widgets do not have as nice an appearance or quite the 
functionality of commercial widget sets, they are adequate for 
the somewhat limited needs of the program. Further discussion of 
the widgets used in xplot and Xt in general is given below. 

As previously stated, the X environment is event-driven. 
Rather than performing certain actions in a specified order, as 
is customary in traditional programs, X applications must respond 
to events, which are usually user inputs to the program, either 
directly or indirectly. For example, at the Xlib level, a plot 
must be redrawn not only when created but also whenever it is 
exposed (from underneath other windows) or resized. All of the 
techniques discussed for drawing plots really refers to the 
actions taken whenever one of them must be updated . 


A simple layout was chosen for the plots generated by xplot. 
Only simple, rectangular 2D scatter or line plots are provided. 
Many simplifying assumptions about the kind of data being 
displayed can be made since the application is designed 
specifically to work with the output of the Caveny program. The 
independent variables are either axial distance along the grain 
or simulation time. The number of variable data traces displayed 
are limited to sixteen for time -dependent and eight for 
distance-dependent plots. The general layout characteristic of 
these plots is shown in Figure VI. 4. 
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Figure VI. 4: Plot Layout Design 

The size of a plot depends on both the size of the window 
enclosing it and the borders chosen for it. Only the data curves 
are drawn in the data area; all other plot components are drawn 
in the borders around it. The width and height of the plot then 
can be expressed as: 

Wplot = Wwindow - Bleft - Bright 
Hplot = Hwindow - Btop - Bbottom 

The elements of the plot are drawn in the following order: 
plot title, plot legend, x-axis title, axis labels, axes, 
gridlines, and then the plot data. This order was chosen to 
insure that various elements of the plot neither obscure nor 
erase the others . 

The plot title is centered over the main plot area. The font 
used for the title is Helvetica 18 bold. The width of the text 
string is computed using the XStringWidth function. Once the 



width of the text string is known, the x-component of the title 
can be calculated using the following equation: 

Xtitle = Bleft + (Wplot - Wtitle) / 2 


The plot legend is placed across the bottom of the plot 
area. Each element of the legend is drawn as a small rectangle 
having the same color as the data it represents, followed by a 
string label corresponding to the same data. Up to sixteen sets 
of data can be shown in the plot legend area, with each line 
containing as many as four identifying markers. 

The plot labels are drawn using successive calls to 
XDrawString . The sprintf function from the C standard library was 
used to print the text to a string in the format chosen. These 
axis labels are drawn at regular intervals along the appropriate 

axis . 

An axis title is drawn for the x-axis only, since X does not 
support drawing text rotated to arbitrary angles. The x-axis 
title consists of the independent variable, either axial location 
or time, drawn in Helvetica 10 Bold, centered below the main plot 
area. It also indicats which slice of either time or space is 
currently being displayed. 

The plot axes are drawn using a 1 point black line along the 
bottom and left edges of the plot area. Since all of the data 
that could be plotted, including dependent and independent data, 
contained only positive values, no provision was made for drawing 
the lines where Y = 0 or X = 0 inside the plot area. 

The grids on the plot are drawn at regular intervals 
corresponding to the number of axis labels. These lines are drawn 
as 0 width lines inside the bounding rectangle that formed the 
main plot area. 

The plot data is the most complicated aspect of the plot to 
draw. Each trace of data has its own attributes, and there can be 
up to 16 traces of data on a single plot. The location of any 
single data point can be calculated from interpolation: 

Xdata [ i ] = Bleft + (X [i] - Xmin) /Xrange * Wplot 

Ydata [ i] = Bbottom - (Y[i] - Ymin) /Yrange * Hplot 

Using these equations, data can lie outside the plot area 
rectangle if the scales have been manually set. Such values of 
dataare either less than the minimum or greater than the maximum 
values of the scales. To insure that the data drawn in the mam 
data area of the plot frame does not exceed the limits of the 
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area, the calls for drawing data were clipped to the plot data 
area The clip mask of the current graphics context is set to the 
rectangle R(Bleft, Btop, Wplot, Hplot) . Explicitly setting the 
clip mask overrides the mask that is set by an expose event. 
Rather then have the program maintain a list of all areas _ 
exposed, combine them with the plot area into a single clip mask, 
and then redraw the plot, xplot simply redraws the entire p ot on 
any expose event . 

Xlito provides line attributes as part of the graphics^ 
context, or current drawing environment . These attributes include 
foreground and background drawing colors, line width, style 
(solid, on-off dashed, double-dashed), join mode (rounded, 
mitered, or beveled), and end cap style (butt or rounded). _ ^ 

Implementing the various types of lines that can be specified in 
xplot only requires calling XSetLineAttributes with the 
appropriate values before any drawing commands . 

Since X is implemented using client-server methodology, _ the 
speed with which it can perform drawing operations is as limited 
by the rate it can receive commands as it is by the speed of the 
processor. To reduce the amount of communication (often network) 
t ra ffi C/ xlib provides routines for drawing arrays of objects as 
well as the routines used to draw single objects. Points, lines, 
line segments, rectangles, and arcs, can all be drawn either 
individually or in groups. To maximize refresh rate, xplot 
assembles a table of all the points for a given trace of data, 
then calls XDrawLines. 

Descriptions of the Athena widgets used in creating the 
xplot application are given in Table VI-1. 
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Table VI-1: Widgets Used in xplot 

Class 

Description 

Use 

' l 

Box 

Geometry management of arbitrary 
widgets in a box of specified 
dimension. 


Command 

A rectangular button used to invoke 


Dialog 

A window, used to hold other 
widgets, that can be placed and 
iconified by the window manager. 

All pop-up selection panels are 
instanced from dialogs. 

Form 

Organizes the layout and size of its 
child widgets. 

Organizing sub-objects in all 
dialogs. 

Label 

A rectangular area containing 
descriptive text. 

All dialogs, to indicate the 
function of an area or an 
associated text field or button. 


An array of text entries, drawn in a 
column, from which one or more can 
be selected. 

Scrolling lists, such as the 
command and history summary and 
the input file variables. 

MenuButto 

A rectangular button from which a 
pull-down menu can be displayed. 

Only used in the Top Level menu 
to contain all menus. 

Text 

An area containing text that can be 
either display-only or editable. 

All dialogs, whenever the user 
is prompted for data entry; in 
dialogs to display large amounts 
of read-only text. 

Toggle 

A button that can be in one of two 
states, either on or off. 

Setting options about the 
appearance of plots. 

Viewport 

An area with scroll bars that allows 
the user to view an object larger 
than the area provided. 

Scrolling lists, such as the 
command and history summary and 
the input file variables. 


The dialogs used in xplot are all constructed using similar 
techniques. A shell widget interacts with the window manager , 
allows a window to be moved, resized, restacked, or iconified, 
and can contain only a single child widget. A dialog shell widge 
wSs SJed for all of the dialogs other than the Top Level A form 
widget was created as the shell’s only child, and 1 an es 
geometry management of all the child widgets. Normally, these 
children were arranged by function into related groups. 

The way the Format dialog is constructed is representative 
of the dialogs used in xplot. A diagram of this dialog and 
^ c rrivP'n as Ficfure VI -5 . A sincfl^ Form widcjGt is the 

only°child of a dialog shell. A number of box widgets are its 
children, arranged in an order that logically divides functions 
of the dialogs. These boxes are attached to one another using 
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the geometry management facilities provided by the form widget, 
so the dialog is only as large as it needs to be to contain all 
of its children. The text fields and command buttons inside the 
boxes are organized in columns to minimize the space they occupy. 

Only a fraction of the widgets created for an Xt application 
need to be modified during program execution. The values of some 
of these widgets will be modified by the user, and some will be 
changed by the application to reflect the current state of the 
data it maintains. In general, widgets that need to be modified 
include text fields (used for data entry and presentation) , 
command buttons (used to reflect currently available or 
acceptable choices) , and toggle buttons (used to reflect a 
current state of a selection) . 


The entire diatog is 
contained wthin* 
DfatogShst widget 



Boms ere used to 
offset function* 
groups of widgets 
Som one enother 


Commands ere 
used to perform 
ecborm wthin too 
detog 


The only ditto of the 
Shot is a Form 
widget menegtoget 
other children 


Figure VI . 5 : Example dialog design 
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In an Xt application, the ID associated with a particular 
widget is returned by the function used to create it. For most of 
the widgets used in xplot, that ID is retained only as long as 
needed For instance, all widget creation routines require the 
widget's parent to be passed as one of the arguments However 
for decorative and geometry-management widgets (like form, , 

or label), once their children have been created, there is no 
reason to maintain their ID, since it is not necessary to modify 
them. For the widgets that did need to have their values read and 
set, the variable-arguments form of the Xt get and set functions, 
XtVaGetValues and XtVaSetValues , are used. 

Xt provides higher-level event processing than is possible 
through Xlib alone. Instead of notifying the client application 
that the user has pressed and released one of the mouse buttons 
at a certain location on the screen, it might instead notify that 
a button has been selected, a menu has popped up, or a scroll bar 
had been paged. Xt provides these higher-level events through a 
mechanism it defines as callback routines. Callback routines tel 
a widget how to respond to certain well-defined events. Most 
events will not be of interest to the application, and the 
widget's default behavior will be the only action that is either 
necessary or appropriate. 

In addition to the graphical interface provided by xplot, a 
command-line language was written. The grammar supported by xplot 
is very small. Since there are less than fifty possible legal 
tokens in the vocabulary, there is no need for a separate lexical 
analyzer or parser (such as those provided by the Unix utility 
programs lex and yacc) . Instead, each line is simply broken down 
into its constituent components by a single function call. This 
function uses string comparison routines to arrive at allowable 
input. To keep the code as straightforward as possible, the 
actions and options that can be set from the command-line are 
mapped to equivalent functions already present from within the 
graphical interface. 

The source code for xplot is organized into files, grouping 
subroutines that deal with related aspects of the application. 
This organization is used for two main reasons. First, to 
encapsulate the data and routines into a similar area, so that as 
little global data and global subroutine calls are necessary. 
Second, to make compilation both more efficient and quicker 
during development. The source code files correspond to the 
various main aspects of the program, i.e. the command line 
interface, the internal data, the plot formats, drawing the 

plots, etc. 

A complete listing of the source files used to build xplot 
and their associated functions is given xn Tables VI 2 a F* d VI 3 * 
All source files have a .c suffix, and all header files have a .h 

suffix . 
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Table VI-2: xplot Source Files 


Filename 

Lines 

Description 

xp commands . c 

248 

Passes input from the command-line, translates 
commands into their graphic equivalents. 


250 

Manipulates all of the internal data. 

vnfnrmal'S . c 

358 

Loads and saves formats from files. 

xpgeneral . c 

137 

Creates graphic contexts, loads colors and 
fonts- — 

xpgeometry . c 

169 

Responds to all conf igure-notif y anc * expose 
events generated bv the window manager. 


716 

Loads all input datasets. ; 



xplot . c 

69 

Initializes all of the Xt interface, calls all 
the subroutines to create all of the widgets, 
then calls XtAppMainLoop. 

xpmenus . c 

582 

Creates menus in the Top-Level dialog and 
responds to commands issued from these menus. 

■HR H 

731 

Draws all of the aspects of the plots. 

— 

xpwidqets . c 

907 

Creates all widgets in the application. 


Table VI-3: xplot Header Files 


Filename 

Lines 

xpcommands . h 

13 

xpdata . h 

25 

xpformats .h 

13 

xpgeneral . h 

23 

xpgeometry. h 

8 

xpinput .h 

9 

xpmenus . h 

14 

xpplots . h 

23 


47 
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VII. XPLOT USER INSTRUCTIONS 


The xplot program is designed to be usable both by those 
familiar with graphical interfaces as 

comfortable with command-lme interfaces. Although it is 
graphical in design, xplot does provide rudimentary support for a 
Sand-line language. A discussion of the graphical interface 
and the commands are presented first, then each of the text 
commands are discussed in terms of their graphic equivalents. 

The top-level dialog of xplot, the window that is displayed 
when the program is first run, contains controls to access all 
capabilities of the program. This dialog is shown in Figure 

VII. 1, 



A menu bar, located across the top of the dialog, provides 
menus for performing certain functions. To the immediate right o 
the dialog is an area for entering the pathnames associated with 
certain files recognizable by xplot, such as datasets, plot 
format files, etc. Below the menu bar are three areas used in 
conjunction with the command-line interface. To the left is a 
scrolling list containing all command verbs that xplot 
recognizes. To the right is a scrolling list containing the last 
tittle n commands entered by the user. Below these two scrolling 
lists is the text field in which all commands must be entered. 
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Selecting either the verbs from the command list or previous 
commands^from the history list will place them 

t-Mt field but it will not execute the command. To execute a 
command? the user must press either return or enter in the text 

fiel<3 The top-level menus -File", -Input-, -Execute", "Dataset-, 
"Format" "Plot", and "Windows" represent the logical divisions 
of the program. The "File" menu allows the user to select files 
for viewing, save the current file, or quit the application. The 
n T nnut 11 "Dataset”, and "Format" menus are identical in layout, 
breach de^f with a different kind of file that can be read by 

xplot . 

The InDUt File dialog allows the user to create and modify 
input files for the Caveny program without actually editing t e 
text files it uses for input. This dialog is shown in 
Figure VI I. 2. 



Fllenane 


Nescafe* 


Cnee and* j Start \ Append | Prev 1 Next [ Write 


Figure VII. 2: xplot Input File Editor Dialog 
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On the right side of the dialog a scrolling list denoted 
variable List Shows all o£ the case description variables and 
their current values. If the user selects one . , , 

elements the four fields of the dialog described below will be 
updated to reflect the properties associated with this variable. 

Once the user has selected a variable to edit, either 

crrn i i ina list or with one of the commands located 

through the scrolling list or w four text fields wi n be 

tlTlll in S wiSh ?he current data' about that variable. The FORTRAN 
name of the input parameter will be displayed in the Variable 
fSld A brief (one-line) description of the < parameter a 
significance to the Caveny program will be visible in the 
Definition field. The Explanation section will Jj°ntain up to six 

1±n f eS iil^iU^siS rcoisLf^f the in?erS“ representation 
variable (i e INTEGER, REAL, etc.), the FORTRAN format 
!n whic^the^ variable will be read (i.e. , F10.4) , .the units that 
it should be expressed in terms of, limits, on the JalMI, it canr, 
take on, etc. All three of these sections are read-only. 

The contents of the Entry section can vary greatly depending 
on the var?ab?e Selected. For simple scalar quantities^ which are 
summarized in the Variable List, a single value will be 
displayed . ^ 

The Filename field is a read-write text field that contains 
the pathname to be used in all read/write operations. When the 
dialog is first raised, this field will be empty. At this point, 
the user can either start a new input case file or enter 
pathname in this field. The value provided in this field will be 
used for all read and write operations. As such it is possible 
for a user to load in a new input dataset W1 ^ 0 ”t s aving the or 
changes made, save the new values m place of the old file. 

Save ?he current values to a new filename. The messages text 
field will contain pertinent information about the operations 
that the dialog performs and diagnostics about errors that it may 
encounter in attempting to carry out the user commands. 

Th ® d°™ a ?fble V VIi a i le F0? more^nformation'abou^the format 
=”caveiy p?og;aTinput F me°and the variables defined in 
namelist NAME, see Section IV. 
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Table VII. 1: Input File Editor Dialog Commands 


iaoit 

Command 

Action 

Start 

Returns to the first variable available and 
selects the first variable name in the scrolling 
list. 

Prev 

Next 

Advances to the previous/next variable available 
and highlights the previous/next variable 
displayed in the scrolling list. When either of 
these commands are selected, the value in the 
Entry field will be substituted for the previous 
value of the current variable, and this value 
will be reflected in the Variable List. 

Write 

Saves the current values to the filename 

Cancel 

Dismisses this dialog, but maintains all of the 
current settings. 


The Execution dialog contains text fields that allow a user 
to configure the specifics of running the Caveny program from 
within xplot. This dialog is shown m Figure VII. 3. 


^ % 

Execution Summary 


Parameters 

Values 

Connands 

Path 


Cancel 

Command 


Rccept 

Input File 


Execute 

r Output File 

U. - — — 



— 1 J — .M T-S 

a 1 Oirr 


Figure VII. 3: xplot Execution Setup Dialog 
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The Execution dialog consists of four text fields. The Path 
field displays the pathname, either global or local, that is t 
directory that will be searched for the Caveny program. The 
Command field displays the name of the Caveny program executable . 
By default, the path is simply the local path ./ and "J® 
executable is named "caveny". The Input File and Output File 
fields name, respectively, the files expected b Y the Caveny 
program as input and the filename that will be used for its 
output. By changing these filenames, various input cases and 
output datasets can be maintained by the. user. The user can also 
run the Caveny program from this dialog, once he has sele 
proper configuration, using the Execute button. T ^ e 
available with the Execution dialog are summarized in Table 14. 


Table VII. 2: Execution Dialog Commands 


Command 

laDlC Vii.4. ? . — — — 

Action 

Cancel 

Dismisses this dialog, saving all of its current 
settinqs. 

Accept 

Dismisses this dialog, saving all of its current 
settinqs to internal values for these parameters. 

Execute 

Immediately runs the Caveny program using the 
parameters currently specified in the dialog; 
saves all current values. . — 


The Plot Manager dialog contains controls that 
user to configure up to eight plots simultaneously, 
is shown in Figure VII. 4. 


allow the 
This dialog 
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Figure VII. 4: Plot Manager Dialog 


The Commands column allows the user to switch the 
independent variable, move either along the grain or m time, or 
animate any of the eight available plots The commands available 
in the Plot Manager dialog are summarized in Table VII.J. 


Table VII. 3: Plot Manager Dialog Commands 


J 

Command 

Action 

Space 

Plot/ 

Time Plot 

Changes the orientation of the currently active plot. 

Toggles . 

between the plot being displayed as variables as 

Increment 

Decrement 

Changes the view to the next/previous time or spatial 
plane location. 

Animate » 
Animate << 

Quickly increments/decrements the view to the 
beqinninq/end of the time or spatial stations. 

Edit 

Format 

Brings up the Format dialog, with the values shown 
corresponding to the currently active plot. If the 
format dialog is already visible, then change the dialog 
so that it displays the values associated with the 
currently active plot. - 

Cancel 

Dismisses this dialog, saving all of its current 
settings . 
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The Format Dialog allows the user to edit the format used to 
display any of the eight plots available. This dialog is shown in 

Figure VII. 5. 
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Figure VII. 5: xplot Format Dialog 
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The format dialog will allow the user to edit the formats of 
plots other than the active plot, and these plots need not be 
either created or visible to edit their format. These formats can 
be both loaded from and saved to external files for future use. 
This dialog works in conjunction with the Data Attributes dialog 
to modify the internal state of the formats of the plots. 

Formats can be saved or loaded just like datasets. When the 
program is begun, a default set of format values will be 
installed for each available plot. The user can edit the format 
to the style preferred, then save and recall these settings using 
the Filename field. The Title field specifies the label that will 
appear centered across the top of the plot. 

There are two columns of toggle buttons labeled Variables. 
The left column contains all of the variables in the current 
dataset that are functions of both position and time. The 
variables listed in the right column are functions of time alone. 
Note that nothing will be displayed on the plot if the user 
selects only time-dependent data from the Format dialog and 
attempts to create a space plot. The organization of data m a 
p]_ot file is discussed in Section VI. 

By selecting one of the toggle buttons in the Options 
column, the user can determine which features to use on the 
current plot. The first option, "Scale Data", forces the Y axis 
to take on minimum and maximum values corresponding to the 
minimum and maximum vales of the data displayed on the plot. The 
second option, "Smooth Data", sets the minimum and maximum values 
and the number of gridlines to even numbers that are more 
appealing. None of the other six options are currently defined. 

Values for both the x- and y-axes can be set using the text 
fields at the bottom of the dialog. If the "Scale Data" option is 
not active, then the minimum and maximum values of either axis 
can be set manually. Note, however, that both minimum and maximum 
must be set for each if any are to be set. The number of 
qridlines and labels for either axis can be set in the same way. 
The minimum number of gridlines is 1, and the maximum is 20. The 
borders affect the way the plot is sized and drawn. The text 
fields in the x-axis area control the size of the left and right 
borders, and the text fields in the y-axis area control the 
bottom and top borders. All plot decorations appear outside the 
main plot area, so the borders may have to be resized to give the 
plot an appropriate appearance. 

The Format field dictates the number of decimal points and 
total field width that is used in displaying the axis labels. 
These values should be specified in terms of standard C 
formatting patterns as used with the standard input/output 
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• rpv,^c format is very similar to the syntax 

£U ^ C fi? n REAL variables in ^TR^J FoAt specifications. The 
commands availlKe in the Format Dialog are summarized in Table 


VII . 4 . 

Tab!** VII. 4: Format Dialog Commands 

Command 

Action 

Prev 

Switches to the previous or next format specification. Any of 
the formats numbered [0-7] can be_edited at any time. _ 

Next 

Load 


OoVvJ 

Cancel 


Accept 

Changes the properties of the plot corresponding to the 
rnrrent format to the values input. — " " 

Revert 

112 to either the previous or next format and 

then switching back. — . — — ““ 


T>Vi«a riata Attributes dialog allows the user to customize the 
appearance of the data displayed in the main content area o 
plot. This dialog is shown m Figure VII. 6. 
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T £i S ."SiSS c Sle 9 lS? 0 ? n ol n S l5Si t SnSiS B 'SS i “ ad 1 in 

variable name will be Ranged « this color whenever the dialog 
is updated. The number of dashes can ^ |P^“ S the nlmber of 

column. The format for entrie d bv the number of pixels 

pixels that should be drawn on by the ^ foll £ wed 

that should be drawn "off . For instance, ^ P 

by 10 pixels off «?“ ld *%?f ^^tL'dlaJog, since it 

onl^ref lect^additional infusion about the format currently 
selected with the Format dialog. 

The File Edit dialog provides an editable text 
which the user can load, view, make changes to, and save ASCII 
text files. This dialog is shown in Figure VII. . 
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The dialog is modeled after the xedit application included 
in the standard distribution of the X Windows from MIT. This 
rather simple dialog contains three action buttons, a single-line 
text field for entering filenames, and a large text field for 
displaying the contents of text files. This dialog is included to 
provide the user a simple capability of editing text files from 
within the application, without having to guit and restart the 
interface program or edit field in another shell tool. The load 
and save buttons read in or write out their contents based on the 
filename specified in the single-line text field. The contents 
either of the text buffer or the external text file will be 
overwritten upon executing one of these commands. The commands 
available in the File Edit dialog are summarized in Table VII. 5. 
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Table VII. 5: File Editor Dialog Commands 


Command 

Action 

Done 

Dismisses the dialog. Maintains the current state 

of the text displayed, but does not save those 
chanqes to a file. 

Save 

Saves the current state of the text to the 
filename given in the Filename text field. — _ 

Load 

Loads the contents of the filename specified in 
the Filename field into the current text buffer. 
Overwrites the contents of the text buffer with 
the contents of the file specified in the Filename 
field. If the file does not exist (i.e., if the 
filename represents a new file), then that file is 
opened (if possible) and the contents of the text 
buffer are cleared. — 


The View Dataset dialog provides a simple read-only text 
field area displaying the currently loaded dataset. This dialog 
is shown in Figure VII . 8 . 


— 

Dataset: . ./source Afcavenv. plot 


8 25 4 

Pressure 
Tenperature 
Mach Nunber 
Burn Rate 
Distance 
Mass Flou 
Surface Tenp 
Reynolds 
Exit Pressure 
Mass Burned 
Mach Nozzle 
Thrust 

O.I 6 G 66 7 332.26 0.0149 0.0000 0.0000 0 j 

57*97 12.53 301.09 0 . 01 G 8 0.0000 0.0000 0 , 

113*28 12 45 299.09 0.0159 0.0000 0.0000 0. _ 

168*58 12.45 299.09 0.0158 0.0000 0.0000 0| 

t - — r 71 _ 


Figure VII. 8: xplot View Dataset Dialog 


This dialog is provided only to allow the user to easily 
view the values as data in the plots created. This dialog 
functions just like the File Edit dialog, except it always 
contains the text of the current dataset, and the text m it can 
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h>© neither edited nor saved. The dialog can be resized according 
to the width of data it displays. There are no commands 
associated with the View Dataset dialog. It must be shown or 
hidden from either the Dataset menu or the Window Manager dialog. 

A simple command-line interface (or CLX) is provided for 
those users more comfortable with the more conventional means of 
entering commands . Note that each of the commands available 
through the CLI corresponds to one of the commands available in 
the graphical interface. A summary of the commands available is 
given in Table VII. 7. In Table VII. 7, the notation $1 refers to 
the first argument, $2 to the second, etc. 
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Command 



Table VII.7: Command-Line Interface Verbs 

Description 


Changes the a ctive plot to the one specified in $1. 

Changes the color associated with the data specified in $1 
to the color specified in $2. Equivalent to changing the 
appropriate "color" field on the Line Format dialog. 

Loads the file specified in $1 as the current datase t_. 

Decrements the plot specified in $1. If no argument is 
given, the active plot is decremented. Equivalent to 
selecting "Decrement" in the Plot Manager dialog. 

Loads the file specified in $1 as the text to be edited in 
the File Edit dialog and raises this dialog. If no file is 
given, then the previous filen ame is used. 

Executes the Caveny program from the path specified in $1 
and the application name specified in $2. If neither 
argument is given, the defaults are taken from the values 
most recently specified in the E xecute dialog. 

Changes the format being edited in the Format dialog 
to the value specified in $1. • 

Hides the plot specified in $1. Equivalent to selecting 
"Hide" for the appropriate plot in the Plot Manager dialog. 

Increments the plot specified in $1. If no argument is 
given, the active plot is incremented. Equivalent to 
selecting "Increment" in the Plot Mana g er dialog. 

Loads the file specified in $1 as the input dataset to be 
used in the Input File Editor dialog and raises this Dialog. 
If no file is specified, then the dialog is raised with a 
new input dataset. 

oath Loads a dataset from the full pathname given in $1, If no 

P argument is given, it attempts to re-load the dataset most 

recently used with a "load" command from th e Dataset — menu. 

Sets a time/space plane on the currently active plot to $1. 

If the value specified in $1 is less than one or more than 
the time/space planes available, no action is taken- 


Creates the plot specified in $1 


Terminates execution of xplot 


Saves the current format specified in $1. 

Shows the plot specified in $1. Equivalent to selecting 
"Show" for the appropriate plot in the Plot Manager dialog_. — 

Toggles the plot specified in $1 to being a space plot. * f 
no argument is given, the active plot is toggled. Equivalent 
to toggling the "Space plot/Time Plot button on the Plot 
Manager Dialog. — 

Toggles the plot specified in $1 to being a time plot. If no 
argument is given, the active plot is toggled. Equivalent to 
toggling the "Space plot/Time Plot" button on the Plot 
Manager Dialo 
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The xplot application has a number of external files that 
must be present to run properly. In addition, a number of 
jj- t- i nnal files can be used to configure the application 
according to the user's preferences, or to provide additional 
support for interactions with other programs. 

In each of the following example listings, a vertical bar 
and all tJxt following it should be interpreted as comments and 
should not be inserted into the input files. 

XPlot clr This file specifies the colors available to be 
used in the line format dialog. The file is composed of a number 
of lines containing a number and a valid color name as specified 
in the X11R4 color database . Blank lines in the input file are 
ianored* These°number-color pairs allow the user to define the 
colors that will be available for the various traces of data, 
example file excerpt is given as Table VII. 8: 


Table VII. 8: 


1 

2 


White 

Black 


Xplot. clr. File Format 

Color #1 = White (255, 255, 255) 
Color #2 = Black ( 0, 0, 0) 


18 Red 

19 Green 

20 Blue 


Color 

Color 

Color 


#18= Red 
#19= Red 
#20= Red 


(255, 
( 0, 
( 0 , 


0 , 

255, 

0 , 


0 ) 

0) 

255) 


The XPlot.var file contains entries for all of the data to 
be placed in a Caveny program input dataset generated by t 
xplot program. Each variable entry in this file consists of a 
number of parameters surrounded by a leading and a trailing 
double-quote. The separate fields are terminated by a trailing 
slash (/) The fields of a variable entry are (in order) . the 
variable list index, the FORTRAN name of the variable, A short 
definition of the variable, its default value (in appropriate 
units), and an optional long description. An example format of 
the entries is given in Table VII. 9. 


Table VII. 9: 


0 / 

NDELX/ 

Number of spacewise steps (dx) 
along the propellant grain/ 
20 / 

Maximum value = 29 

Format: Integer Units: None 


XPlot.var File Format 
| Leading double-quote 

I Variable list index 

| FORTRAN variable name 

| Short variable description 

| Default value 

| Long variable description 


Trailing double-quote 
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override these settings . 

Plot format files. Unlike the input files already described, 

ran be any number of plot format files, and they can have 

any naSe The default format file, which is read upon Program 
any name, ine f Th: i s fii e is a simple text file 

start-up, 1S f .?^ Xn ^ nt i ons selected in a plot format dialog. Each of 

enumerating the options^elected.i ^ che file , in the 

the fields store . ^ dialog. Although there is no 

reason'why Se'SErlSSd nSt edi? thil file using either the 
input File dialog^r^ IfP-^J^^/^creaied by'the Plot 

dialog directly. A sample listing for the XPlot.def file 
is given as Table VII. 10. 


Table VII. 10: 


Plot #0 
0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

64 54 32 32 

10.4 10.4 
10 10 


0 

0 

0 


Sample XPlot.def File 
Plot title 
Space data shown 
Time data shown 
Options selected 
Border (L, T, R, B) 

Axis label formats 

Number of axis gridlines/labels 


The example plots shown in this section were all produced by 
xplot T with data generated by the Caveny pr. ogr, am • Al 1 ® 

nlots used the default values for format attributes from G ® 

XPlot clr and Xplot.def files' The setup required for the plots 
XPlot.cir ana Ay following the plot. The data for these 

pL« S Sas generated by the Caveny code from the sample input frle 
given as Appendix A. 

Fiaures VII 9-VII.12 show the development of pressure in ^ he 

ianiter motor pressure wave down the length of the grain, 
igniter mo p shown to the same scale to 

IE Is :s esus s s, 

down the chamber . 
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Figure VII. 9: Pressure versus Location at TIme=30 ms 
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Figure VII. 11: Pressure versus Location at Time=90 ms 
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Figure VII. 12: Pressure versus Location at Time-120 ms 
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These four plots are typical of the default type produced by 
xplot. All of the data are presented as a function of distance 
along the motor geometry. Figures VII. 9 - VII. 12 represent four 
time slices of the entire motor. Each plot contains a maximum of 
30 points per curve, since this number is the maximum that the 
Caveny program allows for axial stations. Pressure was the only 
variable selected from the format dialog. The time was set using 
the Increment and Decrement buttons on the Plot Manager dialog. 
The scaling and number of gridlines were set manually. The number 
of decimal places for the x-axis were left at their default value 
of zero, but the y-axis labels were set to a 0.3 format. 

Figures VII. 13 and VII. 14 show the data generated by the 
Caveny program in a different fashion. These plots illustrate the 
time history of a scalar quantity at two different axial 
distances along the motor. In this particular example, the 
chamber pressure is displayed as a function of time at the head 
end and aft end stations. Each of these plots contains 126 data 
points, since the Caveny input file was configured using its 
print interval and total simulation time to generate this number. 
This viewing option is set using the Space Plot /Time Plot toggle 
button on the Plot Manager dialog. With this option, the user can 
quickly switch back and forth in the manner chosen to view the 
data. The station is selected using the Increment and Decrement 
buttons on the same dialog. The scales, grids, and title were 
each set manually. 
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Figures VII. 15 and VII. 16 show that multiple traces of data 
can be displayed on a single plot. In these two figures, the Mach 
number and grain burn rate are shown as functions of distance at 
two different times in the simulation (90 and 180 milliseconds) . 
Of interest in the first plot is the spike in Mach number and the 
sharp drop-off in burn rate near the aft end of the motor. The 
spike again corresponds to the igniter discharge, and the burn 
rate reflects the fact that the aft end of the motor has not yet 
begun to burn. By 180 milliseconds, both curves have flattened 
out to more stable values all along the motor. In each plot, the 
legend reflects the colors associated with the traces for the two 

variables . , , , , 

To configure these two plots, Mach number and burn rate were 

the only variables selected using the Format dialog. To make the 
traces print more distinctly, the colors and dashes for burn rate 
were changed from their default values. Using the Format Lines 
dialog the burn rate color was set to black, and the dashes were 
set to' 20 on 10 off (by entering "20 10" in the variable's 
appropriate field) . In these two figures, the plots were manually 
titled and automatically scaled. 
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Figures VII. 17 and VII. 18 are presented as examples of 
variables that are functions only of time and 

the grain. In these two figures, the Reynolds number at the exit 
of the motor and the thrust generated by the motor ^Plotted. 
The curve for the thrust produced is in accordance with the plot 
of aft end pressure shown in Figure VII. 14. The thrust is zero 
until the igniter gases begin to be discharged at approximately 
ef milliseconds. The Reynolds number at exit is a predominantly 
linear function with the exception of the ignition spike from 
to 100 milliseconds. 

To create plots of variables that are time-dependent only, 
the Time Plot option must be selected from the Plot Manager 
dialog If thiO option is not chosen, no data will be dravm in 
the Dlot frame Note that the Increment and Decrement buttons 
this dialogic no? function if only time -dependent variables are 
selected In all other aspects, these variables are configured 
and displayed in exactly the same way as the other variables when 

plotted as functions of time. 


VII-22 



Figure VII. 17: Reynolds Number at Exit versus Time 
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APPENDIX A 

Sample Caveny Program Input File 


A-l 



FLAME SPREAD CONTROL AT HEAD END , DYNAMIC THRUST CHECKOUT 

= 0 . 002 , 

= 0 . 001 , 

= 0.00016, 

= 24, 

= 298.89, 

= 0.84717, 

= - 2 . 0 , 

= 2685.4, 

= 3353.0, • 

=28.18, 

= 1.136, 

= 2325.81, 

= 2.67, 

= 1329.972, 

= 0 . 01 , 

= 0 . 0011 , 

= 1.758, 

= 0.3, 

= 0.0009939, 

= 1.054, 

= 68.08, 

= 0 . 0 , 

= 0.33, 

= 145.64, 

= 0.97844, 

= 12.31, 

= 1.1348, 

= 0.0027, 

= 0.00906, 

= 0 . 01 , 

= 19, 

= 25, 

= 3, 

= 0 . 002 , 

= 0.1, 0.1, 0.1, 27*1.0, 

= 1.4, 1.3, 1.2, 1.1, 26*1.0, 

= 0.010, 0.007, 0.005, 27*1.0, 

= 0.94, 

= 850.0, 

= 0 . 1 , 

= 10 , 

= 1 , 

= 20 


46 

fie END 


47 

48 

2 

0.0000 

0.0000 

49 


1.0000 

1.0000 

50 

51 

2 

0.0000 

0.0000 

52 


1.0000 

1.0000 

53 

54 

2 

0.0000 

0.0000 

55 


1.0000 

1.0000 

56 

57 

2 

0.0000 

0.0000 

58 


1.0000 

1.0000 

59 

60 

2 

0.0000 

0.0000 

61 


1.0000 

1.0000 

62 

63 

2 

0.0000 

0.0000 

64 


1.0000 

1.0000 

65 

66 

2 

0.0000 

0.0000 

67 


1.0000 

1.0000 

68 

2 




1 

MM- 1 17DEC77 

2 

fie NAME 

3 

TMAX 

4 

TPRINT 

5 

DELTAT 

6 

NDELX 

7 

TPI 

8 

PAM 

9 

UNIT 

10 

TIGN 

11 

TFREF 

12 

W 

13 

GAMA 

14 

AT 

15 

XP 

16 

XE 

17 

RUFSUR 

18 

FKPR 

19 

ROPR 

20 

CPR 

21 

SIGP 

22 

RREF 

23 

PREF 

24 

EBC 

25 

BREXP 

26 

DE 

27 

CM 

28 

ALFAD 

29 

GAMAN 

30 

PISUBK 

31 

DELTTF 

32 

DFSDT 

33 

NIGTAB 

34 

NAPDVX 

35 

NDATA 

36 

TPRINT 

37 

BF1 

38 

CHC 

39 

DELTFF 

40 

DDHC 

41 

TPSCRI 

42 

FCRIT 

43 

NBL 

44 

NPUNCH 

45 

NBPTAB 
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69 


0.0000 

0.0000 

70 


1.0000 

1.0000 

71 

2 



72 


0.0000 

0.0000 

73 


1.0000 

1.0000 

74 

2 



75 


0.0000 

0.0000 

76 


1.0000 

1.0000 

77 

2 



78 


0.0000 

0.0000 

79 


1.0000 

1.0000 

80 

2 



81 


0.0000 

0.0000 

82 


1.0000 

1.0000 

83 

2 



84 


0.0000 

0.0000 

85 


1.0000 

1.0000 

86 

2 



87 


0.0000 

0.0000 

88 


1.0000 

1.0000 

89 

2 



90 


0.0000 

0.0000 

91 


1.0000 

1.0000 

92 

2 



93 


0.0000 

0.0000 

94 


1.0000 

1.0000 

95 

2 



96 


0.0000 

0.0000 

97 


1.0000 

1.0000 

98 

2 



99 


0.0000 

0.0000 

100 


1.0000 

1.0000 

101 

2 



102 


0.0000 

0.0000 

103 


1.0000 

1.0000 

104 

2 



105 


0.0000 

0.0000 

106 


1.0000 

1.0000 

107 


0 . 

1 . 

108 


0.02 

18.11 

109 


0.028 

46.95 

110 


0.044 

406.05 

11 JL 


0.052 

484.35 

112 


0.064 

514.35 

113 


0.07 

516.69 

114 


0.088 

510.06 

115 


0.18 

427.47 

116 


0.24 

392.42 

117 


0.288 

371.00 

118 


0.344 

348.43 

119 


0.356 

339.86 

120 


0.38 

313.75 

121 


0.424 

275.18 

122 


0,46 

253.38 

123 


0.508 

232.73 

124 


0.556 

219.5 

125 


0.566 

0.0 

126 


2.670 

2538.440 

127 


57.970 

2543.620 

128 


113.280 

2561.700 

129 


168.580 

2706.870 

130 


223.890 

2612.100 

131 


279.190 

2645.000 

132 


334.500 

2658.850 

133 


389.800 

2705.320 

134 


445.100 

2774.530 

135 


500.410 

2844.420 

136 


555.710 

2915.190 

137 


611.020 

3029.790 


1154.450 

1154.45 

31 . 

1156.490 

1156.49 


1144.950 

1144.95 


278.930 

278.93 


181.180 

181.18 


182.310 

182.31 

329.52 

182.790 

182.79 


184.380 

184.38 


186.720 

186.72 


189.060 

189.06 


191.400 

191.40 


195.130 

195.13 

649.52 


22439 . 


877 . 


2896 . 


59741 . 


20151 . 


20003 . 


2.21 


1.5 


1.5 
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.38 

666.320 

2757.930 

.39 

721.630 

2720.090 

.40 

776.930 

2789.300 

.41 

832.230 

2859.380 

.42 

887.540 

2930.330 

.43 

942.840 

3072.280 

.44 

998.150 

3025.310 

.45 

1053.450 

3224.030 

.46 

1108.760 

3731.700 

.47 

1164.060 

4316.430 

.48 

1219.360 

4946.200 

.49 

1274.670 

5404.350 

.50 

1329.972 

5594.670 

.51 

0.6808 

0.1002 

.52 

4.7651 

0.4383 

.53 

68.08 

1.054 

154 

&NAME 


L55 

TMAX = 

o 

o 

L56 

DELFAC 

= 1.3 

L57 

Sc END 


L58 

&NAME 


L59 

TMAX = 

0.30 

L60 

& END 


L61 

ScNAME 


L 6 2 

TMAX = 

0.40, 

L 6 3 

TPRINT 

= 0.005 

L64 

& END 



186.160 186.16 

184.880 184.88 

187.220 187.22 

189.560 189.56 

191.890 191.89 942. 

196.490 196.49 

194.980 194.98 

201.280 201.28 

216.550 216.55 

232.900 232.9 

249.310 249.31 1293.2 

260.600 260.6 

265.150 265.15 . 
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APPENDIX B 

Caveny Program Makefile 



1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 


# File: makefile 

# This file contains all dependencies necessary 



F77 = f 77 
FFLAGS= -g 

OBJECTS = \ 

ans . o\ 
breal . o\ 
fecal . o\ 
fet . o\ 
f ct2 , o\ 
filt.o\ 
hccal . o\ 
inhomo . o\ 
iterpl . o\ 
lagin2 . o\ 
lbc.o\ 
migdot . o\ 
pnch . o\ 
prepc . o\ 
prepp. o\ 
rbc . o\ 
redpun. o\ 
runge . o\ 
sdata.oX 
sdatac .o\ 
slotk . o\ 
solute .o\ 
solutp . o\ 
srnbug. o\ 
thrusl . o\ 
timest.oX 
trid.oX 
decomp .o\ 
solv . o\ 
gsolv. o\ 
asolv. o\ 
vsolv. o\ 
rep!2 . o\ 
repll . o\ 
sing. o\ 
tscal . o\ 
main.o 


caveny : $ ( OBJECTS ) 

$ ( F77 ) $ (FFLAGS) -o caveny $(OBJECTS) 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 


xpCommands . c 

*/ 

♦include ■xplncs.h* 

♦include "xpDefs.h" 

♦include "xpwidgets . h" 

♦include "xpgeneral . h v 
♦include "xpmenus.h' 

♦include "xpformats . h" 

♦include 'xpplots.h* 

♦define cSpace 32 
♦define cTab 9 

/* $1: Global variables */ 

char cli [64] = { “ } ; a , . „ „ . „ „ . . 

char cliString[16] (64] = {""#**#* " * * * * 
int clilndex=0 , cliNum=0 ; 

f* $2: Local variables - (none) */ 

char coxnmandString[16] [32] , *commandStack[16J ; 
char inputstring [64 ] ; 

/* $3: Global functions */ 

void ListCallbackO ; 
void SaveCommand { ) ; 
void RecallCommand ( ) ; 
void Echo ( ) ; 

/* $4: Local functions */ 

void Push ( ) ; 

void Breakdownstring 0 ; 

void DoParse{); 

/* $5: Implementation e======s===*”-====* S!S "== =:3s========== 

void ListCallback (w, client_data, call_data) 

Widget w; XtPointer client_data T call_data; 

XawListReturnStruct *r = (XawListRetumStruct*) call_data; 
char s [ 64 ] ; 

strcpy (s , ( *r ) .string) ? 

XtVaSetValues (cliText, XtNst ring, s, NULL) ; 
XawTextSetlnsertionPoint (cliText , str len (s) ) ; 

) 

void SaveCommand (w, even t, par ams , numParams ) 

Widget w; XEvent ‘event; String ‘params; Cardinal ‘numParams? 

( 

XawTextPosition n; 
char s [ 64 ] , *ss; 
int num; 

XtVaGetValues (cliText, XtNstring, &ss,NULL) ; 
n=strlen (ss) ; 
if (n>0 ) { 

if (n>63) n= 6 3 ; 
stmcpy (s, ss # n) ; 
s [n] = ' \0 ' ; 
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69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 


} 


Push(s) ; 
cliNum++; 
clilndex=0 ; 

XtVaSetValues (cliText, 

XtNstring, * * , XtNinsertPosi tion, strlen 
printf ( ■ >cli : %3d %s\n" , cliNum, s) ; 


) , NULL) ; 


strcpy ( inputstring , s ) ; , . 
BreakdovmStr ing ( inputstring , commands tring , fcnum) 
Do Par se (num, commandstring} ; 


void Recal 1 Command (w, event , par ams ,numParams) *^^' 0 *^*™* * 

Widget w; XEvent ‘event; String ‘params; Cardinal numParams, 

if (*numParams>0) 

if (params [0] [0]==|1 , ){ 

if ( (cliIndex<cliNum) && ( clilndex<15 ) ) { 

clilndex++ ; 


} 


XtVaSetValues (cliText, XtNstring , cliString [clilndex] . NULL) ; 
XawTextSet Insert ionPoint (cliText, strlen ( cl iStnng [ cl i Index] ) ) , 


} 


if (params [0] [0] == ' 2 ’ ) { 

if (clilndex>0) { 

XtVaSetValues (cliText / XtNstring, cliString [clilndex] ,NULU ; 

XawTextSet Insert ionPoint (cliText, strlen (cliStr mg [clilndex] ) ) , 

} 

} 


void Echo(s) 

String s; 

int sLength=strlen (s) , code; 

XawTextBlock text; 

text . f irstPos=0 ; text . length=sLength; text .ptr^s; text . f ormat-FMT8BIT 

code = Xav/TextReplace(logText, startPos, endPos,&text) ; 

startPos = startPos+sLength; 

endPos = startPos; ^ . 

XawTextSetlnsertionPoint ( logText , startPos) ; 


void Push(s) 
String s; 

{ 

int i ; 


for ( i = (NUMcliString-1 ) ;i>l;i — ) 

strcpy (cliString [ i ] , cliString [ i-1 ] ) ; 
strcpy (cliString [1] ,s) ; 
strcpy (cliString [0] , *■) ; 

for ( i=0 ; i<NUMcliString; i++) . 

commandStack [NUMcliString-i-1] =cliString[i] ; . 

XawListChange(commList, commandStack, 15, 0, False) ; 


void Breakdownstring (old, new, count) 

char old[); char newt) 132]; int ‘count; 

int i=0 , m=0 , n=0 ; 


while (old [ i ] ! - 1 \0 1 ) ( 

while ( (oldti] ==cSpace) 


(oldti) ==cTab) ) 


while” ( (oldti) !=cSpace)&&(old[i] !=cTab)tS.(old(i) 


\0 ’ ) && (n<32 ) ) ( 
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138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 
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155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 
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173 

174 

175 

176 

177 

178 

179 

180 
181 
182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 


new[m] [n] = old[i]; 
i++ ; 
n ++ ; 

} 

new [m] [n] « 1 \ 0 ' ; 
m+ + ; 
n=0 ; 

) 

*count=m; 


void DoParse ( argc , argv) 

int argc; char argv [3 [32]; 

{ 

int l,m,n; 
if (argc>0) { 

if (strcmp (argv [0] , "quit" ) ==0} { 

exit (0) ; 

else if (strcmp {argv [0] , "change" ) ==0) { 

printf { * {<change> un implemented) \n" ) ; 

else if (strcmp (argv[0] , "decrement* ) ==0) { 

printf (•{<dec[rement]> un implemented) \n" ) ; 

else if (strcmp (argv[0) , "edit" ) ==0) { 

printf { ■ (<edit> unimplemented} \n" ) ; 

else if (strcmp (argv [0] , "execute" ) ==0) { 

printf C {<exec [ute] > un implemented) \n" ) ; 

else if (strcmp (argv [0] , "hide")==0) { 

if (argc>2) { 

if (strcmp (argv [1] , "plot* ) ==0) { 

n=atoi (argv[2J ) ; , r 

if ( (n>=0) && (n<=MAXPLOTS) && (plotExists [n] ) ) 
ShowHide2 (dummy , wPlotO+n, False, 0) ; 

) 

) 

else if (strcmp (argv [0] , • increment" ) ==0) { 

printf { “ (<inc [rement] > un implemented) \n“ ) ; 

else if (strcmp (argv [0] , "kill* ) ==0) { 

if (argol) { 

if (strcmp (argv[l] , "plot ") ==0) { 

if (argc>2 ) { 

n=atoi (argv[2] ) ; 

if ( (n>=0) && (n<MAXPLOTS)&& (plotExists [n] ) ) { 

KillPlot (n) ; 

) 

) 

) 

) 

else if (strcmp (argv [ 0] , * load" ) ==0 ) { 

if (argol) { 

if (strcmp (argv[l] , "data" ) ==0) { 

LoadData (dsName) ; 

) 

) 

else if (strcmp(argv[0] , "make")===0) { 

if (argc>2) { 

n=atoi (argv [2 ] ) ; 

if ( (n>=0 ) && (n<MAXPLOTS) && ( IplotExists [n] ) ) { 

MakePlot (n) ; 

) 
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207 

208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 
1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 


} 


else if ( strcmp (argv [ 0 ] , "print* ) ==0) { 

printf ( * {<print> un implemented } \n 9 ) ; 


else if (strcmp (argv[0] , 'save* ) -=0) 
printf ( " {<save> un implemented} \n* 

else if (strcmp (argv [0] , “set* ) ==0) 
if (argol) { 

if (strcmp (argv [1] , "data* ) --0) 
if (argc>2) { 

strcpy ( dsName , argv [ 2 ] ) ; 
SetStatusLabel ( ) ; 


{ 

) ; 
{ 

{ 


} 

else Selectstring ( 1 ) ? 


if (strcmp (argv [1] , " format* ) -=0) { 

if (argc>2) ( 

strcpy (pfName , argv [ 2 ] ) ? 
SetStatusLabel { ) ; 


else Selectstring ( 2 ) ; 

) 

) 

else if (strcmp (argv [ 0 ] , * show" ) ==0) { 

if (argc>2 ) { , 

if (strcmp (argv [ 1 ] , “plot ■ ) ==0) ( 

n=atoi (argv [2] ) ; _ v , 

if ( (n>=0 ) && (n<=7 ) && (plotExists [n] ) ) 
ShowHide2 (dummy , wPlotO+n, True , 0) ; 


} 

} 

else if (strcmp (argv [0] , "view*)==0] I { 
printf ( • {<view> unimplemented} \n“ ) ; 

else printf ('{unrecognized command [%s] [%d] . . . )\n* , argv[0] , argc-1 ) 


/* 

xpData . c 

*/ 


♦include *xplncs.h' 

♦include "xpDefs.h" 

♦include "xpplots.h" 

♦include "xpwidgets .h - 

/* $1: Global variables */ 

Boolean dataLoaded=False; 

char columnNames [2*HAXCOLS] [24] ; 

int numPlanes / numCols , numRows # numTimeDep # 

int currentPlane [MAXPLOTS ]={0, 0,0, 0,0, 0,0,0}; 
int currentRow [MAXPLOTS ]={0, 0,0, 0,0, 0,0,0); 

float x [MAXROWS ] ; 
float z [MAXPLANES] ; 


float plotData [MAXPLOTS] [2*MAXCOLS] [MAXPLANES] ; 
float dMin [MAXPLOTS] [2 *MAXCOLS] ; 
float dMax [MAXPLOTS ] [2*MAXC0LS] ; 
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28 

29 
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50 

51 
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54 

55 

56 

57 

58 

59 
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61 
62 
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64 

65 

66 

67 
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70 

71 

72 

73 
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81 
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float xMin , xMax , xRange ; 
float zMin , zMax , zRange ; 

float yMin [MAXPLOTS ] ,yMax[MAXPLOTS] ,yRange[MAXPLOTS] ; 

/* $ 2 : Local variables */ 

float d [MAXPLANES J [MAXCOLS+1] [MAXROWS] ; 

/* $3: Global functions */ 

void LoadData ( ) ; 
void PrintDataO; 
void CopyDataO; 

/* $4: Local functions */ 


void StatlndO; 
void StatDataO; 
void PrintPlotData { ) ; 

/* $5: Implementation ====== 


void LoadData (infile) 
char *infile; 

( 

FILE *f i ; 
int i, j , k,n=0? 
float dummy; 
char s [ 64) ; 


printf(“>Now loading dataset from file: <%s>\nMnfile) ; 
fi = fopen(infile, *r* ) ; 
i f ( f i ! =NULL) { 

strcpy ( s , * Dataset : • ) ; 

strcat (s , inf ile) ; _ . 

XtVaSetValues (dataShell , XtNti tie , s , NULL) ; 

XtVaSetValues (dataText , XtNstring , inf ile, NULL) ; 

for ( i=0 ; i<8 ; i++ ) . , vtttt t \ 

XtVaSetValues (PMitemU) [i] , XtNsens it ive, True, NULL) , 

dataLoaded=True; . _ . 

fscanf (fi, - %d%d%d%d\n" , tnumCols , fcnumRows f &numTimeDep) ; 

for (k=0 ; k<numCols ; k+ + ) { 

fgets (columnNames [k] , 25, fi) ; 

for (i=0; i<strlen (columnNames [k] ) ; i++) , , tV _, 

if (columnNames [k] [i]= =, \n') columnNames [k] [i]= \0 
XtVaSetValues (FEItem(l) [k] , XtNlabel , columnNames [k] , 

XtNsens itive, True, NULL) ; 

XtVaSetValues (CFItem[0] [k] , XtNlabel , columnNames [k] , 

XtNsensitive , True, NULL) ; T . 

XtVaSetValues (DAItem[0) [k] , XtNlabel , columnNames [k] , NULL) , 

} , i r 

for (k=0;k<numTimeDep; k+ + ) ( 

fgets ( columnNames [ MAXCOLS+k ] , 2 5 , f i ) , 

for ( i=0 ; i<strlen (columnNames [MAXCOLS+k] ) ; i++) 

if (columnNames [MAXCOLS+k] [i] == 1 \n ' ) 

columnNames [MAXCOLS+k] [ i ]=’ \0 ' ; , , 

XtVaSetValues (FEItem[2) [k] , XtNlabel , columnNames [MAXCOLS+k] , 

XtVaSetValues*CFItemTo^HAXCOLS+k] .XtNlabel , columnNames [MAXCOLS+k] 

XtVaSetviluestDAIt^^ts+k] .XtNlabel .columnNames [MAXCOLS+k] .NULL) 


if { (numRows>0) && (numCols>0)) { 
k=0 ; 
do { 

n=fscanf (fi, *%e B ,&z[k] ) ; 
if (n ! =EOF) { 

for ( i=0 ; i<numRows ; i++) ( 

if (k==0) 
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131 
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134 
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136 

137 
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141 
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143 
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fscanf { f i , " %e" , &x [ i ] ) ; 
else 

fscanf (fi, "%e" , &dummy) ; 
for { j=0; j<numCols; j++) 

fscanf (fi, "%e" , &d [k] ( j 3 [i] ) ; 

> • ^ 
for ( i=0 ; icnumTimeDep ? i++ ) 

fscanf ( f i , -%e',£cd[k] [MAXCOLS] [i]) ; 

} 

} while { (n ! =EOF } && ( k< (MAXPLANES-1 ) ) ); 

} 

numPlanes =k-l ; 
fclose (f i) ; 

StatlndO ; 

else printf{ "Unable to load dataset (does not exist) \n"); 


void PrintData (outf ile) 
char * outf ile ; 

{ 

FILE * for- 
int i,j,k; 


printf ( " >Now printing dataset to file: <%s>\n" , outf ile) 

fo = f open (outf ile , "w" ) ; * 

fDrintf (fo, *%4d %4d %4d\n", numPlanes ,numCols , numRows) , 
if ( (numRows>0) S.& (numCols>0) <numPlanes>0) ) 

for (k=0;k<numPlanes;k++) { 

fprintf (fo, -%12.6f\n" ,z[k] ) ? 
for ( i=0 ; i<numRows ; i++) { 

fprintf (fo,*%12.6f*»x[iJ); 
for ( j=0 ; j<numCols ; j++) r 

fprintf (fo, •%12.6f",d[k] [j] [i]) ; 
fprintf ( fo, " \n* ) ; 


} 


} 


fclose ( fo) 


void CopyData (plotNum) 
int plotNum; 


{ 


int i , j ; 

if (planar [plotNum] ) 

for ( j=0 ; j<numCols ; j++) 

for (i = 0; icnumRows; i++) „ , , 

plotData [plotNum] [j] [ i ] =d [currentPlane [plotNum] ] [ 3 ] [ 1 ] 


else { 

for ( j=0 ; jcnumCols ; j++) ^ 
for ( i=0 ; icnumPlanes ; i++ ) 

plotData [plotNum] [ j ] [ i ] =d [ i ] 
for ( j =0 ; j <numTimeDep ; j ++ ) 
for ( i=0 ; i<numPlanes ; i++) 

plotData [plotNum] [MAXCOLS+j] 


[j] [currentRow [plotNum] ] ; 
[ i] =d [ i] [MAXCOLS] [j] ; 


) 

} 


void StatData (p) 
int p; 

{ 

int i, j ,m,n, test; 
float Max, Min; 
int count; 


count 

test 


(planar [p 3 ? numRows : numPlanes) ; 
(planar [p] ? MAXCOLS : 2 * MAXCOLS); 
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for ( j = 0 ; j <test ? j++) 
if (useCol [p] [ j ] ) ( 

Min= 1000000; 

Max- -1000000 ; 

for ( i=0; i<count; i+ + ) ( 

if (plotData Ip] [j) [ i ] <Min) 
if (plotData [p] [j] t i] >Max) 

) 

dMin [p] [ j ] =Min; 
dMax [ p ] [ j ] =Max ; 


Min=plotData [p] t j] [i] ; 
Max=plotData [p] Cjl C i 3 ; 


Min= 1000000; 

Max=- 1000000 ; 

for ( j=0; j<test? j++) 

if (useCol [p] [jl) ( , f 

if (dMinlp] t j]<Min) Min=dMm[p] 13 J ; 
if {dMax [p] [ j] >Max) Max=dMax[p] [j] ; 

} 

yMintp] =Min; 
yMax [ p ] =Max ; . 

yRange [pj =yMax(p] -yMm[p3 ; 
if (yRange [p] <=0 . 001) { 

yRange [p] =1 . 0 ; 
yMax [ p ] =yMin[p] +1. 0; 

} 


void StatlndO 
{ 

int i ; 

xMin= 1000000; 

xMax = -1000000 ; 

for ( i=0 ; i<numRows ; i++) { 

if (x ( i 3 <xMin) xMin=x{i]; 
if {x [ i] >xMax) xMax=x[i] ; 

} 

xRang e = xMax - xM i n ; 

zMin= 1000000; 

zMax= - 1 0 0 0 0 0 0 ; 

for ( i=0 ; i<numPlanes ; i++ } { 

if ( z [ i ] <zMin) zMin=z[i] ; 
if (z [ i ] >zMax) zMax=z [ i ] ; 

} 

zRange= zMax- zMin ; 


void PrintPlotData (p) 
int p; 

{ 

int i,j,k; 


, curr„»«x» [pi . 

if (planar [p] ) { , 

if ( {numRows>0) && (numCols>0) ) 

printf ( “%10 . 4f \n* , z [currentPlane [p] 3 ) ; 
f or ( i=0 ; i<numRows ; i++ ) { 

printf ("%10.4f" f x[i] ) ; 
for { j=0 ; j<numCols ; j++) 

printf ("%10.4f " f plotData [p] [ 3 ] [ 1 ] ) ; 
printf C\n-) ; 


pr int f { " xMin : %10 . 4f \n" , xMin) ; 

printf ("xMax: %10 . 4f \n" f xMax) ; 

printf ( -xRange : %10 . 4f \n" , xRange) ; 


) 

else { 
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) 


if ( (numPlanes>0) && (numCols>0 ) ) 

printf p%10.4f\n" , x [currentRow [p] ] ) ; 
f or ( i=0 ; icnumPlanes ; i++ ) ( 

printf <"%10.4f,2[i]) ; 
for ( j=0; j<numCols ; j++) 

printf ( "%10 . 4f " , plo tData [plColCij)' 
printf ( *\n') ; 

printf CzMin: %10 . 4f\n" , zMin) ; 

printf ("zMax: %10 . 4f \n* , zMax) ; 

printf ( ’ zRange : %10 . 4f \n" , zRange) ; 

printf ( "yMin : %10 . 4f \n" , yMin[p] ) . 

printf CyMax: % 10 . 4f \n’ , yMax [p] ) ; 

printf ( • yRange : %10 . 4f \n' ,yRange [p] ) ; 


/* 

xp Formats . c 

*/ 

# include "xplncs.h" 

# include *xpDefs.h" 

# include "xpgeneral . h* 

# include *xpwidgets .h* 

# include ^menus^" 
finclude ■xpdata.h" 

♦ include ■xpplots.h’' 

♦define feLoad 1 
♦define feSave 2 
♦define feCancel 3 
♦define feAccept 4 
♦define feRevert 5 

/* $1: Global variables */ 

PlotFormat f[MAXPLOTS]; 
int activeFormat=0 ; 

/* $2: Local variables - (none) */ 

l* $3: Global functions */ 

void InitFormat ( ) ; 
void SetFEValues ( ) ; 
void GetFEValues { ) ; 
void RecalcFormat ( ) ; 
void FEToggleCallback { ) ; 
void FECommandCallback ( ) ; 
void SetFEValues () ; 
void GetFEValues ( ) ; 

/* $4: Local functions - (none) */ 


void PrintFormat ( ) ; 

*/ 

/ * $5: implementation ======== 

void InitFormat (p) 
int p; 

{ 

int i,j; 
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sprintf <f[p] .title, "Plot #%d",p) ; 
for { i = 0 ; i<8 ; i + + ) { 
f [p] . color [i] =i; 

f [p] . color [MAXCOLS+i ] =MAXCOLS+i ; 
f [pj .numDashes [ i ] = 0 ; 
for < j=0; j<8 ; j++) 

f [p] .dashes [ i 3 (j]=0; 
f [p] . opt [ i] =False; 

f [p] ,borderL=64; f [p] .borderR=32; f [p] .borderT=32 ; f [p] ,borderB=64 
f [p] .xNum=10; f[p].yNum=10; 

f [p] ,xMin=f [p] .xMax=0; f [p] ,yMin=f [pi .yMax=0; 
strcpy ( f [p] . xFormat , ■ .0*) ; strcpy(flp). yFormat, ’ .0') ; 
f [p] .opt [aScaleAxes] =True; f [p] . opt [aSmoothNum] =False; 
f [p] . xTick=4 ; f [p] .yTick=4; 
f [p] .width=640; f [p] . height=480 ; 

f [p] . xlnc= { f [p] .width -f [p] . borderL-f Ip] .borderR) / f [p] .xNum; 
f [p] .yinc=(f [p] .height- ftp] .borderT-f [p] .borderB) /f [p] .yNum; 
f [p] . nGrid=f Ep3 . xNum+f [p] .yNum+2; 
f [p] . xLength=f [p] .xNum*f [p] .xlnc; 
f[p] . yLength=f [ p ] ,yNum*f [p] .ylnc; 

) 

void SetFEValues (p) 
int p; 

{ 

char s[16][32],t[16]; 
int i,j; 

strcpy (s [0] , f [p] . title) ; . 

XtVaSetValues (FEItem[FEBTitlel [0] ,XtNstring, s [0) , NULL) ; 

for (i=0;i<(2*MAXCOLS) ?i++) { 

sprintf (s [i] , "%d" , f [p] .color [i] ) ; 

XtVaSetValues (DAI tem[l] { i ] , XtNstring, s [ i ] , NULL) ; 

XtVaSetValues (DAI tem[0] [i] , 

XtNbackground, lineColor [ f [p] .color [i] 3 , 

XtNf oreground, rColor [Black] .pixel, NULL) ; 

} 

for ( i=0 ? i< ( 2 *MAXCOLS) ; i++) ( 
strcpy (s [ i ] , " * ) ; 

for (j=0; j<f[p] .numDashes [i] ; j++) { 

sprintf (t, -%d ",(int) f [p] .dashes [i] [j ]) ; 
street (s [ i ] , t ) ; 

XtVaSetValues (DAI tern [2] [i] , XtNstring , s [ i ] ,NULL) ; 

) 

sprintf (t, " %%%s f " , f [p] .xFormat) ; 
sprintf (s[0],t,f[p] .xMin ) j 
sprintf {s[l],t,f[pl .xMax) ; 

sprintf (s [2] , ^d - , f [p] .xNum) ; 
sprintf (s [3 j # *%d" , f [p] .xTick) ; 

strcpy (s [4] , ftp] .xFormat) ; 
sprintf ( s [ 5 ] , ■ %d" , f [p] . borderL) ; 
sprintf (s[6] , ’%d“,f[p] .borderR ) ; 
for (i=0; i<7 ; i^+) 

XtVaSetValues (FEItem [FEBX] [i] , XtNstring, s [ i] , NULL) ; 

sprintf (t, -%%%sf,f [p] .yFormat) ; 
sprintf (s[0] ,t,f[p] .yMin ) ; 
sprintf (s [1] , t, f tp] .yMax) ; 
sprintf (s [2] , *%d" , f [p] .yNum) ; 
sprintf (s [3] , B %d“,f [p] .yTick) ; 

strcpy (s [4] ,f [p] .yFormat) ; 
sprintf (s[5] , "%d“,f[pl .borderT) ; 
sprintf (s [6] , "%d*,f [pi .borderB) ; 
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for ( 1=0 ; i <7 ; i++ ) , . t « : i t, T t t \ 

XtVaSetValues (FEItem[FEBY] [ i ] , XtNstring, s [ i ] , NULL) 


for ( i= 0 ; i<MAXCOLS; i++ ) 

1 ^XtVaSetValues (FEItem[FEBVars] [i] .XtNsensitive, True, NULL) ; 

^XtVaSetValues (FEItem[FEBVars] [i] , XtNsensitive, False, NULL) ; 
for (i= 0 ; i<MAXCOLS; i++) 

if ( i<numTimeDep) . . T x . 

XtVaSetValues (FEItem[ 2 ] [i] .XtNsensitive, True, NULL) , 

el XtVaSetValues(FEItem[ 2 ] [ i) .XtNsensitive, False, NULL) ; 
f °XtVaSetValues (FEItem[FEBVars) [i] , XtNstate, useCol [p] [i] .NULL) ; 
f 0 xiias 4 tVa^esTF?itei[ 2 ) [i J , XtNstate, useCol [p] [MAXCOLS.i] .NULL) 
f °XtVaSe tValies ( FEI tem [FEBOpts ] [ i ] , XtNstate , f (p ] . opt [ i ] , NULL ) ; 


void GetFEValues (p) 
int p; 

{ 

String t; 

char s 1 16] [ 16] ; 

int i, j, count, u [8] ; 

XtVaGetValues (FEItem[0] [0] .XtNstring, it, NULL) ; 
strcpy(f [pl .title,t) ; 


for (i = 0;i<(2*MAXCOLS) ?!+♦) ( . xtttt t \ • 

XtVaGetValues (DAItemfl] [x] , XtNstring, &t, NULL) , 


strcpy (s [ i ] , t ) ; 
j * atoi (s [i] ) ; 

f [p] . color [i] * ( { j>=0) fitfit (3 <= 31) 


? j : 


i) ; 


} 


for <i = 0;i<(2*MAXCOLS) ;i + +) { 

XtVaGetValues (DAI tem[2] [i] , XtNstring, &t , NULL) ; 

count*sscanftsIi] , "%d %d %d %d %d %d td id' . 

&u ( 0 ) , &u 1 1 1 , fcu [ 2 ] , &u [ 3 ] ,fcu[0] ,&u[l] ,tu[2] ,&u[3) ) , 

f [p] . numDashes [ i] =count ; 
for { j=0; j<8; j++) 

f [p] .dashes ( i ] [j]=(char) u [ 3 ] ; 

) 


for ( i=0 ; i<7 ; i++) { 

XtVaGetValues (FEI tem [FEBX] [i] 

strcpy ( s [ i 3 , t) ; 

f [p] ,xMin=atoi (s [0] ) ; 
f[p] .xMax=atoi (s [1] ) ; 
f [p] . xNum=atoi ( s [2 ] ) ; 
f [p] . xTick=atoi (s [3 ] ) ; 
stmcpy ( f [p] ■ xFormat , s [4 ] ,15) ; 
f [p] . borderL=atoi (s [ 5 ] ) ; 
f [pj .borderR=atoi (s [6] ) ; 


XtNstring, Set, NULL) ; 


for ( i=0 ; i<7 ; i++ ) { 

XtVaGetValues (FEItem[FEBY] 
strcpy (s [i] , t) ; 


[i] ( XtNstring, &t, NULL) ; 


f [p] .yMin=atoi ( s [0] ) ; 
ftp] .yMax=atoi (s 1 1 ] ) ; 
f [p] ,yNum=atoi (s [2] ) ; 
ftp] .yTick=atoi(s(3] ) ; 
strnepy ( f tp] .yFormat , s [ 4 ] , 15 ) ; 
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f [p] . borderT=atoi ( s [ 5 ] ) ; 
ftp] .borderB=atoi (s [6] ) ; 


for ( i=0 ; i<numCols ; i++) 

XtVaGetValues (FEItemtFEBVars ] [i] , XtNstate, fcuseCol [p] [i] , NULL) ; 
for (i=0; i<numTimeDep; i++) 

XtVaGetValues (FEItem[2] [ i ], XtNstate, fcuseCol [p] [MAXCOLS+i] , NULL) ; 


f or { i=0 ; i<8 ; i++ ) 

XtVaGetValues (FEItem[FEBOpts 3 [i] 


XtNstate, &f [p] . opt [ i ] , NULL) ; 


void RecalcFormat (p) 
int p; 

{ f [p] .xlnc={f[p] .width -f [p] .borderL-f Cp] .borderR) /f [p] .xNum; 
f [p] .ylnc= ( f [p] .height-f [p] .borderT-f [p] . borderB) / f [p] .yNum; 
ftp] . nGrid=f [ p ] .xNum+f [p] .yNum+2; 
ftp] . xLength=f [p] .xNum*f [p] .xlnc; 
ftp] .yLength=f [p] .yNum* ftp] .ylnc; 

/*PrintFormat (p> ; */ 

) 


/* 

void PrintFormat (p) 


int p; 


printf { " \n" 

) ; 

printf ( " 

width: 

printf ( " 

height : 

printf ( “ 

xlnc : 

printf C 

ylnc : 

printf ( * 

xLength 

printf ( • 

yLength 


} 

*/ 


%d\n" , f [p] .width) ; 
%d\n" , f [pj .height) ; 
%d\n" , f C p] .xlnc) ; 
%d\n* , ftp] .ylnc) ; 
%d\n" , f [p] . xLength) ? 
%d\n" , f [p] .yLength) ; 


void LoadFormat { ) 

{ 

FILE *fi; 

char s [64] , *t; 

int i, j , a=activeFormat ; 

XtVaGetValues (FEItemt 9] [0] , XtNstring ( it, NULL) ; 
strcpy (s, t) ; 


f i = fopen (s , "r ■ ) 
if (f i ! -NULL) { 
fgets ( t , 63 , f i) ; 
strcpy (s, t) ; 
if (strlen(s) >1) 
fgets (t, 63 , fi) ; 
strcpy (s, t) ; 
if (strlen(s) > 1 ) 
fgets (s , 63 , f i) ; 
strcpy (s, t) ; 
if (strlen (s) > 1 ) 


strncpy(f [a] . title, s, strlen (s) -1) ? 
strncpytf [a] . xFormat, s , strlen (s ) -1 ) ; 
strncpy (f [a] .yFormat, s, strlen (s) -1) ; 


for (i=0;i<8;i++) ( 

fscanf { f i , " %d* , &f [a] . color [ i] ) ; 
printf ( "%d - , f [a] .color [i] ) ; 

) 

printf ( * \n* ) ; 
for { i = 0 ; i<8 ; i + + ) { 
fscanf (fi, "%d",&j) ; 

f [a] . opt [ i ) = { ( j-=0) ? False : True ); 

fscanf (fi, "%d %d* , &f [a] . xNum, &f [a] .yNum) ; 

fscanf ( fi, *%e %e %e %e" , &f [a] .xMin, &f [a] .xMax, &f [a] .yMin, &f {a] .yMax) 
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fscanf (fii “%d %d- , &f (a] .xTick, la] .yTick) ; 

fS &m{ £ b;r5er^?? a ] d birderR,tf [a] -borderT, *f [a] .borders) ; 
fclose ( f i) ; 

SetFEValues (a) ; 


void SaveFormatO 

{ 

FILE *fo; 
char s [64] , *t; 

int i , a=act iveFormat ; 

XtVaGetValues ( FEI tem [ 9 ] [0] , XtNstring, &t,NULL) ; 
strcpy (s, t) ; 

if (s[0]!=NULL) { 

fo = fopen ( s , “w") ; . . 

forint f(fo f "%s\n",f[a] .title) , , 

fprintf (fo, ■%s\n%s\n- ,f [a] .xFormat, f [a] .yFormat) ; 

for ( i=0 ; i<8 ; i++ ) 

fprintf (fo, ■ %d *,f[a] .color[i]) ; 
fprintf ( fo , ‘ \n" ) ; 

for ( i=0 ; i<8 ; i++) 

fprintf (fo, "%d ■ , f [a] . opt [ 1 J ) ; 
fprintf ( fo, " \n" ) ; 

fprintf ( fo! "%f ' %f f [a^ •yM* x ) 

fprintf < fo, -%d %d\n* , f [a] .xTick, f [a] .yTrck) ; 

f.ftJS&.f W -borderT, f [a] -borders, ; 

fclose(fo) ; 


void FEPageCallback (w, client_data, call_data) 
Widget w; XtPointer client_data, call_data; 


( 


int i = (int) client_data; 
char s [64] ; 

switch (i) { 

case Q j 

if (activeFormat< (MAXPLOTS-1) ) { 

activeFormat++ ; 

SetFEValues (act iveFormat) ; 

sprintf (s, "Format #%d" # act iveFormat ) ; 
XtVaSetValues (FETop [FEBPage] , XtNlabel , s , NULL) , 

) 

break; 
case 1: 

if ( ac t iveFormat >0) { 

act iveFormat-- ; 

SetFEValues (act iveFormat) ; 
sprintf (s, -Format #%d" , act iveFormat) ; 
XtVaSetValues (FETop [FEBPage] , XtNlabel ,s, NULL) , 


} 

break; 


) 


void FECommandCal lback ( w , c 1 ient_data , cal 1 data) 
Widget w; XtPointer client_data, call_data, 
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int i = (int) client_data; 

switch (i) { 

case feLoad: 

LoadFormat { ) ; 
break ; 

case feSave: 

SaveFormat ( ) ; 
break; 

case feCancel : 

ShowHide2 {dummy , wFormat , True, 0) ; 
SetFEValues (activeFormat) ; 
break; 

case feAccept: 

GetFEValues (activeFormat) ; 
if (plotExists [activeFormat] ) { 

RecalcFormat (activeFormat) ; 
StatData (activeFormat) ; 
SetupGrid (activeFormat) ; 
SetupAxisLabels (activeFormat) ; 
SetupPoints (activeFormat) ; 
RedrawPlot2 (activeFormat) ; 

} 

break; 

case feRevert: 

SetFEValues (activePlot) ; 
break; 

} 

) 


xpGeneral . c 


*/ 

♦include "xplncs.h" 
♦include "xpDefs.h* 

♦include ’xpwidgets .h" 
♦include *xp commands . h* 
♦include *xpdata.h“ 
♦include "xpplots.h* 

/* $1: Global variables */ 

Screen ‘screen; 

Display ‘display; 

Boolean windowUp [ 48 ] ; 


GC gc [MAXGCS ] ; 

XColor exact , color [ 16 3 , rColor [16); 

Font f 6x10 , f 6x12 , f 6x13 , fHelvl8b; 

XFontStruct ‘f s6x!0 , *f s6x!2 , *fs6xl3 , *fsHelvl8b; 


XawTextPosition startPos=0 , endPos=0 ; 

char inName [ 6 4 ] = {"'} , dsName ( 64 ) = { ■ • ) , P f Name [ 6 4 ) = { - - 


) , plName [4] = { 




/* $2: Local variables - (none) */ 


/* $3: Global functions */ 


void LoadFonts ( ) ; 
void LoadColors ( ) ; 
void InitData ( ) ; 
void ProcessKey ( ) ; 
void Idle ( ) ; 


/* $4: Local functions - (none) */ 
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void LoadFonts { ) 

( 

f 6x10 
f 6x12 
f 6x13 
fHelvl8b 


XLoadFont (display, "6x10" ) 

XLoadFont (display, "6x12") 

XLoadFont (display, "6x13") . t * ifl-180- 

XLoadFont (display, ■ -* -helvetica-bold-r 18 180 


f ") 


f s6xl0 
f s6xl2 
f s6xl3 
fsHelvl8b 

} 


XQueryFont (display , f 6x10) ; 

XQueryFont (display, f 6x12) ; 
XQueryFont (display, f 6x13) 
XQueryFont (display , fHelvl8b) ; 


void LoadColors ( ) 

{ 

Colormap cmap; 

Status r; 

FILE *fi; 

char s [ 64 ] , t [ 64 ] , *u , *v; 
int i, count =0; 

cmap = XDefaultColormap( display, DefaultScreen (display) ) ; 


fi = fopen ( "XPlot .clr ’ , "r • ) ; 
if (fi!=NULL) { . , 

while ( fscanf ( f i, "%d %s' , Sa , s) . =EOF) { 
pr int f (" — %4d %s\n",i,s); 

lf rixAllocNamedColor (display , cmap, s , fcexact , fccolor [ i ] ) 


} 

} 

fclose ( f i) ; 


} 

r =XAl locNamedColor (display, cmap , " Black ] 
r=XAllocNamedColor (display, cmap, "White 
r=XAl locNamedColor (display , cmap , *Red" 
r=XAl locNamedColor (display , cmap , * Green 

r=XAllocNamedColor (display , cmap, "Blue 
r=XAllocNamedColor (display , cmap , " Yellow 


, &exact, fcrColor [Black) ) ; 

, &exact, trColor [White] ) ; 
,&exact,&rColor [Red] ) ; 

, &exact , &rColor [Green] ) ; 

, texact, &rColor [Blue] ) ; 
,texact,&rColor [Yellow] ) ; 


for (i= 0 ;i<( 2 *MAXCOLS) ;i++) 
lineColor [i]=color [i] .pixel; 


) 

void InitData ( ) 

{ 

int i , j , k; 

for ( k=0 ; k<MAXPLOTS ; k++ ) 

for ( j=0 ; j< (2*MAXCOLS) ? j++) 
useCol[k] (j]=False; 

for ( i=0 ; i<NUMcliString ; i++) 

strcpy (cliString [ i ] , mm ) ; 


for ( i=0 ; i<MAXPLOTS ; i++) 
plotExists [ i] =False ? 


) 

void XnitGCs { ) 

Window wind = XtWindow ( topLevel) ; 
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int i ; 

for (i=0; i<4; i++) 

gc [ i ] =XCreateGC (display , wind ,0,0); 

XSetForeground (display, gc[0] , rColor [Black] .pixel) ; 
XSetFont (display, gc[0] , £6x10) ; 

XSetForeground ( display , gc [ 1 ] , rColor [Black] .pixel) ; 
XSetFont (display, gc[l] , f6xl2) ; 

XSetForeground (display, gc [2 ] , rColor [Black] .pixel) ; 
XSetFont (display, gc [2] , fHelvl8b) ; 


) 

void DoKey ( w , event, params, numParams) 

Widget w; XEvent *event; String *params; Cardinal *numParams? 

XKeyEvent *e = (XKeyEvent*) event? 

} 

void Idle (w, event, params, numParams) 

Widget w; XEvent * event; String *params; Cardinal *numParams; 

{ 

) 

/* 


xpGeometry . c 


*/ 

♦include "xplncs.h" 

♦include “xpDefs.h* 

♦include "xpwidgets .h* 

♦include "xpgeneral . h" 

♦include "xpcommands .h* 

♦include ■ xp input. h" 

♦define pi 3.1415926535 

/* $1: Global variables - (none) */ 

/* $2: Local variables */ 

int activeGeometry=-l ; 
float geomVars[8]; 

static String MGInputLabels [ 6 ] [ 8 ] = { 

( "Web Thickness Internal Radius* , "Half-Angle" , "Epsilon* , 
•Fillet Radius*," ",* * , " " ) 

); 

static String MGOutputLabels [ 6 ) [ 8] = { 

{ "Initial Area*, "Final Area",* " > " / " " > " • " ' " ^ 

) ; 

/* $3: Global functions */ 

void MGCallback ( ) ; 
void ChsmgeMGLabels ( ) ; 
void ResizeGeometry ( ) ; 
void Expo seGeome try ( ) ; 

/* $4: Local functions */ 
void RecalcGeometry {) ; 
void Redr awGeome try ( ) ; 
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/* $5: Implementation ==== — -= — ======= 

void MGCallback(w, client_data, call_data) 

Widget w; XtPointer client_data, call__aata; 

* int n = (int) client_data; 
int m; 

char s[64],*t; 
switch (n) { 

C3 m = (int) XawToggleGetCurrent (MGItem[2 ] [0] ) ; 
if (m! =NULL) { 

activeGeometry=m-l ; 

RecalcGeometry 0 ; 

} 

break; 

Ca showHide2 (dummy, wMotorCon figuration, True, 0) ; 

break; 

Ca showHide2 (dummy, wMotorGeome try, False, 0) ; 

break; 


void ChangeMGLabel s (n) 
int n; 

{ 

int i ; 
n = 0; 

f °XtVaSetValues I MGI tem [ 0 ] [i] .XtNlabel.MGInputLabels [n] [i] , NULL) 

f °XtVaSetValues*MGItem[3 J [i] .XtNlabel.MGOutputLabels [n] [i] .NULL) 


void RecalcGeometry 0 
{ 

int i , j , k ; 
char s [64] , *t ; 
float re, ri, rr ,w,bp; 


for (i=0; i<8;i + +) ( . t \ . 

XtVaGetValues (MGItem[l] [ i] , XtNstring, &t,NULL) , 

strcpy (s , t) ; 

sscanf (s , *%e* , tgeomVars [ 1 J ) # 

} 

switch (activeGeometry) { 

case 0 : . , . • . 

re=geomVars [0] ; ri=geomVars [ 1] ; w=re-ri, 

bp=2 *ri*pi ; 
break; 

Ca re=geomVars[0] ; ri=geomVars [1] ; rr=geomVars [ 2 ] ; w=re-ri; 
bp= 2 *ri*pi+ 2 *rr*pi; 
break ; 
case 2 : 
break ; 

sprintf (s, *%.4f*,bp) ; 
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XtVaSetValues (MGItem[4] [0] , XtNstring , s , NULL) ; 
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) 


void Expo seGeome try (w, event , params , numParams ) 

Widget w; XEvent * event; String *params; Cardinal *numParams; 

XExposeEvent *myevent = (XExposeEvent* ) event; 
static int i=G; 


if ( (*myevent) .count==0) { 

i++ ; 

print f ( • : Expos eGeome try (%d) . . . \n" , i) ; 
Redr awGeome try { ) ; 

) 

} 


void ResizeGeometry (widget , event, parains , numParams) 

Widget widget; XEvent *event; String *params; Cardinal 

static Dimension oldWidth=360 , oldHeight=360 ; 

Dimension w,h; 


♦numParams ; 


XtVaGetValues (MGPict , XtNwidth , &w , XtNheight , &h , NULL} ; 
if { (w<=oldWidth)&&(h<=oldHeight) ) { 

pr int f ( “ : ResizeGeometry calling Redr awGeome try . . . \n ); 
RedrawGeometry ( } ; 


oldWidth=w; 

oldHeight=h; 


void RedrawGeometry ( ) 

( 

Dimension w,h; 

Window wind = XtWindow (MGPict ) ; 


printf ( - : RedrawGeometry . . .\n") ; 

XtVaGetValues (MGPict , XtNwidth, &w, XtNheight , &h, NULL) ; 
XClearWindow (display , wind) ; 


XSetLineAttributes (display, gc { 0 ] , 4 , LineSolid, CapButt , JoinRound) ; 

XDrawArc(display, wind, gc[0] ,4,4, w-8,h-8, 0,360*64) ; 

XSetLineAttributes (display , gc ( 0] , 0, LineSolid, CapButt, JoinRound) ; 


switch (activeGeometry) { 
case 0: 
break; 

} 


) 

/* 

xplnput . c 

*/ 

# include "xplncs.h" 
♦include "xpDefs.h* 

♦include “xpwidgets .h* 
♦include "xpgeneral . h" 
♦include ■xpmenus.h" 

♦define MAXVARS 44 
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♦define NDELX 0 
♦define TPI 1 

♦define PAM 2 

♦define NINERT 3 
♦define UNIT 4 
♦define AT 5 

♦define XP 6 

♦define XE 7 

♦define XG 8 

♦define GAMA 9 
♦define W 10 

♦define TIGN 11 
♦define TFREF 12 
♦define RUFSUR 13 
♦define DDRG 14 
♦define DDHC 15 
♦define FKPR 16 
♦define ROPR 17 
♦define CPR 18 
♦define TOREF 19 
♦define SIGP 20 
♦define TPSCRI 21 
♦define RREF 22 
♦define PREF 23 
♦define BREXP 24 
♦define EBC 25 
♦define EBEX 26 
♦define DE 27 

♦define CM 28 

♦define ALFAD 29 
♦define EROAT 30 
♦define EROEXP 31 
♦define TMAX 32 
♦define DELTAT 33 
♦define TPRINT 34 
♦define PZONE 35 
♦define LAMBDA 36 
♦define NIGTAB 37 
♦define NAPDVX 38 
♦define NDATA 39 
♦define NPPR 40 
♦define ITABLE 41 
♦define GTABLE 42 
♦define BTABLE 43 


/* $1: Global Variables */ 

char varName [ MAXVARS ] [32] , varDesc [MAXVARS] [80] ; 
char varNumb [MAXVARS] [12] , var Inf o [MAXVARS] [800] 
float varValue [MAXVARS] ; r _ AW _. 

float iTable [ 30 ] [2 ] , gTable [ 30] [8] ,bTable[30] [2] 

/* $2: Local Variables */ 

int numVars , inputlndex=0 ; r-iAAirArt1 

char *IEStrings [100] , IEEntnes [100] [20]; 

/* $3: Global Functions V 


void LoadlnputData { ) ; 
void IECallback ( ) ; 
void IEListCallback ( ) ; 
void Wr i telnputFile ( ) ; 
void AcceptlnputEntry { ) ; 

/* $4: Local Functions */ 

void Get InputValues ( ) ; 
void Set InputValues ( ) ; 
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void IncDecInput () ; 

/* $5: Implementation ===========================«==========”===— */ 

void LoadlnputData { ) 

{ 

FILE* fi; 
int i,j=0; 
char s [64 ] ; 

fi=fopen( "XPlot.var" , "r" ) ; 
do { 

i=f scanf (f i , ■%*[ \t\n] " > ; 

if { i ! =EOF) { 

i = f s can f(fi, •\ , '%*d/\n") ? 
i=f scanf (fi, - %[ A /] /\n‘ , varName [ j 3 ) ; 
i=fscanf (fi, "%[ A /1 /\n“ , varDesc [ j 3 ) ; 
i= f scanf (fi, *% [ A /] /\n w , varNumbt j 3 ) ; 
i=fscanf (fi, ■% ( A \ ' 3 \ ’ ’ , varlnfo [ j 3 ) ; 

} 

j++; 

} while ( ( i ! =EOF) && { i ! =0) && ( j<MAXVARS) ); 

numVars= j -1 ; 


for ( j=0 ; j<numVars ; j++) 

i = sscanf ( varNumb [ j ] , * %e " , &varValue [ j 3 ) ; 
printf ( * \n" ) ; 


for ( i =0 ; i<=BTABLE ; i++ ) { 

if ( i cITABLE) 

sprintf ( IEEntries ( i ] , " %-6s 


%10.4f ■ ,varName[i3 ,varValue[i] ) ; 


sprintf (IEEntries [i] , ■ %c%s " # varName [ i ] [ 0 ] , * TABLE “ ) ; 

IEStrings [ i 3 = IEEntries [ i 3 ; 
printf ("%2d %s\n*,i, IEEntries [ i 3 ) ; 


} 

printf ( * \n* ) ; 


XawLis tChange ( lEList , IEStrings , BTABLE+1 , 0 , False) ; 

} 


void DoInputO 
{ 

int i,j,k; 

char s [ 64 3 , * t , *u, *v; 

ShowHide2 (dummy, 1000, True, 0) ; 

} 

void AcceptlnputEntry (w, event, params , numParams } 

Widget w; XEvent *event; String *params ; Cardinal ‘numParams; 

{ 

if (numParams >0 ) { 

if (params [0] [0 ]=='+' ) 

IncDecInput (True) ; 
if (params [0] [ 0 ] == ' - ' ) 

IncDecInput (False) ? 

} 

} 

void GetlnputValues ( i ) 
int i ; 

{ 

int j,k; 

char s [25003 , *t; 

if ( (i>=0) && (i<ITABLE) ) { 
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XtVaGetValues (IEItem[3] [0] , XtNstring , &t , NULL) ; 
strcpy (s , t) ; 

sscanf ( t , " %e " , kvarValue [ i ] ) ; . r-i\ 

sprintf (IEEntries [i] , *%-6s = %10 . 4f " , varName [ i] , varValue [ i ] ) 

XawListUnhighlight ( IEList} ; 


} 


void SetlnputValues (i) 
int i; 

( 

int j,k; 

char s[2500] , t [ 8 0 ] ; 

if { { i> = 0 ) { i<ITABLE) ) { t __ .-I 

XtVaSetValues ( IEI tem [ 0 ] 1 0 ] , XtNstring , varName i , NULL 
XtVaSetValues ( IEItem[ 1] [0] , XtNstring, varDesc i , NULL 
XtVaSetValues (IEI tem [2] [0] , XtNstring, varlnfoli] , NULL) 
sprintf (s, "%.2f" ,varValue[i] ) ; ( . 

XtVaSetValues (IEItem[33 [0] XtNstring, s,NULL) , 
XawListHighlight { IEList , i) ; 

} 

if ( i = -ITABLE ) { . _ r ' i ktttt t \ 

XtVaSetValues {IEI tem [0] [0] , XtNstring, varName i 1 .NULL) 
XtVaSetValues { IEI tem [ 1 1 [ 0 ) , XtNstring , varDesc i , NULL 
XtVaSetValues ( IEI tem [ 2 1 [ 0 ] , XtNstring , var Inf o [ i ] , NULL ) 
strcpy (s, " " ) ; 
k=varValue [NIGTAB] ; 

f °sprintfM t^"%lo! 4f%10 . 4f \n" , iTable [ j ] [0] , iTablel j ] [ 1] ) 
strcat(s, t) ; 


} 


XtVaSetValues (IEI tem[ 3 ] [0] , XtNstring, s 
XawListHighlight (IEList, i) ; 


, NULL) 


if ( i ==GTABLE ) { , „ vtttt t \ 

XtVaSetValues ( IEI tem [ 0 ] [ 0 ] , XtNs tr mg , varName i , NULL 
XtVaSetValues ( IEI tem [ 1 1 [0] , XtNstring, varDesc i -NULL 
XtVaSetValues (IEItem[2] [0] , XtNstring, var Info[i] .NULL) 
strcpy (s , *"); 
k=varValue [NAPDVX] ; 
strcpy (s, ■ •) ; 
for ( j=0; j<k; j++) { 

sDrintf ( t , •%10.4f%10.4f%10.4f%10.4f" , 

gTable [ j ] [0] ,gTable[j] [1] ,gTable[j] [2] ,gTable[:)} [3]) 

s treat {s , t ) ; 
if (gTable [ j ] [4 ] >0 ) { 

sprintf (t, "%10.4f%10.4f%10.4f%10.4f", 

gTable ( j ] [4] , gTable [ j] [5] ,gTable{ j] [6] , gTable [ 3 ] [7] ) 

streat (s, t) ; 

} 

streat (s , " \n" ) ; 

XtVaSetValues (IEI tem l 3] [0] , XtNstring, s , NULL) ; 

XtVaSetValues (IEItem[3] [0] , XtNstring , s , NULL) ; 
XawListHighlight ( IEList , i ) ; 


) 


if { i ==BTABLE ) { „ r *i mttt t ^ 

XtVaSetValues (IEItemtO) [0] .XtNstring, varName [ 1 ] .NULL) 
XtVaSetValues ( IEI tem [ 1 ] 1 0 ) , XtNs tr ing , varDesc i , OTJLL 
XtVaSetValues ( IEItem[2 ] [0] , XtNstring, varlnfoli] .NULL) 


strcpy (s, " " ) ; 
k=varValue [NIGTAB] ; 
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} 


for ( j=0; j<k; j++) { 

sprint f ( t , *%10.4f%10 
strcat (s , t) ; 


. 4f \n" , bTable [ j ] [0] ,bTable[ j] [1] ) ; 


XtVaSetValues { IEItem[ 3 ] 1 0 ] , XtNstring , s , NULL) ; 
XtVaSetValues { IEI tem[ 3 ] [ 0 ] , XtNstring , s , NULL) ; 
XawListHighlight ( IEList , i ) ; 


void IECallback (w, client_data, call_data) 

Widget w; XtPointer client_data, call_data; 

* int n = (int) client_data; 
char s [ 128] , *t; 

switch (n) { 

case 0: /* Start | Restart */ 
inputlndex=0 ; 

SetlnputValues { inputlndex) ; 
break; 

case 1: /* Insert */ 
break; 

case 2: /* Prev */ 

IncDecInput (False) ; 
break ; 

case 3: /* Next */ 

IncDecInput (True) ; 
break ; 

case 4: /* Write */ t 

XtVaGetValues ( IEItem( 4 ] [ 0 ) , XtNstring, &t , NULL) ; 
strcpy(s f t); 

if (strcmp(s, - •) 1=0) WritelnputFile (s) ; 
break; 

case 5: /* Cancel V 

ShowHide2 (dummy ,wlnput, False, 0) ; 
input lndex=0 ; 

SetlnputValues (inputlndex) ; 
break; 

} 

} 

void IEListCallback (w, client_data, call_data) 

Widget w; XtPointer client_data, call_data; 

( XawListRetumStruct *r = (XawListRetumStruct* ) call.data; 
char s [64] ; 
int n; 

n = (*r) .list_index; B 

printf ("Element [%02d] selected from IEList\n ,n); 

GetlnputValues (inputlndex) ; 
inputlndex=n; 

SetlnputValues (inputlndex) ; 

} 

void IncDecInput ( inc ) 

Boolean inc; 

{ 

char s [ 1024 ] , *t ; 

if ( ( inc) && (inputlndex< (BTABLE) ) ) { 

GetlnputValues (inputlndex) ; 

input Index* + ; 

SetlnputValues ( inputlndex) ; 

} 

if ( ( i inc) Set ( inputlndex>0) ) { 
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GetlnputValues { inputlndex) ; 
inputlndex--; 

SetlnputValues (inputlndex) ; 


void WritelnputFile (s) 
char *s; 

{ 

FILE *f ; 
int i,j,k; 

f = fopen(s, "w* ) ; 

fprintf (f," &NAME\n“); 
fprintf ( f , " TMAX = 

fprintf (f," LAMBDA - 

fprintf ( f , " TPRINT = 

fprintf ( f , * NPUNCH = 

fprintf (f,‘ NBL 
fprintf (f, # NDELX = 

fprintf ( f , " TPI = 

fprintf (f,“ PAM ” 

fprintf ( f , * UNIT 

fprintf ( f / “ NINERT 

fprintf ( f , * NPNPXT 

fprintf (f," AT 

fprintf (f," XP 

fprintf (f,* XIP 

fprintf (f," XG 

fprintf { f , * XE 

fprintf(f," GAMA 

fprintf (f, a TFREF 

fprintf (f,“ TIGN 

fprintf ( f , " RUFSUR 

fprintf (f,“ DDRG 

fprintf (f , " DDHC 

fprintf ( f , " API AT 

fprintf (fr" AFTOS 

fprintf ( f , * FKPR 

fprintf (f,* ROPR 

fprintf ( f , ’ CPR 

fprintf (f," TOREF 

fprintf (f,“ SIGP 

fprintf ( f , • TPSCRI 

fprintf ( f , " RREF 

fprintf (f,* PREF 

fprintf (f," BREXP 

fprintf (f r “ EBC 

fprintf (f," EBEX 

fprintf ( f , * DE 

fprintf {f , * CM 

fprintf (f," ALFAD 

fprintf (f," EROAT 

fprintf (f," EROEXP 

fprintf (f,* NIGTAB 

fprintf (f , * NAPDVX 

fprintf (f,* NDATA 

fprintf ( f , “ NPPR 

fprintf (f,* PZONE 

fprintf (f,* NRESRT 

fprintf ( f , * DELTAT 

fprintf (f# * PCI 

fprintf <f," GAMAN 

fprintf (f,* PISUBK 

fprintf (f,* DELTTF 

fprintf (f, * DFSDT 

fprintf (f , p TUN 

fprintf (f," FCRIT 


%10 . 4f \n w , varValue [ 0}); 
%10 . 4f \n" , varValue [ 1]); 
%10.4f\n", varValue [ 2] ) ; 
%10 . 4f \n" , varValue[ 3] ) ; 
%10 . 4f \n" , varValue[ 4] ) ; 
%10 . 4f \n 9 , varValue[ 5] ) ; 
%10 . 4f \n" , varValue [ 6] ) ; 
%10 ,4f\n* .varValue [ 7]); 
%10 . 4f \n" , varValue [ 8] ) ; 
%10.4f\n" , varValue [ 9]); 
%10 . 4f\n" , varValue [10] ) ; 
%10 . 4f \n* f varValue [11] ) ; 
%10 . 4f \n # , varValue [12] ) ; 
%10 . 4f\n* , varValue [13] ) ; 
%10 . 4f \n" , varValue [14] ) ; 
%10 . 4f \n" , varValue [15] ) ; 
%10 - 4f \n“ , varValue [16] ) ; 
%10 . 4f \n* , varValue[17] ) ; 
%10 . 4f \n* , varValue[l8] ) ; 
%10 . 4f\n* , varValue [19] ) ; 
%10 . 4f \n a , varValue[20] ) ; 
%10 . 4f \n a , varValue [21] ) ; 
% 10. 4f\n a , varValue [22] ) ; 
%10 . 4f \n a , varValue [23 ] ) ; 
%10 . 4f \n* , varValue[24] ) ; 
%10 . 4f\n a , varValue[25] ) ; 
%10 . 4f \n" , varValue [26] ) i 
%10.4f\n“ , varValue [27] ) ; 
%10 . 4f \n" , varValue[28] ) ; 
%10 . 4f \n* , varValue [29] ) 
%10 . 4f\n* , varValue [30] ) 
%10 . 4f \n* , varValue [31] ) 
%10 . 4f \n" , varValue [32 ] ) 
%10 . 4f \n“ , varValue [33] ) 
%10 . 4f \n* , varValue [34] ) 
%10 . 4f\n" , varValue [35] ) 
% 1 0 . 4 f \n " , varValue [36]) 
%10.4f\n“ , varValue [37] ) 

= %10 . 4f \n“ , varValue [38] ) 

= % 1 0 . 4 f \n a , varValue [39]) 

= % 10 . 4f \n a , varValue[40] ) 

= % 10 . 4f\n a , varValue [41] ) 
= %10 . 4f\n" , varValue [42 ] ) 
= %10 . 4f\n a , varValue [43 ] ) 
= %10.4f\n a ,yarValue[44] ) 
= %10 . 4f\n a , varValue [45] ) 
= % 10 . 4f\n“ , varValue [46] ) 
= %10.4f\n", varValue [47] ) 
= %10.4f\n", varValue [48 ] ) 
= %10 . 4f\n a , varValue [49] ) 
= %10 . 4f \n a , varValue [50] ) 
= %10 . 4f \n" , varValue[51] ) 
= % 10 . 4f\n" , varValue [52] ) 
= % 10 . 4f\n a , varValue [53] ) 
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384 
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10 
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fprintf ( f , * 
fprintf ( f , * 
fprintf ( f , ' 
fprintf ( f , * 
fprintf { f , ■ 
fprintf { f , " 


POPEN = %10 . 4f \n* , varValue [54] ) ; 
DELFAC - %10.4f\n", varValue [55] ) ; 
CHC = % 10 . 4f\n" , varValue [56] ) ; 
BF1 = %10 . 4f\n" , varValue [ 57] ) ; 
DELTFF = %10 . 4f\n" , varValue [58 ]) ; 
&END\n* ) ; 


for ( i=0 ; i<varValue [NIGTAB] ; i++) , 

fprintf (f , ■%10.4f%10.4f\n", iTable [ jl (0 3 , iTable [ 3 3 [ 13 ) ; 


for (i=0; i<varValue [NAPDVX] ; i + +) ( 

forint f (f "%10.4f%10.4f%10.4f%10.4f" , 

^ P gTable [ j ] [0] , gTable [j ] [ 1] , gTable [j-].[2] , gTable t ) 

if (gTable [ j ] [ 4 ] ! =0) fprintf ( f, "%10 . 4f" , gTable [ 3 3 [4] ) , 
if (gTable [ j ] [5] ! =0) fprintf ( f, ’%10 . 4f , gTable [ j ] [5 }) ; 
fprintf (f , ‘ \n" ) ; 

for (i=0; i<varValue [NDATA] ; i++) p ., rdll . 

fprintf if, *%10 . 4f%10 . 4f \n" , bTable[jl (03 , bTable [ 3 ] (13) ; 

fclose(f ) ; 


) 


/* 

xplot . c 

*/ 


♦include ■xplncs.h" 
♦include •xpDefs.h* 


♦include 

♦include 

♦include 

♦include 

♦include 

♦include 

♦include 

♦include 

♦include 


"xpwidgets . h* 
■xpgeneral .h" 
•xpmenus .h" 

■ xpcommands . h • 
"xpinput .h" 
"xpdata.h" 
"xpformats .h* 
"xpplots . h" 

" xpgeometry . h " 


void main{argc,argv) 

int argc; char **argv; 

{ 

F I LE * f i ; 

static XLActionsRec xplotActions ( 3 = { 
("idle*, Idle}, 

{"doKey" ,DoKey} , 

( "redrawPlot" ,RedrawPlot) , 

{ " saveCommand" , SaveCommand} , 

{ ■ recallCommand" f Recall Command) , 

{ "resizePlot* ,ResizePlot} , 

{ " resizeGeometry ■ , ResizeGeometry) , 

{• expos eGeome try " , Expos eGeome try} , 

{ • accept InputEntry* , Accept InputEntry) , 

}; 

XtAppContext app_ context; 

int i ; 


fi = fopenCXPlot.dat*, "w") ; „ T _ 

topLevel = XtVaAppInitialize (&app_context , ’XPlot ,NULL,0, 
aargc , argv , NULL , NULL) ; 

X t App AddAc t i ons (app_context,xplotActions,XtNumber (xplotActions) } 


CreateWidgets () ; 
XtRealizeWidget (topLevel) ; 
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screen =XtScreen(topLevel) ; 
display=XtDisplay ( topLevel ) ; 

LoadFonts ( ) ; 

LoadColors ( ) ; 

InitData ( ) ; 

InitGCs ( ) ; 

for (i = 0; i<MAXPLOTS; i>+) 
InitFormat (i) ; 
LoadlnputData ( ) ; 

SetFEValues (0) ; 
SetStatusLabel ( ) ; 
ChangeMGLabe 1 s ( 0 ) ; 


printf ( * \n* ) ; 

XtlnstallAccelerators ( rootPane , cliText) ; 
XtlnstallAccelerators (cliLabel, cliText) , 

XtAppMainLoop (app_context) ,* 


) 


/* 

xpMenus . c 

*/ 

♦include "xplncs.h" 

♦include *xpDefs.h“ 

♦ include -xpwidgets . h* 

♦ include “xpgeneral .h* 

♦ include " xpcomnvands . h • 

♦include "xpdata.h" 

♦ include "xpformats . h' 

♦include -xpplots.h* 

/* $i; Global variables - (none) */ 
/* $2: Local variables - (none) */ 
/* $3: Global functions */ 


void MenuSelect ( ) ; 
void CoromandSelect ( ) ; 
void VerbCallback ( ) ; 
void SFCallback ( ) ; 
void SFGetFile ( ) ; 
void Selectstring 0 ; 
void ShowHide2 () ; 
void Loadlnput ( ) ; 
void FileCallback ( ) ; 
void ExecCallback ( ) ; 


f* $4: Local functions */ 
void SetStatusLabel ( ) ; 
void RunExec ( ) ; 


/* $5: Implementation ========== 


5id MenuSelect (w, client_data, call_data) 
widaet w; XtPointer client_data , call_data; 


int menuData = (int) client_data; 
int menuNum, itemNum; 
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char s [64] , *t ; 

FILE * f ; 

menuNum = menuData/HIWORD; 
itemNum = menuData%HIWORD; 


switch (menuNum) { 
case IDMenuFile: 
swi tch { itemNum) { 
case 0 : 

ShowHide2 (dummy, wFile, True, 0) ; 
break; 
case 1: 

XtVaGetValues (menuText , XtNstring , tt , NULL) ; 


strcpy (s , t) ; 
f = fopen(s, "r* ) ; 
if (f!=NULL) 

XtVaSetValues ( f ileText , XtNstring , s , XtNtype, XawAscnFile , NULL) 


break; 
case 2 : 

ShowHide2 (dummy, wFile, True, 0) ; 
break; 
case 3 : 
break; 
case 4 : 

XtVaGetValues (menuText, XtNstring, &t, NULL) ; 
strcpy (s , t) ; 
f = fopen (s , "w" ) ; 

XtVaGetValues { f ileText , XtNstring, tt, NULL) ; 
fprintf ( f , " %s " , t) ; 
fclose(f); 
case IDmFileQuit: 
exit (0) ; 


) 

break ; 

case IDMenuInpu : 

switch (itemNum) { 
case IDmlnpuSet: 

SelectString(O) ; 
break; 

case IDmlnpuLoad: 

Loadlnput ( ) ; 
break ; 

case IDmlnpuLoadSF : 
SelectString(O) ; 

Loadlnput ( ) ; 
break; 
case 4 ; 

ShowHide2 (dummy, wlnput, True, 0) ; 
break ; 


} 

break; 

case IDMenuExec : 

switch (itemNum) { 
case IDmExecSetup : 

ShowHide2 (dummy, wExecute, True, 0) ; 
break; 

case IDmExecRun : 

RunExec ( ) ; 


break; 


} 

break; 

case IDMenuData: 
switch (itemNum) { 
case IDmDataSet: 
Selectstring (1) ; 
break; 

case IDmDataLoad: 
LoadData (dsName) ; 
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break; 

case IDmDataLoadSF : 

Selectstring (1) ; . 

LoadData (dsName) ; 
break; 

case IDmDataView: 

ShowHide2 ( dummy , wData , True , 0 ) ; 
break; 

case IDmDataPrint : 

PrintDataCofile.dat") ; 
break; 

} 

break; 

case IDMenuForm: 
switch ( itemNum) { 

case IDmFormSet: 

Selectstring (2) ; 
break; 

case IDmFormLoadSF : 

Selectstring (2) ; 
break; 

} 

break; 

case IDMenuPlot ; 

switch (itemNum) { 

case IDmPlotSet: 

Selectstring (3) ; 
break; 

case IDmPlotlnc : 

IncDecPlot (True) ; 
break; 

case IDmPlotDec; 

IncDecPlot (False) ; 
break; 


} 

break ; 

case IDMenuWind: 

switch (itemNum) { 

case 0: ^ A . 

ShowHide2 ( dummy , wWindowManager , True , 0 ) ; 

break ; 
case 1: 

ShowHide2 (dummy ,wFile , True , 0) ; 
break; 
case 2: 

ShowHide2 (dummy, wlnput , True, 0) ; 
break; 
case 3: 

ShowHide2 (dummy, wExecute, True, 0) ; 
break ; 
case 4: 

ShowHide2 (dummy, wData, True, 0) ; 
break ; 
case 5; 

ShowHide2 (dummy, wFormat, True, 0) ; 
break; 
case 6: 

ShowHide2 (dummy, wFormatAxes , True, 0) ; 

break ? 

case 7 ; m , 

ShowHide2 ( dummy, wForma tL me s, True, 0) ; 

break; 

case 8; . m A . 

ShowHide2 (dummy, wFormatCurveFits, True, 0) 
break ; 
case 9: 

ShowHide2 (dummy , wPlotManager , True, 0 ) ; 
break; 
case 10 : 
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ShowHide2 (dummy, wMotorGeometry, True, 0) ; 
break; 
case 11: 

ShowHide2 (dummy , wMotorConf igurat ion, True, 0) 
break; 

} 

break; 

) 

} 

void CommandSelect (w, client_data, call^data) 

Widget w; XtPointer client_data, call_data; 

int commNum = (int) client_data; 
char s [64] ; 

static Boolean up(l] = {False}; 

switch (commNum) 

{ 

case 0: 

ShowHide2 (dummy , wWindowManager , True, 0) ; 
break ; 
case 1 : 

ShowHide2 (dummy , wFile, True , 0) ; 
break; 
case 2 : 

ShowHide2 (dummy, wlnput, True, 0) ; 
break; 
case 3: 

ShowHide2 (dummy, wExecute, True, 0) ; 
break ; 
case 4 : 

ShowHide2 ( dummy , wData , True , 0 ) ; 
break; 
case 5: 

ShowHide2 (dummy, wFormat, True, 0) ; 
break; 
case 6 : 

ShowHide2 ( dummy, wFormatAxes, True, 0) ; 
break; 
case 7 ; 

ShowHide2 (dummy , wFormatLines , True, 0) ; 
break; 
case 8: 

ShowHide2 (dummy , wFormatCurveF its , True , 0) ; 
break; 
case 9: 

ShowHide2 (dummy, wPlotManager , True, 0) ; 
break; 
case 10: 

ShowHide2 (dummy , wHo tor Geometry , True, 0) ; 
break; 
case 11: 

ShowHide2 (dummy, wMotorConf igurat ion, True, 0) ? 
break; 

) 

} 

void ShowHide2 (w, id, show, code ) 

Widget w; int id; Boolean show; int code; 

{ 

char s [ 64] ; 

if { ( id>=0 ) && ( id<8 ) && (plotExists [ id] ) ) { 

if (show) { 

windowUp [ id] =True ; 

XtPopup (plotShell [ id] ,XtGrabNone) ; 

} 

else ( 
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windowUp [ id] =False ; 

XtPopdown (plotShell [id] ) ; 

) 

} 

switch (id) { 

case wWindowManager : 
if ( ! windowUp [ id] ) { 

strcpy(s, "Hide Window Manager*) ; 

XtPopup (commands he 11 , XtGrabNone) ; 

} 

else ( 

s trcpy ( s , * Show Window Manager ) ; 

XtPopdown ( commandShell ) ; 

XtVaSetValues (c [0] , XtNlabel , s,NULL) ; 
XtVaSetValues (mEntry [ 6 ] [ 0 ] , XtNlabel , s / NULL ) , 

windowUp [ id] =! windowUp [id] ; 

break; 
case wFile: 

if (! windowUp [ id] ) { 

strcpy ( s , "Hide File Editor*) ; 

XtPopup ( f ileShell / XtGrabNone) ; 

) 

else t 

strcpy { s , " Show File Editor ); 

XtPopdown ( fileShell) ; 

XtVaSetValues ( c [ 1 ] , XtNlabel , s , NUbL ) ; 
XtVaSetValues (mEntry [6] [ 1 ] r XtNlabel f s , NULL) ; 
windowUp [ id] = ! windowUp [id] ; 
break; 

case wlnput: 

if ( i windowUp [ id] ) { 

strcpy (s , "Hide Input File Editor ); 

XtPopup ( IEShell , XtGrabNone) ; 

} 

else { . . 

strcpy (s , * Show Input File Editor ); 

XtPopdown ( IEShell ) ; 

XtVaSetValues (c [2] f XtNlabel , s f NULL) ; 
XtVaSetValues (mEntry ( 6] [2] , XtNlabel, s, NULL) ; 
windowUp [ id] = [windowUp [ id] ; 
break ; 

case wExecute: 

if (! windowUp ( id] ) { 

strcpy {s , "Hide Execution Setup*); 

XtPopup (exec Shell , XtGrabNone) ; 

} 

6156 t _ . # T 

strcpy ( s , "Show Execution Setup ) ; 

XtPopdown (execShell) ; 

XtVaSetValues ( c [ 3 ] , XtNlabel , s , N^L) • 
XtVaSetValues (mEntry [6] [3] , XtNlabel, s, NULL) ; 
windowUp [ id] = ! windowUp [ id] ; 
break; 
case wData: 

if ( [windowUp [id] ) { 

strcpy (s # "Hide Current Dataset*); 

XtPopup (da t a She 11 , XtGrabNone) ; 

) 

else { 

strcpy (s, * Show Current Dataset ) ; 

XtPopdown (dataShell) ; 

XtVaSetValues ( c ( 4 ] , XtNlabel , s , NULL) ; 
XtVaSetValues (mEntry [6] |4] .XtNlabel, s, NULL) ; 
windowUp(id) = iwindowUptid] ; 
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break; 

case wFormat ; 

if ( l windowUp [id] ) { 

strcpy (s, "Hide Format*}; 

XtPopup (FEShell , XtGrabNone) ; 

} 

else { 

strcpy ( s , " Show Format"); 

XtPopdown (FEShell ) ; 

XtVaSetValues (c[5] , XtNlabel , s , NULL) ; 
XtVaSetValues (mEntry [6 ] [5] , XtNlabel, s, NULL) ; 
windowUp [ id) =! windowUp [ id] ; 

break; 

case wFormatAxes: 
break; 

case wFormatLines : 

if (! windowUp [id] ) { 

strcpy (s , "Hide Format Lines"); 

XtPopup (DAShell, XtGrabNone ) ; 

} 

else { _ 

strcpy (s , "Show Format Lines"); 

XtPopdown (DAShell) ; 

XtVaSetValues <c [7 ] , XtNlabel , s , NULL) ; 
XtVaSetValues {mEntry [6] [7] , XtNlabel , s , NULL) ; 
windowUp [ id] =! windowUp [id] ; 
break; 

case wFormatCurveFits : 
if (! windowUp [ id] ) { 

strcpy (s, "Hide Format Curve Fits"); 

XtPopup (CFShell, XtGrabNone) ; 

} 

else { 

strcpy (s, "Show Format Curve Fits"); 
XtPopdown (CFShell) ; 

XtVaSetValues (c [8] , XtNlabel , s , NULL) ; 
XtVaSetValues (mEntry [6 ] [8] , XtNlabel , s, NULL) ; 
windowUp [ id] = i windowUp [ id] ; 
break; 

case wPlotManager : 

if (! windowUp [ id] ) { 

strcpy(s, "Hide Plot Manager"); 

XtPopup (PMShell, XtGrabNone) ; 

} 

else { 

strcpy (s, "Show Plot Manager"); 

XtPopdown (PMShell) ; 

XtVaSetValues (c [9] , XtNlabel , s , NULL) ; 
XtVaSetValues (mEntry [6] [ 9 ] , XtNlabel , s , NULL) ; 
windowUp! id] =! windowUp [ id] ; 
break; 

case wMotorGeometry : 
if (! windowUp [ id] ) { 

strcpy (s, "Hide Motor Geometry"); 

XtPopup (MGShell, XtGrabNone) ; 

) 

else { 

strcpy (s, "Show Motor Geometry"); 

XtPopdown (MGShell ) ; 

XtVaSetValues (c [ 10] , XtNlabel, s ,NULL) ; 
XtVaSetValues (mEntry [ 6 ] [ 10 ] , XtNlabel , s , NULL) 
windowUp [ id] =! windowUp [ id] ; 
break; 

case wMotorConf iguration: 
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if { !windowUp[id] ) ( . 

strcpy(s, 'Hide Motor Configuration ); 
XtPopup (MGPictShell , XtGrabNone) ; 

) 

€ 1 S 6 { 

strcpy (s, "Show Motor Configuration"); 
XtPopdown {MGPictShell) ; 


XtVaSetValues ( c [ 1 1 ] , XtNlabel , s , NULL) ; 
XtVaSetValues (mEntry [6] [11] , XtNlabel , s , NULL) 

windowUp [ id] = ! windowUp [id] ; 
break; 


) 


void VerbCallback (w, client_data, call data) 

Widget w; XtPointer client_data, call_data, 

{ XawL i s tRe turns t rue t *r = (XawListReturnStruct*) 

strcpytcli, (*r) .string) ; 

XtVaSetValues (cliText.XtNstring, cli.NULL) , 
XawTextSetlnsertionPoint (cliText , strlen (cli) ) > 


raid SFCallback (w, client_data, call_data) 
widaet w; XtPointer client_data, call_aata; 


int action = (int) client_data; 
char filename [32 ] ,s[64] ; 

String f? 


f * XawDialogGetValueString (sf Dialog) ; 
strcpy { filename , f ) ; 

XtPopdown (sfGetFileShell) ; 


switch (action) { 
case SFCommCancel : 
break; 

case SFCommLoad: 

sprintf (s , "The filename selected was: %s\n 
break; 

case SFCommSave: 
break; 

} 

} 


void SFGetFileO 

XtPopup (sfGetFileShell , XtGrabNone) ; 

) 

*/ 


void Selectstring (code) 
int code; 

( 

char s [64 ] ,*ss; 
int i ; 

XtVaGetValues (menuText.XtNstring, tss.NULL) ; 

strncpy ( s , ss , 63 ) ; 
s [ 63 ] = ' \0 ' ; 
if (s[0] ! = ' \0 ' ) { 

switch (code) { 
case 0: 


call_data; 


, filename) 
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strcpy (inName, s) ; 
break; 
case 1 : 

strcpy (dsName, s) ; 
break; 
case 2 : 

strcpy ( p f Name , s ) ; 
break; 
case 3 : 

i=atoi { s ) ; 

if ( (i>=0}fit&{i<=7) && (plotExists [ i ] ) ) { 

activePlot=i 
sprintf {plName , • %d" , i > ; 

> 

break; 

} 

SetStatusLabel ( ) ; 

} 

) 

void SetStatusLabel ( ) 

{ 

char s [256] ; 

strcpy (s , " I : <" ) ; 
s treat (s , inName) ; 
streat (s , •> D: < m ) ; 
s treat ( s , dsName) ; 
streat (s, ■> F: <* ) ; 
streat (s,pfName) ; 
streat (s, *> P: <" ) ; 
streat (s , plName) ; 
streat (s, *>■ ) ; 

XtVaSetValues (statusLabel ,XtNlabel , s,NULL) ; 

) 

void FileCallback (w, client_data ( call_data) 

Widget w; XtPointer client_data # call_data; 

{ 

int n = (int) client_data; 

FILE *f ; 
char s [ 128] , *t ; 

XtVaGetValues (fileEntry [3] , XtNstring , &t , NULL) ; 
strcpy (s ; t) ; 

XtVaGetValues ( f ileText , XtNstring , &t , NULL) ; 

switch (n) { 
case 0 : 

ShowHide2 (dummy, wFile, False, 0) ; 
break; 
case 1: 

f * fopen (s , ) ; 
fprintf ( f , “%s" , t) ; 
fclose ( f ) ; 
break; 
case 2: 

f = fopen (s, “r ■ ) ; 
if ( f ! =NULL) 

XtVaSetValues ( f i leText , XtNs tring , s , XtNtype , XawAsc i iFi le , NULL ) 
break; 

} 


void LoadlnputO 

{ 
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FILE *fi; 
char s [12 8] , *t; 

XtVaGetValues ( f ileEntry [3 ] , XtNs tr ing , &t , NULL) ; 
strcpy (s , t) ; 


f i = fopen (s, "r" ) ; 

1 f XtVaSetValues ( f ileText , XtNs tr ing , s , XtNtype, XawAsciiFile , NULL) ; 
} 


void ExecCallback (w, client_data, call_data) 
Widget w; XtPointer client_data, call_data; 

int n = (int) client_data; 

switch (n) { 

case 0: 

ShowHide2 (dummy, wExecute, False, 0) ; 
break; 
case 1 : 
break; 
case 2 : 

RunExec ( ) ; 
break; 

} 


) 


void RunExec ( ) 

char s [128] , t [128] , *u; 

XtVaGetValues (execltem[ 1 ] [0] , X tNstr ing , &u, NULL) ; 
strcpy(s,u); . 

XtVaGetValues (exec I tem[ 1] [1] , XtNs tr ing, &u, NULL) ; 
strcpy (t, u) ; 
strcat (s, t) ; 

strcat (s," &"); , . , 

XtVaGetValues (execltem[ 1] [2 j , XtNstnng , &u, NULL) ; 
strcpy (t,u) ; 
if (strlen (t ) ! =0) { 

strcat ( s , “ < • ) ; 
strcat (s , t) ; 

XtVaGetValues (execltem[ 1 ] [3] , XtNstring, NULL) ; 
strcpy (t,u) ; 
if (strlen(t) !=0) { 

strcat ( s , " >’ ) ; 
strcat ( s , t) ; 

} 

system (s) ; 

} 

/* 

xpPlots.c 


*/ 


tinclude *xplncs .h" 
♦include “xpDefs .h" 


♦include 

♦include 

♦include 

♦include 


"xpwidgets . h" 
"xpgeneral .h" 
■xpdata .h" 
■xpformats -h* 
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/* $1: Global Variables */ 

unsigned long lineColor ( 64 ] ; 
unsigned long curvColor ( 32 3 ; 


int activePlot=-l ; 

Boolean plotExists [MAXPLOTS] ; 

Boolean useCol [MAXPLOTS] [2*MAXCOLS] ; 

Boolean pi anar [ MAXPLOTS ] = { True , True , True , True , True , 


True , True , True } 


/* $2: Local variables */ 


XSegment Grid [MAXPLOTS] [ 50 ] ; 

XPoint pt [MAXPLOTS] [2*MAXCOLS] [MAXPLANES] ; 
YDninh nt2 f MAXPLOTS 1 [2*MAXCOLS] [128] ; 


char xLabels [MAXPLOTS] [24] [16] ; 
char yLabels [MAXPLOTS] [24] [16] ; 
char zLabels [MAXPLOTS] [24] [16]; 


nimension oldWidth [MAXPLOTS] ={640,640,640,640,640,640,640,640}; 
Dimension oldHeight [MAXPLOTS] ={480,480,480,480.480,480,480,480}; 


/* $3: Global functions */ 


void MakePlotO; 
void KillPlot ( ) ; 
void ActivatePlot { ) ; 
void RedrawPlot { ) ; 
void RedrawPlot 2 { ) ; 
void ResizePlot ( ) ; 
void PMActiveCallbackO ; 
void PMCreateCallbackO ; 
void PMShowHideCallback ( ) ; 
void PMPlanarCallback ( ) ; 
void PMActionCallback ( ) ; 
void IncDecPlot ( ) ; 


/* $4: Local functions */ 
void RecreatePlot ( ) ; 


void SetupGridO; 
void SetupAxisLabels ( ) ; 
void SetupPoints { ) ; 

void DrawGridO; 
void DrawAxisLabels ( ) ; 
void DrawPlotLabels { ) ; 
void DrawLegend { ) ; 
void DrawPoints ( ) ; 

void SmoothAxis Values ( ) ; 


void GESCPO; 
void poly ( ) ; 
void spline ( ) ; 

/* $5: Implementation ======*=========«============ 

void PMActiveCallback (w, client_data , call_data) 
Widget w; XtPointer client_data, call_data; 

{ int n * (int) XawToggleGetCurrent (PMitem[0] [0] } ; 

if ( (n>=mMn<=8) ) 
activePlot=n-l ? 

} 
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void PMCreateCallback(w, client data, call^data) 

Widget w; XtPointer client_data, call_data, 

* int p = (int) client_data; 

1 ^ XtVaSetValues^ PMi tem [ 0? J^XtNsens j t ive • ■ N jjull ) • 

XtVaSetValues (PMitem[2] £! '.XtLenlitive.True.itNlabel, -Hide' . NULL, ; 
MakePlot (p) ; 

ShowHide2 (plotShell [p] ,p,True, 0) ; 
else if (plotExists [p] ) { 

XtVaSetValues (PMitem[0] [pi , XtNsensitive, False, NULL) ; 

tSi .«««.> 

} 

} 

void PMShowHideCallback(w,client_data, call data) 

Widget w; XtPointer client_data, call_data, 
f int p = (int) client_data; 

1 ^XtVaSetValues 1 PMi tem [ 2 ] Ip] , XtNlabel, ‘Show' .NULL) ; 

ShowHide2 (plotShell [p] ,Pi False , 0) , 

) 

el XtvisetValues(PMitem[2] [p] . XtNlabel , "Hide" .NULL) ; 

ShowHide2 (plotShell [pi ,p,True,0) ; 

} 

} 

void PMPlanarCallback(w, client_data, call_data) 

Widget w; XtPointer client_data, call_data; 

( int which = (int) client_data; 

if ( (activePlot>=0)Sc&(activePlot<MAXPLOTS) , { 

if (planar [act ivePlot ] ) { 

plot-, will., 

} 

el se { 

XtVaSetValues^PMiteml 3? [ 0 ] .XtNlabel, "Space Plot" .NULL) ; 

CopyData (activePlot) ; 

StatData (activePlot) ; 

/* PrintPlotData (activePlot) ; */ 

SetupGrid (activePlot) ; 

SetupAxisLabels (activePlot) ; 

SetupPoints (activePlot) ; 

RedrawPlot2 (activePlot) ; 

) 

) 

void PMActionCallback(w, client_data, call_data) 

Widget w; XtPointer client_data, call_data, 

int conrcnNuiu = (int) client_data, 

switch (cominNum) { 
case 1: 

IncDecPlot (True) ; 
break; 
case 2 : 
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IncDecPlot (False) ; 
break; 
case 3 : 

if ( (activePlot> = 0 ) && (activePlot<=7 ) ) { 

if (planar [activePlot J ) 

while (currentPlane [activePlot] < (numPlanes-1) ) 
IncDecPlot (True) ; 

else 

while (currentRow[activePlot] < (numRows-1) ) 
IncDecPlot (True) ; 

) 

break; 
case 4; 

if ( (activePlot> = 0)&& (activePlo.t<=7) ) { 

if (planar [activePlot ] ) 

while (currentPlane [activePlot ] >0) 

IncDecPlot (False) ; 

else 

while (currentRow[activePlot] >0) 

IncDecPlot (False) ; 

} 

break; 
case 5: 

ShowHide2 { dummy , wFormat , True , 0 ) ; 
break; 
case 6 ; 

ShowHi de 2 (dummy , wPlotManager , True, 0) ; 
break; 

} 

) 

void MakePlot(p) 
int p; 

{ 

plotExists [p] =True; 

if (dataLoaded) { 

CopyData (p) ; 

Statlnd( ) ; 

StatData (p) ; 

/* PrintPlotData (p) ; */ 

SetupGrid(p ) ; 

SetupAxisLabels (p) ; 

SetupPoints (p) ; 

ActivatePlot (p) ; 

} 

else 

printfCNo dataset currently loaded\n"); 


void KillPlot(p) 
int p; 

{ 

ShowHide2 (dummy, wPlotO+p, False, 0) ; 
windowUp [wPlotO+p] =False; 
plotExists [p] =False; 
if (p==activePlot) { 

XawToggleUnsetCurrent (PMitem[0] [0] ) ; 
activePlot=-l ; 

) 


void ActivatePlot (p) 
int p; 

{ 

if ( (p>=0)&&(p<MAXPLOTS)&& (plotExists [p] ) ) { 

activePlot = p; 

XawToggleSet Current (PMitem[0] [0] , (caddr_t) p+1) ; 
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1 ^XtVaSetValuesTpMi tern [ 3 ] [0] .xtNlabel. 'Space Plot '.NULL) ; 
el XtVaSetValues(PMitem[3] [0] .XtNlabel, 'Plane Plot' , NULL) ; 


) 

} 


void RedrawPlot (w, event, params , numParams) 

Widget w; XEvent 'event; String params; Cardinal 


C XExposeEvent 'myevent = (XExposeEvent*) event 
Window wind=XtWindow (w) ; 
static int i=0; 
int p = atoi (params [0 ]) ; 


♦numParams ; 


if ( { *myevent) . count==0 ) { 

i++; 

XClearWindow (display , wind) ; 
DrawGrid(p); 

DrawAxisLabels (p) ; 

DrawPlotLabels (p) ; 

DrawLegend (p) ; 

DrawPoints (p) ; 

printf ( * : RedrawPlot ( %d) - #%d\n l ,p,i); 

} 

) 


void RedrawPlot2 (p) 
int p; 

Window wind=XtWindow (plot [p] ) ; 


XClearWindow (display, wind) ; 

DrawGrid(p) ; 

DrawAxisLabels (p) ; 
DrawPlotLabels (p) ; 
DrawLegend (p) ; 

DrawPoints (p) ; 


roid ResizePlot (widget , event , params , numParams ) 
wiHoet widaet; XEvent 'event; String params 


Dimension w,h; 

int p = atoi (params [0] ) ; 


Cardinal *numParams; 



f [p] .width = w; 
f [pj .height = h; 

if ( (w!=oldWidth[p] ) | | (h!=oldHeight[pl ) ) ( 

RecalcFormat(p); 

SetupGrid (p) ; 

SetupPoints (p) ; 

if ( (w<=oldWidth (p] ) Sit (h<=oldHeight [p] ) ) 
RedrawPlot2 (p) ; 


oldWidth[p] =w; 
oldHeight [p] =h; 

} 

void IncDecPlot ( inc) 
Boolean inc; 

* Boolean changed=False ; 
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if ( (activePlot>=0) && (activePlot<MAXPLOTS) ) 
if (planar [activePlot] ) ( 

if ( ( inc) && (cur rent Plane [activePlot] < (numPlanes -1 ) ) ) ( 

current Plane [activePlot ] ++ ; 
changed=True ; 

if ( ( ! inc) ScSt (currentPlane [activePlot] >0) ) ( 

currentPlane [activePlot ] 

changed=True ; 

} 

} 

else { 

if ( (inc) & & (currentRow [activePlot] <'(nuinRows-l) ) ) { 

currentRow [activePlot ] ++ ; 
changed=Truey 

if ( ( ! inc) && (currentRow [activePlot] >0) ) { 

currentRow [activePlot] 
changed=True ; 

} 

) 

if (changed) { 

CopyData (activePlot) ; 

StatData (activePlot) ; 

/* PrintPlotData (activePlot) ; */ 

SetupGrid(activePlot) ; 

SetupAxisLabels (activePlot) ; 

SetupPoints (activePlot) ; 

RedrawPlot2 (activePlot) ; 

} 

} 

void RecreatePlot (p) 
int p? 

{ 

CopyData (p) ; 

StatData (p) ; 

PrintPlotData (p) ; 

SetupGrid(p) ; 

SetupAxisLabels (p) ; 

SetupPoints (p) ; 

RedrawPlot2 (p) ; 

} 

void SetupGrid (p) 
int p; 

{ 

int dx=f [p] .xlnc, dy=f [p] .ylnc; 

int numx=f [p] . xNum, nuiny=f [p] .yNum; 

int lengthx=f [p] . xLength, lengthy=f [p] .yLength; 

int bL=f [p] .borderL , bR=f [p] .borderR; 

int bT=f [p] . borderT , bB= f [p ] .borderB; 

int i ; 

for { i = 0 ; i<=numy; i++) { 

Grid [p] [ i ] .xl = bL-f [p] .yTick ? 

Grid [p] [ i ] .yl = bT+i*dy; 

Grid [p] [ i ] . x2 = bL+lengthx; 

Grid(p] [i] .y2 = bT+i*dy; 

} 

for (i=0; i<=numx; i++) { 

Gridlp] [i^l+numy] .xl = bL+i*dx; 

Gridlp] [ i^l+numy ] .yl = bT; 

Grid [p] [ i+l+numy] .x2 = bL+i A dx; 

Grid[p] [i+l+numy] .y2 = bT+lengthy+f [p] .xTick; 

} 

} 
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void DrawGrid(p) 
int p; 

Window wind=XtWindow (plot [p] ) ? 


XSetForeground (display , gc [ 0 ] , rColor [Black] .pixel) ; 
XDrawSegments (display,wind,gc[0] ,Gnd(p] , f [p] • 


void SetupAxisLabels (p) 
int p; 

( float temp , minx , miny , maxx , maxy , minz , maxz , rangex , r angey , r angez 
int numx=f [p] . xNum , numy = f [p3 .yNum; 
int i; 

char s[ 16 ],t[ 16 ]; 

sprint f (s, *%%%sf,flpl .xFormat) ; 
sprintf (t, -%*%sf , f tpl .yFormat) ? 

if ( If [p] -opt[oAutoScale] ) { 

minx=f [p3 .xMin; 
maxx=f [pi .xMax; 
rangex=f C p ] . xMax-f tpl .xMxn; 
minz=f [p] . xMin; 
maxz=f [p] . xMax; 
rangez=f [pi . xMax- f (p ] .xMin; 
miny-f [pi .yMin; 

maxy = f [ p 1 • yMax ; 

rangey=f [p] . yMax-f [p] .yMin; 

} 

else { 

if (planar (pi) { 
minx=xMin; 
maxx=xMax; 
rangex =xRange ; 

} 

else ( 

minz=zMin; 
maxz=zMax; 
rangez=zRange ; 

) 

miny-yMin[p] ; 
maxy = yMax [p] ; 
rangey=yRange [p] ; 


if (planar [p] ) ( 

for (i = 0 ; i<=numx; i-»-+) { 

temp=minx+i* rangex /numx; 

sprintf (xLabels [p] ( i] , s , temp) ? 

for (i=0 ; i<=numy; i+^) { 

temp=miny + i * r angey /nximy ; 
sprintf (yLabels[p3 [il , t f temp) ; 

) 

) 

else { 

f or (i=0 ; i<=numx; i++) { 

temp=minz+ i * r angez / numx ; 
sprintf (xLabels [pi [ i] # s , temp) ; 

for (i=0; i<=numy; i++) { 

temp=miny + i * r angey /numy ; 
sprintf (yLabelslp] [il , t,temp) ; 

) 

} 


) 
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void DrawAxisLabels (p) 
int p; 

Window wind=XtWindow (plot [p] ) ; 
int dx=f [p] . xlnc , dy=f [p] .ylnc; 
int maxx=f [p] . xNum,maxy=f [p] .yNum; 
int i,w,x,y; 


x=f [p] • borderL; 

y=f [p] .borderT+f [p] .yLength* f [p] .xTick+12; 

XSetForeground (display , gc [ 0 ] f rColor [Black] .pixel ) ; 
for ( i=0 ; i<=maxx; i + +) { 

w=XTextWidth ( f s6x!0 , xLabels [p] [i] , strlen {xLabels [p] [i]))/2; 
XDrawString (display , wind, gc [0] ,x+i*dx-w,y, 
xLabels [p] [i] , strlen {xLabels [p] [i] ) ) ; 

x=f [p] .borderL-f [p] .yTick-6; 
y=f [p] . borderT+f [p] .yLength; 

for(i=0;i<=maxy;i*+) { f , , . 

w=XTextWidth ( f s6xl0 , yLabels [p] [i] , strlen {yLabe Is [p3 [iJ ) ) ; 
XDrawString (display, wind, gc [0] , x-w,y-i*dy, 
yLabels [p] [i] , strlen (yLabels [p] [i] ) ) ; 

) 


void DrawPlotLabels (p) 
int p; 

Window wind=XtWindow(plot [p] ) ; 
char t [ 64 ] ; 
int w , n , h , v ; 


} 


XSetForeground {display, gc [0] , rColor [Black] .pixel) ; 
strcpy{t, f [p] .title) ; 

w=XTextWidth ( f sHelvlSb, t , strlen (t) ) ; 
h=f [p] « borderL^ ( f [p] . xLength-w) / 2; 
v=f [pj .borderT-8; 

XDrawString (display , wind, gc [2 ] , h, v, t, strlen (t) ) ; 


if (planar[p]) 

sprintf (t, ■Distance (inch) [Time = %6.0f : %3d/%3d] 
z [current Plane [p] ] , currentPlane [p] +1 , numPlanes) ; 


sprintf ( t , "Time (msec) [Dist - %6.0f : %3d/%3d]", 
x [currentRow[p] ] , currentRow[p] +l,numRows) ; 
w = XTextWidth { f s6xl0 , t , strlen ( t ) ) ; 
h=f [p] .borderL* ( f [p] . xLength-w) / 2 ; 
v=f [pj ,borderT*f [p] .yLength+30; 

XDrawString (display, wind, gc [0] ,h,v, t,strlen(t) ) ; 


void DrawLegend { p ) 
int p; 

Window wind=XtWindow (plot [p] ) ; 
int j,k=-l; 

XSegment Legend [1]; 

int len=15 , xDist=100 , yDist=14 , x, y; 

int test; 

test = ( planar [p] ? MAXCOLS ; 2*MAXCOLS ); 

XSetLineAttributes (display, gc [0 ] , 8 , LineSolid, CapButt , JoinRound) 
for ( j=0; j<test; j*+) 
if (useCol [p] [ j ] ) ( 

k++; 

x=f [p] .borderL* (k%4) *xDist; 

y=f [p] .borderT+f [p] .yLength* (k/4) *yDist+45; 

Legend [0] . xl=x; 
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} 


Legend [0] .yl=y-4; 

Legend [0] . x2=x+len? 

XSetForeground (display , gc [ 0 ] , rColor [Black] .pixel) , 

XDrawStr ing ( di splay , wind , gc [ 0 ] , x+ len+5 , y , 

columnNames [ j ] , strlen (columnNames n 3 ) ) ; .-mu. 

XSetForeground (display , gc [0] , ImeColor [ f t P 1 .color u 3 3 ' 

XDrawSegments (display , wind, gc [0] , Legend, 1) ; 

xsetForearound (display , gc [0] , rColor [Black] .pixel ) , , 

XSetLineAttributes (display , gc [ 0] , 0 , LineSolid , CapButt , JomRound) 


void SetupPoints (p) 
int p; 

{ int bL=f [p] .borderL, bR=f [p] .borderR; 
int bT=f [p] .borderT, bB=f [p] .borderB; 

float minx=xMin,miny=yMin[p] ,minz=zMm; 
float rangex=xRange , rangey=yRange [p] ,rangez-zRange, 
int lengthx=f [p] .xLength, lengthy=f [p] .yLength; 
int i . i ; 


if ( ! f [p] .opt [oAutoScale] ) { 

minx=f [p] • xMin; 
rangex=f [p] . xMax-f [p] .xMin; 
minz=f [p] .xMin; 
rangez=f [pi .xMaoc-f [p] .xMin; 
miny=f [p] .yMin; 
rangey=f [p] . yMax-f [p] .yMin; 

) 

else { 

if (planar [p] ) { 

minx=xMin; 
rangex=xRange ; 

} 

else { 

minz=zMin; 

rangez=zRange; 

} 

miny=yMin [p] ; 
rangey=yRange [p] ; 


( 


if (planar [p] ) ( 

for ( j=0 ; j<numCols ; j++) 
if (useCol [p] ( j ] ) 

f or ( i=0 ; i<numRows ; i++) { 

pt[p] [j] [i] . x=bL+ ( (x [i] -minx) /rangex) *lengthx; 
pt[p] [ j] C i ] .y=bT+ lengthy- { (plotData[p] [^1 [i] -mmy) / 
rangey) * lengthy; 


} 

} 

else 


) 


{ 


{ 


for ( j=0; j< (2*MAXCOLS) ; j++) 
if (useCol [p] [j] ) 

for ( i = 0 ; i<numPlanes ; i++) { 

pt[p] [j] [i] . x=bL+ ( ( z [ i 1 -minz) /rangez) *lengtnx; 
pt[p] [ j ] { i 1 .y=bT+lengthy- ( (plotData[p] []] [i] mmy) / 
rangey) * lengthy; 

) 

) 


void DrawPoints (p) 
int p; 

{ 
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Window wind=XtWindow(plot [p] ) ; 
static char dashPattern [ ] = {16,8}; 
int j , test; 

XRec tangle r; 

int bL=f [p] . borderL , bR= f E p 3 .borderR; 
int bT=f [p] . bor derT , bB= f [ p ] .borderB; 
int lengthx=f Epl .xLength, lengthy=f [p] .yLength; 

XSetLineAttributes (display, gc [0] ,4 , LineOnOf fDash, CapRound, JoinRound) ; 
XSetDashes (display, gc [0] , 0 , dashPattern, 2) ; 

r . x=bL-2 ? r . y=bT-2 ; r . width=lengthx+4 ; r . height=lengthy+4 ; 
XSetClipRectangles (display, gc [0] , 0, 0, r, 1, Unsorted) ; 

test= (planar [p] ? MAXCOLS : 2*MAXCOLS); 
for { j=0; jctest; j++) 
if (useCol [p] [ j ] ) { 

if (f Ep) .numDashes (j ] >0} { , , _ 

XSetLineAttributes (display, gc [ 0] , 4 , LineOnOf fDash, CapRound, JoinRound) 
XSetDashes ( display, gc [ 0 ] , 0 , f Ep] . dashes [ j ] , f E p 3 • numDashes E 3 J ) * 

} 

else ■ , 

XSetLineAttributes (display, gc (0] , 4 , LineSolid, CapRound, JoinRound) 
XSetForeground (display , gc [ 0 ] , lineColor [ f Epl . color [ j ] 3 ) ; 
if (planar[pj) _ _ ^ . . . 

XDrawLines (display, wind, gc [0 3 ,pt[p] [ j ] , numRows , CoordModeOngm) ; 

else 

XDrawLines (display , wind, gc [0] , pt Ep] E j 3 • numPlanes , CoordModeOrigin) 

XSetForeground (display, gc [0] , rColor EBlack] .pixel) ; 

XSetLineAttributes (display, gc [0J , 0 , LineSolid, CapButt , JoinRound) ? 
XSetClipMask (display , gc (0 3 ,None) ; 


void SmoothAx is Values (vMin, vMax) 
float vMin, vMax; 

{ 

float vRange ; 

int a, al , a2 , b, bl , b2 ; 

vRange = vMax-vMin; 

a = loglO (vRange) ; 
al= a-1 ; 
a2= a-2 ; 


b » vRange/a; 
printf ( ’%d' ,b) ; 


The following routines apply only to curve-fitting (not implemented) 


void GESCP (n, a,b) 

int n; float a [ 10] [ 10 ] , b [ 10) ; 

{ 

int i , j , k, p, nrow[ 10] , ncopy? 
float max,m, sum, s [10) ; 

for ( i=0 ; i<n; i++ ) { 
max = -1.0E+10; 
for (j=0; j<n; 

if ( fabs (a [ i ] [ j J ) >max) max = fabs (a [i] E j 3 ) ; 
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printf (“No unique solution exists\n"); 


s [ i ] -max ; 
nrow( i ] =i ; 

for ( i=0 ? i< (n-1 ) ; i++ ) { 
max=-l . 0E+10 ; 

f0 if <3 r<fabsulirow[jl] [i])/s[nrow[j]]) > max) { 
max = f abs ( a [nrow [ j 1 ] [ i ] ) / s [nrow [ ] ] J ; 

P = j ? 

if (a[nrow[p]l [i)==0.0) 
if (nrow[i] !=nrow[p] ) ( 

ncopy=nrow[i] ; 
nrow [ i ] =nr ow [ p 3 ; 

nrow[p] =ncopy; 

) . % , 
for ( j= (i + 1) ; 3 <n; 3+ + ) l 

m = a[nrow[ j ] ] [i] /a [nrow[i] ] [i] * 

for (k=0;k< (n+1) ;k++) ....... 

a [nrow[ j ] ) [ k] =a [nrow[ j ] 3 [k] -m*a[nrow[ i] 3 [k] ; 

} 

if (a[nrow[n-U ] [n-1] ==0 . 0) printf (‘No unique solution exists\n-) 
b[n-13 =a [nrow [n-1] 3 [n] / a [nrow [n 13 3 [n-1] , 
for ( i= (n-2 ) ; i>-0; i ) ( 
sum=0 ; 

for ( j=i+l; j<n; j++) 

sum=sum+a [nrow[ i) 3 [ 3 3 bu]; 
b [ i ] = (a [nrow[i] 3 [n] -sum) /a [nrow [ i] J [ i] ; 

) 


void poly (n,m,x,y) 

int n , m; float x[3 <!/[) » 


float* a[10 ] [10] ,b[10] ,p[256] ,e[256] ,sum; 
for ( i=0 ; i<n; i++ ) { 
for (j=0; j<n; j++) { 
sum=0 ; 

f 0 sum=sumM m ( i+ j ) ==0 ? 1.0 : (float) pow( (double) x [k] , (double) l+j) 
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if ( { i + j ) ==0 ) sum=m; 
a [ i ] [ j]=sum; 


for ( i=0 ; i<n; i+ + ) ( 
sum=0 ; 

for (k=0 ; k<m; k++) 

sum=sum+y [k] * ( i==0 ? 1.0 
a [ i 3 [n] =sum; 


( float) 


} 

GESCP (n, a,b) ; 

for (i=0; i<n; i++) , ..... 

printf ( "x [%d] = %12 . 6f \n* , i ,b[i] ) ; 
for ( i=0 ; i<m; i++) { 
sum=0 ; 

for (k=0 ; k<n; k++) 

sun\=sum+b [k] * (k==0 ? 1.0 : (float) 


pow{ (double) 


pow( (double) 


x [k] , (double) 


x [ i] , (double) 


p [ i ] =sum; 
e[i] =y [ i 1 -pti) 
printf ( "%12 . 6f 


% 12 . 6 f% 12 . 6 f% 12 . 6 f\n" 


x[ i] ,y [i] » P C il »e[i) ) ; 


) 


sum=0 ; 

for ( i = 0 ; i<m;i++) 
sum+=( float) pow 
printf ( " \n\nTotal 


((double) e [ i) , (double) 
error = %12 . 6f \n* , sum) ; 


2 ) ; 


i) 


k) 


) ; 
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void spline (n, x,y) 

int n; float x[] ,y[] ; 

{ 

int i,j; 

float a [256] ,b[256] ,c[256] ,d[256] , 

h [256] , alpha [256] ,1[256] ,u[256J , z[256] ; 


} 


for { i=0 ; i<=n; i++ ) 
a[i]=y[i] ; 
for ( i = 0 ; i<n; i++) 
h[i] =x[i+l] -x[i] ; 
for (i=l; i<n?i++) 

~ i r ^ l ' ' ' ’ - 


(h[i-l]*h[i] ) ; 

1 [ 0 ] = 1 ; 
u [ 0 ] = 0 ; 
z [0 ] =0 ; 

for (i=l;i<n;i++) { 

l[il*2*(xli+l]-x[i-ll ) -h[i-l] *u [i-l] ; 
u [ i ] =h [ i 3 / 1 [ i ] ; 

z [ i j = ( a lpha [i]-h[i-l]*z[i-l] )/l[i] ; 


} 

1 [n] =1.0 
z [n] =0 . 0 
c [n] =0 . 0 


for ( j= (n-1) ; j>=0; j — ) { 

c[j]=z[j]-u[j]*c[j+l] ; 

b[j] = (a[j + l]-a[j])/h[j]-htj]Mc[j + l]+2*c[j]>/3; 
d[j]»(c[j+l]-c[j] )/ ( 3*h [ j ] ) ; 


) 


/* 


xpWidgets . c 


*/ 

♦include "xplncs.h" 

♦include •xpDefs.h" 

j * following files have to be included for callback routines */ 

♦include "xpmenus.h" 

♦include •xpcommands . h* 

♦include "xpinput.h" 

♦include *xpformats . h" 

♦include ■xpplots.h" 

♦include ■xpgeometry . h" 

/* $1: Global variables */ 

Widget topLevel; 

Widget rootPane; 

Widget menuPane, menuText ; 

Widget statusLabel; 

Widget coironPane , verbViewport , verbList , coiranViewport , commList ; 

Widget cliPane, cliLabel , cliText ; 

Widget inMenuBut ton [NUMMenus] , mS imp leMenu [NUMMenus] ; 

Widget m£n try [NUMMenus ] [MAXENTRIES] ; 

Widget commandShell, commandBox, c [NUMCommands] ; 

Widget PMShell, PMForm, PMBox [4] , PMTop [4] , PMitem[4] [8] ; 

Widget infoShell , inf oCore; 
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Widget logShell , logText ; 

Widget f ileShell, filePane, f ileMenuPane, fileText, fileEntry[4] ; 

Widget IEShell , IEForm, IEBox[8] , IETop[8] , IEItem[8] [2] , 
IEViewport, IEList; 

Widget DAShell , DAForm, DABox [ 5 ] ,DATop[5] ,DAItem[5] 116] , 

DALabel [ 2 ] ; 

Widget CFShell , CFForm, CFBox [4 ] ,CFTop[4] ,CFItem[4] [16] ; 

Widget MGShell , MGForm, MGBox [ 6 3 ,MGTop[6] ,MGItem[6] [12} ; 

Widget MGPictShell , MGPict ; 

Widget FAShell , FAForm, FABox [ 5 ] ,FATop[5] ,FAItem[5] [16] ; 

Widget execShell , execForm, execBox [4 ] ,execTop[4] ,execltem[4] [8] ; 

Widget plotShell [MAXPLOTS] , plot [MAXPLOTS] ,* 

Widget dataShell, dataText; 

Widget forma tShell , forma tText; 

Widget FEShell , FEForm, FEBox [ 16 ] ,FETop[16] , FEItem[16] [16] ; 

Widget entry , dummy ; 

Widget sfGetFileShell , sfDialog, sfCancel , sfLoad, sfSave; 

/* $2: Local variables */ 


static char inMenuButtonNames [ ] [ 12 ] = { 

■mBarFile* , "mBarlnpu" , "mBarExec* , "mBarData* , 

■mBarForm" , "mBarPlot" , "mBarWind - , } ; 
static char mSimpleMenuNames [ ] [12] = { 

"menuFile" , "menulnpu" , " menuExec 1 , "menuData 1 ’ , 

■menuForm" , "menuPlof , "menuWind" , } ; 

static char mEntryNames [NUMMenus] [12 ] [24] = { 

{ •mFile_00* , "mFile_01- , *mFile_02 - , • / * , "mFile.OO - , 

•mFile_04" , "mFile_05' , ■ / * , •mFile_06' , *mFile_07" , 

■/ - , m wFile_Q8 '} , 

{ "mlnpu_00" , *mlnpu_01 * , "mlnpu_02* , ■mlnpu_03 " , ’mlnpu_05" , 
-/ " , ■mlnpu_06" , ■mlnpu_07" , " " , ■ 

■ ■ ■ * } 

{ "mExec_00" , *mExec_01") , 

{ "mData_00 * , *mData_01 • , *mData_02 * , * / * , "mData_03 ’ , 

*mData_04 * , "mData_05 ’ , "mData_06 ■ , "mData_07 * , "niData_08 * , 

w 9 f 9 " } , 

{ "mForm__00 ■ ] "mForm_01“ f “mForm_02" , "mForm_03 " , "mForm_04 - , 
- / - , •mForm_05 " , - mForm_06 " , “mForm_07" , "mForm^OS" , 

“mForm_09" , "mFonn^lO " } , 

{ "mPlot_00* , *mPlot_01* f -inPlot_02 1 , - / " > ’ , mPlot_03 ■ , 

■mPlot_04 " , "mPlot_05 * , *mPlot_06 ■ , "mPlot_07 ■ , " ■ 

■ ■ ■ 9 } t 

{ "mWind_00" -inWind_01* , *mWind_02* , "mWind_03" , *mWind_04- , 
"mWind_05" , *mWind_06 * , "mWind_07 B , •mWind_08' , “mWind_09 ■ , 
*mWind_10“ , "mWind_ll" }, 


static int numMenuEn tries [NUMMenus ] = { 

NUMFileEntries , NUMInpuEntries , NUMExecEntries , NUMDataEntr les , 
NUMFormEntries , NUMPlotEntries f NUMWindEntries } ; 


static String coirunElements [ ] = { 

* 1 , 1 ", * " i * ’ r * > * * • 1 * 


, NULL) ; 


static String verbs [] ={ . , „ , ..... 

■active" , 'change*, 'color*, "data* , ’decrement* , "edit* , "execute , fit , 

* format* , "hide* , "increment* , "input" , "kill* , "list" , "load" , "make" , 
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■plane* , 'plot' , -print’ , -quit - , -save" , 
■set - , "show - , -space" , -view* , NULL, } ; 

static String IEentries [] = ( 

■ - , NULL } ; 

*/ 

♦define NUMFEOptions 2 
♦define NUMFECommandString 4 
♦define NUMFEAxes String 4 

static char s[64],t[64]; 


/* $3: Global functions */ 
void CreateWidgets ( ) ; 

/* $4: Local functions - (none) */ 
/* $5: Implementation ============ 



*/ 


void CreateWidgets ( ) 

static char mMenuButtonNames ( } [ 10] = { 

■mBarFile" , -mBarlnpu - , -mBarExec" , -mBarData - # 

-mBarForm- , -mBarPlot" , "mBarWind" , } ; 
static char mSimpleMenuNames [ ] [10] = { 

■menuFile- , -menulnpu" , "menuExec" , ’menuData’ , 

"menuForm* , "menuPlot" , "menuWind" , } ; 
static char mEntryNames [NUMMenus ] [ 12 ] [ 12 ] = { 

{ -mFile_00- , ’mFile_01" , "mFile_02 ■ , * / " , "mFile_03" , 

"mFile_04 " , ’mFile_05" , ■ / " * "mFile_06" , -mFile_07" , 

•/ , "mFile_08 ’ } , 

{ "mlnpu_00’ , -mlnpu_01- , "mlnpu_02" , *mlnpu_03* , ‘mlnpu_05’ , 

■ / ■ , -mlnpu_06" , ’mlnpu_07" , " " - * 

■ ■ * * } r 

{ "roExec_00 * , *mExec_01 " } , 

{ "mData_00* , *mData_01" , *mData_02* , "/ •,*mData_03 , 

*mData_04 ' , *mData_05 * , "mData_06* , ‘mData_07 ' , ■mData_08 ’ , 

■ *” * * " } f 

{ - mForm_00 " , -mForm_01 - , ’mForm_02 - , "mForm_03 - , "mForm_04 " , 

- / ■ # ■mForm_05 - , - mForm_06 " , "mForm_07 - , - mForm_08" , 

•mForm_09 - , ’mForm_.10’ } , „ 

{ -mPlot_00" , "mPlot__01 “ , -mPlot_02 * , ’ / ", -mPlot_03 - , 

"mPlot_04 ■ , ’mPlot_05 " , "mPlot_06 - , "mPlot_07 ■ , • 

■ ■ ■ ■ ) , 

{ ’mWind_00’ \ -mWind_01 - , -mWind_02-, "mWind_03-, ’mWind_J)4 - , 

"mWind_05 ’ , -mWind_06" , - mWind_07- , - mWind_08- , ’mWind_09 - , 

"mWind_10 " , ’mWind_ll " }, 

static int numMenuEntries [NUMMenus] = { 

NUMFileEntr ies , NUMInpuEntries , NUMExecEn tries , NUMDataEn tries , 
NUMFormEntries , NUMPlotEntries , NUMWindEntries } ; 
static String commElements E ] = { i 

- •/* V *" *'" * ' " 

static String verbs [] = { „ . „ „ 

■active" , "change" , - color’ , "data" , -decrement’ , - edit’ , ’execute , fit , 
■format’ , -hide" , "increment’ , ’input’ , -kill" , "list" , "load" , "make - , 
■plane" , -plot’ , -print’ , -quit" , -save’ , "set- , 

■show" , -space" , -view’ ( NULL, } ; 
static String IEentries [3 = ( 

■ - , NULL } ; 
int i, j, total; 

/* Generate the pane that will serve as the root pane for the window */ 
rootPane =XtVaCreateManagedWidget ( “rootPane - ,panedWidgetClass, topLevel, 
XtNorientation, XtorientVertical , 

NULL) ; 
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' \lnultnt ^XtVaCreateManagedwTdgeC^menuPane^ . ^edWidgetClass , rootPane . 
XtNorientation, XtorientHorizontal , 

XtNshowGr ip, False, 

NULL) ; 

/* .Generate the menu buttons, placing them one after another */ 
for ( i = 0 ; i<NUMMenus ; i + + ) 

mMpnuButton [ i ] = XtVaCreateManagedWidge t ( 

XnuButtonNames [ i 1 , menuButtonWidgetClass , menuPane , 

XtNborderWidth, 0 , 

XtNshowGr ip, False, r • i 

XtNmenuName , mSimpleMenuNames [ i ] / 

NULL) ; 

for ( j=0; j<NUMMenus; j++) ( . ... 

mS Si^ SJe^aLs [ j 1? Bi£l2enSwidgetClass , mMenuButton [ j ] , 

NULL) ; . ... , 

for ( i=0 ; i<numMenuEntries [ 3 J ; i ++ J l 

if (mEntryNames [ j ] [ i] [0] == ' / ' ) . 

mEntrv HI [ i] = XtVaCreateManagedWidget ( .... 

mEntryNames [ j ] [i] , smeLineObjectClass ,mSimpleMenu [ 3 , 

NULL) ; 

el mPntrvfi] [i] = XtVaCreateManagedWidget ( 

^mEntryNames [ 3 ] [i] , smeBSBObjectClass.mSimpleMenuH ) . 
XtAddCallback (mEntry [ j ] [i] ,XtNcallback,MenuSelect, 

(XtPointer) ( j *HIW0RD+i) ) ; 

) 

} 

menuText = XtVaCreateManagedWidget ( 

■menuText • , asciiTextWidgetClass , menuPane, 

XtNborderWidth, 0 , 

XtNshowGr ip , False , 

XtNeditType , XawtextEdit , 

NULL) ? 

statusLabel = XtVaCreateManagedWidget ( 

■ statusLabel ■ , labelWidgetClass , rootPane , 

XtNwidth, 600, 

XtNshowGr ip , False , 

XtNskipAdjust , True, 

NULL ) ; 

commPane = XtVaCreateManagedWidget ( 

■commPane" , panedWidgetClass , rootPane, 

XtNorientation, XtorientHorizontal , 

XtNshowGr ip , False , 

verbViewport * XtVaCreateManagedWidget ( 

• verbViewport " , viewportwidgetclass , commPane , 

XtNmin, 100, 

XtNal 1 owVert , True , 

XtNf orceBars , True , 

XtNshowGrip , False , 

verbList = XtVaCreateManagedWidget} 

■verbList' , listWidgetClass, verbViewport, 

XtNverticalList,True, 

XtNcolumnSpacing , 0 , 

XtNrowSpacing, 2, 

XtNdef aultColumns , 1 , 

XtNf orceColumns , True , 

XtNlist, verbs, 

XtAddCallback (verbList, XtNcallback, VerbCallback, 0) ; 
commViewport = XtVaCreateManagedWidget ( 
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■commViewport" , viewportWidgetClass , commPane , 

XtNa 1 lowVert , True , 

XtNf orceBars , True , 

XtNshowGrip , False , 

NULL) ; 

commList = XtVaCreateManagedWidget ( 

"commList" , listWidgetClass , commViewport , 

XtNverticalList , True , 

XtNcolumnSpacing , 0 , 

XtNdef aultColumns , 1, 

XtNf orceColumns , True, 

XtNlist, commElements , 

XtNnumber Strings , 15, 

NULL) ; 

XtAddCallback(commList , XtNcallback, ListCallback, 0 ) ; 

cliPane = XtVaCreateManagedWidget < 

"cliPane" , panedWidgetClass , rootPane , 

XtNskipAdjust , True, 

XtNorientation, XtorientHorizontal, 

XtNshowGrip , False , 

NULL ) ; 

cliLabel = XtVaCreateManagedWidget ( 

"cliLabel " , labelWidgetClass, cliPane, 

XtNlabel , " ", 

XtNshowGrip , False , 

XtNjustify,XtJustifyLef t , 

XtNmin, 100, 

XtNresize, False, 

NULL) ; 

cliText = XtVaCreateManagedWidget ( 

"cliText • , asciiTextWidgetClass, cliPane, 

XtNeditType , XawtextEdit, 

NULL) ? 

/* _Create the widgets used to make the command palette work */ 

commandShell = XtVaCreatePopupShell ( 

• commandShell " , transientShellWidgetClass , topLevel , 

XtNtitle, "Window Manager", 

NULL) ; 

commandBox ^XtVaCreateManagedWidget ( "commandBox" , boxWidgetClass , 
commandShell , NULL) ; 

for ( i=0 ; i<12 ; i++) { 
sprintf (s, "c%02d" , i) ; 
c I i] = XtVaCreateManagedWidget { 
s , commandWidgetClass , commandBox, 

XtNborderWidth, 0 , 

XtNwidth, 200, 

XtNresize, False, 

XtNj us t i f y , Xt Jus t i f yLe f t , 

NULL); , % >t 

XtAddCallback ( c [ i ] , XtNcallback, CommandSelect , (XtPointer) i) ; 

} ; 

I* ^create the widgets used to make the command line log work */ 

/* 

logShell = XtVaCreatePopupShell ( 

" logShell " , transientShellWidgetClass , topLevel , 

XtNti tie, "Command Line Log", 

NULL) ; 

logText = XtVaCreateManagedWidget ( 

"logText" , asciiTextWidgetClass, logShell, 

XtNwidth, 480, 

XtNheight , 240 , 

XtNeditType , XawtextAppend , 

XtNscrollVertical , XawtextScrol lAlways , 

XtNshowGrip , False , 

NULL) ; 
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/* _Create the widgets used to view plots (up to MAXPLOTS-8) */ 
for (i = 0; i<MAXPLOTS; i++) { 
char name [ 64 ], label [ 64 ] ; 
sprintf (name , "plotShell__%d* , i) ; 
sprint f ( label , "xplot : Plot #%d’,i); 
plotShell ( i ] = XtVaCreatePopupShell ( 

name, transientShellWidgetClass, topLevel, 

XtNminWidth, 320, 

XtNminHeight , 240 , 

XtNsaveUnder f False , 

XtNtitle , label, 

NULL); . 

sprintf (name , “plot %d* , i ) ; 

plot [ i] = XtVaCreateManagedWidget { 

name , widgetClass , plotShell { i ] , 

XtNwidth, 640, 

XtNheight, 480, 

NULL) ; 

) 

/* .create the widgets used to view the current dataset */ 
dataShell = XtVaCreatePopupShell ( 

•dataShell" , transientShellWidgetClass, topLevel, 

XtNtitle, ’Dataset’ , 

dataText = XtVaCreateManagedWidget ( 

"dataText" , asciiTextWidgetClass , dataShell, 

XtNwidth, 480, 

XtNheight, 240, 

XtNedi tType , XawtextRead , 

XtNtype.XawAsciiFile, 

XtNscrollVertical.XawtextScrollAlways, 

XtNstring, ’XPlot. dat’, 

NULL) ; 

/* _ C reate the widgets used for a standard get|put file dialog 
/* 

sfGetFileShell = XtVaCreatePopupShell ( 

■ sfGetFileShell ■ , transientShellWidgetClass , topLevel , 
XtNtitle, "Filename Dialog", 

NULL) ; ^ ^ 

sf Dialog = XtVaCreateManagedWidget ( 

• sf Dialog* ,dialogWidgetClass, sfGetFileShell, 

XtNlabel, "Load Dataset (Named) 

XtNvalue , " ’ , 

NULL); a .. V 

sf Cancel = XtVaCreateManagedWidget ( 

■sfCancel" , commandWidgetClass , sf Dialog, 

XtNlabel, "Cancel" , 

NULL); ^ ^ 

sfLoad * XtVaCreateManagedWidget ( 

"sfLoad" , commandWidgetClass, sfDialog , 

XtNlabel, "Load" , 

NULL); 

sfSave * XtVaCreateManagedWidget ( 

" sfSave " , commandWidgetClass , s fDialog , 

XtNlabel, "Save" , 

XtAddCal lback ( sfCancel, XtNcallback, SFCallback, (XtPointer) 0) 
XtAddCallback (sfLoad , XtNcallback, SFCallback, (XtPointer) 1), 
X tAddCal lback { sfSave , XtNcallback, SFCallback, (XtPointer) 2) 

*/ 

/* _Create the widgets used in the dialog to edit inputs */ 
IEShell = XtVaCreatePopupShell ( 

■IEShell" , transientShellWidgetClass, topLevel, 

XtNtitle, "Input File Editor", 

NULL) ; 
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IEForm = XtVaCreateManagedWidget ( 

" IEForm" , f ormWidgetClass , IEShell, 

NULL.) ; 

for { i=0 ; i<8 ; i++) { 

spr intf { s , ■ IEBox_%d " , i ) ; 

IEBox [ i ] = XtVaCreateManagedWidget ( 
s , boxWidgetClass , IEForm, 

XtNhSpace , 1 , 

XtNvSpace , 1, 

XtNorientation,XtorientHorizontal , 

NULL) ; 

sprint f (s , " IETop_%d* , i ) ; 

IETop[i] = XtVaCreateManagedWidget { 
s , labelWidgetClass , IEBox ( i ] , 

XtNresize, False, 

NULL) ; 

) 

for { i=0 ; i<6 ; i++) { 

sprintf (s, ■ IEItem_%d_0" , i) ; 
if ( (i<3) | ] { i==5) ) 

IEItemfi] [0] = XtVaCreateManagedWidget ( 

s , asciiTextWidgetClass , IEBox [i] , 

XtNdisplayCaret , False , 

NULL) ; 

else . 

IEItem[i] [0] = XtVaCreateManagedWidget { 
s , asciiTextWidgetClass , IEBox[i] , 

X tNedi tType , Xawtex tEdi t , 

NULL) ; 

} 

for ( i=0 ; i<6 ; i++) { 

sprintf (s, • IEItem_6_%d“ ,i) ; 

IEItem[6] [i] = XtVaCreateManagedWidget ( 
s , commandWidgetClass , IEBox [ 6 ] , 

NULL) ; 

XtAddCallback ( IEItem( 6 } [i] , XtNcallback, IECallback, (XtPointer) l) ; 

) 

IEViewport = XtVaCreateManagedWidget ( 

"IEViewport* , viewportWidgetClass , IEBox [7] , 

XtNal lowVer t , True , 

XtNf orceBars , True , 

XtNtop , XtChainTop , XtNf romHoriz , IEBox [0] , 

NULL) ; 

IEList = XtVaCreateMemagedWidget ( 

■ IEList" , listWidgetClass , IEViewport, 

XtNverticalList , True, 

XtNcolumnSpacing, 0, 

XtNrowSpacing , 2 , 

XtNdefaultColumns, 1, 

XtNf orceColumns , True , 

XtNlist , IEentries , 

NULL) ; 

XtAddCallback (IEList, XtNcallback, IEListCallback, 0) ; 


XtVaSetValues ( IEBox [ 0 ] , XtNlef t , XtChainLef t , XtNtop , XtChainTop , NULL) ; 
for (i=l; i<7 ; i++) 

XtVaSetValues (IEBox [i] ,XtNfromVert, IEBox [l-l] , 

XtNlef t, XtChainLef t, NULL) ; 

XtVaSetValues { IEBox [ 7 ] , XtNf romHoriz , IEBox [ 0 ] , XtNtop , XtChainTop , 
XtNorientation,XtorientVertical,NULL) ; 

/* ^Create the widgets used in the dialog for the colors and hashes */ 

DAShell = XtVaCreatePopupShell ( 

■DAShell" , transientShellWidgetClass , topLevel, 
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XtNtitle, “Data Attributes Editor" , 
NULL); 

DAForm = XtVaCreateManagedWidget! 

* DAFonn’ , f ormWidgetClass , DAShel 1 , 
NULL) ; 


/ * 

DALabel [ 0 ] = XtVaCreateManagedWidget ( 

" DALabel_0 " , labelWidgetClass , DAForm, 
XtNf romHoriz , DABox [ 0 ] , 

NULL); 

DALabel [1] = XtVaCreateManagedWidget < 

" DALabel_l " , labelWidgetClass , DAForm, 
XtNf romHoriz , DALabel [ 0 ] , 

NULL) ; 




for ( i=0 ; i<5 ; i++) ( 

sprint f (s, "DABox_%d" , i) ; 

DABox [i] = XtVaCreateManagedWidget ( 

s , boxWidgetClass , DAForm, 
XtNhSpace, 1, 

XtNvSpace , 1 , . _ 

XtNorientation, XtorientVertical , 

XtNtop , XtChainTop , 


NULL) ; 

sprint f ( s , "DATop_%d" , i) ; 

DATop [ i ] = XtVaCreateManagedWidget ( 

s, labelWidgetClass, DABox [i] , 
XtNresize, False, 


NULL) ; 


for ( i=l ; i<5 ; i++) 

XtVaSetValues (DABox 


i] , XtNf romHoriz, DABox [i-1] , NULL) ; 


for (i* 0 ;i< 16 ;i++) { 

sprint f (s, "DAI tem_0_%0 2d" , i) ; . 

DAItem[0][i] = XtVaCreateManagedWidget! 
s, labelWidgetClass, DABox [0] , 
XtNresize, False, 

XtNlabel, " ", 

NULL) ; 


for (j=l; j<5; j++) 

for ( i=0 ; i<16 ; i++) { , 

sprintf (s , "DAItem_%d_%02d " . 
DAItem[ j ] [ i ] = XtVaCreateManagedWidget ( 

s , asciiTextWidgetClass , DABox [j ] , 

XtNeditType , XawtextEdit , 

NULL) ; 


/* _Create the widgets used in the dialog to edit formats */ 


FEShell = XtVaCreatePopupShell ( 

■FEShell " , transientShellWidgetClass , topLevel , 
XtNtitle, "Format Editor", 

NULL); 

FEForm = XtVaCreateManagedWidget ( 

"FEForm" , f ormWidgetClass , FEShell , 

NULL) ; 


for (i=0;i<10;i++) { 

sprintf (s , "FEBox_%d" , i ) ; 

FEBox [ i ] = XtVaCreateManagedWidget ( 

s , boxWidgetClass , FEForm, 

XtNhSpace , 1 , 

XtNvSpace, 1, 

NULL) ; 

sprintf (s, "FETop_%d" , i) ; 

FETop [ i ) = XtVaCreateManagedWidget ( 
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s , label WidgetC lass , FEBox [ i ] , 

XtNresize, False, 

NULL) ; 

) 

XtVaSetValues (FEBox [91, XtNlef t , XtChainLef t , XtNtop , XtChainTop , 
XtNorientation, XtorientHorizontal , NULL) ; 

XtVaSetValues (FEBox [ 0 ] , XtNlef t , XtChainLef t , XtNf romVert , FEBox [ 9 ] , 
XtNorientation, XtorientHorizontal , NULL) ; 

XtVaSetValues (FEBox [1] , XtNlef t , XtChainLef t , NULL) ; 

XtVaSetValues (FEBox [ 5 ] , XtNlef t , XtChainLef t , NULL ) ; 
for (i=l; i<5; 

XtVaSetValues (FEBox [ i 3 , XtNf romVert , FEBox [ 0 ] , 

XtNf romHoriz , FEBox [ i-1 ] # NULL) ; 
for ( i=5 ; i<9 ; i++) 

XtVaSetValues (FEBox[i) , XtNf romVert , FEBox [ 1 ] , 

XtNf romHoriz, FEBox [i-1] , NULL) ; 

XtVaSetValues (FEBox[lj , XtNf romHoriz, NULL, NULL) ; 

XtVaSetValues (FEBox [5] , XtNf romHoriz , NULL, NULL) ; 


FEItem[ 9] [ 0 ] = XtVaCreateManagedWidget ( 

*FEItem_9_0 * , asciiTextWidgetClass , FEBox [ 9 ] , 

XtNedi tType , XawtextEdi t , 

NULL) ; 

FEItemfO] [0] = XtVaCreateManagedWidget ( 

■ FEItem_0_0 ■ , asciiTextWidgetClass , FEBox [ 0 ] , 

XtNedi tType , XawtextEdi t , 

NULL) ; 

for (i=0;i<8;i++) { 

sprint f (s, •FEItem_l_%d" , i) ; 

FEItemtl] [i] = XtVaCreateManagedWidget ( 
s, toggleWidgetClass, FEBox [1] , 

XtNres i ze , Fa 1 se , 

XtNs ens i t i ve .False, 

XtNlabel, ■ ' , 

NULL) ; 

} 

for ( i=0 ; i<8 ; i++) { 

sprintf (s , B FEItem_2_%d* , i) ; 

FEItem[2 ] [i] = XtVaCreateManagedWidget ( 

s , toggleWidgetClass , FEBox [ 2 ] , 

XtNresize, False, 

XtNsensitive, False, 

XtNlabel, - ■ , 

NULL) ? 

) 

for ( i = 0 ; i<8 ; i++) { 

sprintf (s, *FEItem_3_%d- , i) ; 

FEItem[3][i] = XtVaCreateManagedWidget ( 
s , toggleWidgetClass , FEBox [31, 

XtNresize, False, 

NULL) ; 

} 

for ( i=0 ; i<2 ; i++) { 

sprintf (s, *FEItem_4_%d" , i) ; 

FEItem[4] [i] = XtVaCreateManagedWidget ( 
s , commandWidgetClass , FEBox [4 ] , 

NULL) ; 

XtAddCal lback ( FEI tern ( 4 ] [i] , XtNcallback, FEPageCallback, (XtPointer) i) 

) 

for ( i=0 ; i<7 ; i++ ) { 

sprintf (s, -FEItem_5_%d- , i) ; 

FEItem[5][i] * XtVaCreateManagedWidget ( 
s , labelWidgetClass , FEBox [ 5 ] , 

XtNjustify , XtJustifyRight , 

NULL) ; 

} 

for ( i=0 ; i<7 ; i++) { 
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sprintf ( s , "FEltem_6_%d* , i) ; 

FEItem[ 6] [ i] - XtVaCreateManagedWidget ( 

s , asci iTextWidgetClass , FEBox [ 6 ] , 

XtNeditType , XawtextEdit , 

NULL) ; 

for (i=0;i<7;i>+) ( 

sDrintf (s, ■FEItem_7_%d" , i) ; 

FEItem[7]tiJ = XtVaCreateManagedWidget ( 
s , asci iTextWidgetClass, FEBox [7 ] , 

XtNeditType , XawtextEdit , 

NULL) ; 

) • 

'* FEItem[81 [0] = XtVaCreateManagedWidget^ 

■FEItem_8_0" , asci iTextWidgetClass, FEBox [8] , 

XtNed i tType , XawtextEdit , 

NULL) ; 

*/ 

for (i=l; i<6;i++) ( . 

SDrintf (s , ■FEItem_8_% < 3 , l) ; 

FEItem[8) [i] = XtVaCreateManagedWidget < 

s , commandWidgetClass , FEBox [ 8 ) , 

XtAddCallback (FEItem[8] [i] .XtNcallback, 

FECommandCallback, (XtPointer) x) ; 

r Create the widgets that I need to make the active plot dialog V 
PMShell = XtVaCreatePopupShelK _ 

■PMShell- , transientShellWidgetClass , topLevel , 

XtNtitle, *Plot Manager*, 

PMForm^ = XtVaCreateManagedWidget ( 

■ pfff'orm" , f orroWidgetClass , PMShell , 

NULL) ; 

for (i=0;i<4;i++) { 

sprintf (s, * PMBox_%d , i) ; 

PMBox [ i ] = XtVaCreateManagedWidget t 

s , boxWidgetClass , PMForm, 

XtNhSpace, 1, 

XtNvSpace , 1 , 

XtNtop, XtChainTop , 

NULL) ; 

sprintf (s, *PMTop_%d’ ,i) ; 

PMTop [ i ] = XtVaCreateManagedWidget { 

s , labelWidgetClass , PMBox ( i j , 

NULL); 

XtVaSetValues (PMBox [0] . XtNlef t , XCChainLef t , HULL I ; 

‘ I PMBox 1 i ] . xt»£ro»Hor lx , PMBox I i-1 ) , HULL) ; 

for ( i=0 ; i<8 ; i++ ) { 

sprintf ( s , ‘ PMToggle_%d* ,i); 
sprintf ( t , "Plot #%d*,i); 

PMitem[0] t i ] = XtVaCreateManagedWidget { 

s , toggleWidgetClass , PMBox [ 0 ] , 

XtNsensitive , False , 

XtNradioData, (caddr_t) 1 + 1 , 

XtNradioGroup, PMitemlO] [0] , 

XtNlabel , t , 

XtAddCallback (PMitem[0] til .XtNcallback, 

PMActiveCallback, (XtPointer) i) ; 

1 ' XawTogg 1 eChangeRadi oGr oup ( PMi t em [ 0 ] t i] . PMitem[0] [0] ) ; 

) 
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for (i=0;i<8;i++) { 

PMitem [ 1 ] [ i 3 = XtVaCreateManagedWidget ( 

■ entry" , commandWidgetClass , PMBox [ 1 J , 

XtNresize, False, 

XtNsensitive, False, 

XtN label, "Create" , 

NULL) ; 

XtAddCallback (PMitem[ 1 ] [i] , XtNcallback, 

PMCreateCallback, (XtPointer) i) ; 

) 

for (i=0;i<8;i++) { 

PMi tem[ 2 ] [ i ) = XtVaCreateManagedWidget ( 

" entry" , commandWidgetClass , PMBox [ 2 ) , 

XtNsensitive, False, 

XtNresize, False, 

XtN label, "Show" , 

NULL) ; 

XtAddCallback (PMitemf 2) [ i ) , XtNcallback, 

PMShowHideCal lback , (XtPointer) i) ; 

) 

for ( i-0 ; i<l ; i++) { 

sprint f (s , " PMitem_3_%d" , i ) ; 

PMitem[3][i] = XtVaCreateManagedWidget { 
s , commandWidgetClass , PMBox [ 3 ] , 

XtNresize, False, 

NULL) ; 

XtAddCallback ( PMi tem[ 3] [ i ] , XtNcallback, PMPlanarCal lback, 0 ) 

) 

for ( i = l ; i<7 ; i++ ) { 

sprintf {s , " PMitem_3_%d" , i) ; 

PMitem[3] [i] = XtVaCreateManagedWidget { 
s, commandWidgetClass , PMBox [3] , 

XtNresize, False, 

NULL) ? 

XtAddCallback {PMi tem[3] [i] , XtNcallback, 

PMActionCal lback, (XtPointer) i); 


execShell = XtVaCreatePopupShell ( 

"execShell " , transientShellWidgetClass , topLevel , 
XtNtitle, "Execution Summary", 

NULL) ; 

execForm = XtVaCreateManagedWidget ( 

• execForm" , formWidgetClass , execShell , 

NULL) ; 

for { i=0 ; i<3 ; i++ ) { 

sprintf (s, "execBox_%d* , i) ; 
execBox [ij = XtVaCreateManagedWidget ( 
s , boxWidgetClass , execForm, 

XtNhSpace, 1, 

XtNvSpace , 1 , 

XtNtop , XtChainTop , 

NULL) ; 

sprintf (s , "execTop_%d " , i ) ; 
execTop [ i ] « XtVaCrea ceManagedWidget ( 
s, labelWidgetClass, execBox[i] , 

NULL) ; 


) 

for ( i=l ; i<3 ; i++ ) 

XtVaSetValues ( execBox [ i ) , XtNfromHoriz , execBox[i-l] , NULL) ? 


for ( i=0 ; i<4 ; i++ ) { 

sprintf ( s , " execltem_0_%d" , i ) ; - — , 
execltem[0) (i) = XtVaCreateManagedWidget ( 

s , labelWidgetClass , execBox ( 0 ] , 
XtNresize, False, 

XtNjustify, XtJustifyRight , 

NULL) ; 
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for (i=0; i<4;i++) ( 

sprintf ( s , - execltem_l_%<*" , i) ; , M 

exedtemUl [il = XtVaCreateManagedWidget ( 
s , asciiTextWidgetClass , execBox [ 1 J , 

XtNeditType, XawtextEdit , 

NULL) ; 

for (i=0;i<3;i++) { 

sprintf (s, -execltem_2_%d ,i); . 

execltem[2] [i] = XtVaCreateManagedWidget { 

s , commandWidgetClass , execBox [ 2 ] , 

XtNresize , False, 

XtNjustify, XtJustifyCenter , 

XtAddCal lback ( exec I tem [ 2 ] [i] ,XtNcallback,ExecCallback, (XtPointer) i) 

) 

f ileShell = XtVaCreatePopupShe 1 1 ( 

■ f ileShell" , transientShellWidgetClass, topLevel, 

XtNtitle, ’Input Data’, 

f ilePane = XtVaCreateManagedWidget ( 

■ f ilePane 1 , panedWidgetClass # f ileShell , 

f i^eMenuPane = XtVaCreateManagedWidget ( 

" f ileMenuPane “ , panedWidgetClass , f ilePane, 

XtNor ientation , Xtor ientHorizontal , 

XtNshowGr ip , False , 

XtNskipAdjust , True, 

f ileEntrv ( 0) = XtVaCreateManagedWidget { 

■ f ileEntry_0 “ , commandWidgetClass , f ileMenuPane, 

XtN label , -Done" , 

XtNshowGr ip , False , 

fil^Ent^yll] = XtVaCreateManagedWidget ( 

* f ileEntry_l “ , commandWidgetClass , filer 
XtNlabel , -Save" , 

XtNshowGr ip , False , 

f ileEntrv [2 ] = XtVaCreateManagedWidget { 

" f ileEntry_2 * , commandWidgetClass, f ileMenuPane, 

XtNlabel, “Load - , 

XtNshowGrip , False , 

NULL) ; 

£ °XtAddciuMck(£ileB>try(il ,X t Nc.llb,ck.Fil e cllb»ck. (Xtpointkr. i) , 
f ileEntrv [3] = XtVaCreateManagedWidget ( 

• f ileEntry_3 " , asciiTextWidgetClass , f ileMenuPane, 

XtNedi tType , XawtextEdit , 

XtNshowGrip , False , 

NULL) ; 

fileText = XtVaCreateManagedWidget < 

- fileText- , asciiTextWidgetClass , f ilePane, 

XtNwidth, 480 , 

XtNheight, 240, 

XtNedi tType, XawtextEdit, 

XtNscrollVertical , XawtextScrollAlways , 

XtNshowGrip , False , 

NULL) ; 

/. create the widgets for the motor geometry dialog */ 

MGShell = XtVaCreatePopupShell ( 

•MGShell 1 , transientShellWidgetClass , topLevel, 

XtNtitle, "Solid Rocket Motor Geometry , 

NULL) ; . , 

MGForm = XtVaCreateManagedWidget ( 


, f ileMenuPane, 
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•MGForm" , f ormWidgetClass , MGShell , 
NULL) ; 

for (i=0; i<6; i++) { 

sprint f (s, "MGBox_%d" , i) ; 

MGBox [ i ] * XtVaCreateManagedWidget ( 
s , boxWidgetClass , MGForm, 

XtNhSpace , 1, 

XtNvSpace , 1 , 

NULL) ; 

sprintf { s , *MGTop_%d* , i) ; 

MGTopCi] = XtVaCreateManagedWidget { 
s , labelWidgetClass , MGBox [ i ] , 


NULL) ; 

XtVaSetValues (MGBox [ 0 ] , XtNlef t , XtChainLef t , XtNtop # XtChainTop , NULL) ; 
for ( i=l ; i<3 ; i++ ) 

XtVaSetValues (MGBox [i] , XtNfromHoriz, MGBox [i-1] , NULL) ; 

XtVaSetValues (MGBox ( 3 ] , XtNf romVert , MGBox [01 , XtNlef t , XtChainLef t , NULL) 
for ( i=4 ; i<6 ; i++) 

XtVaSetValues (MGBox [ i ) , XtNfromHoriz , MGBox ( ll ] , 

XtNf romVert , MGBox [ 0 3 , NULL) ; 


MGPictShell = XtVaCreatePopupShell ( 

•MGPictShell * , transientShellWidgetClass , topLevel , 
XtNminWidth, 240 , 

XtNminHeight , 240 , 

XtNmaxAspectX, 1 , 

XtNmaxAspectY , 1 , 

XtNsaveUnder , False, 

XtNtitle, "XXX", 

NULL) ; 

MGPict = XtVaCreateManagedWidget ( 

•MGPict • , widgetClass , MGPictShell , 

XtNwidth, 360 , 

XtNheight, 360, 

NULL) ; 


for ( i=0 ; i<8 ; i++) { 

sprintf ( s , "MGI tem_0 — %d “ , i ) ; 

MGItem[0][i] = XtVaCreateManagedWidget ( 
s, labelWidgetClass, MGBox [0] , 

XtNresize, False, NULL) ; 

) 

for { i=0 ; i<8 ; i++ ) { 

sprintf (s, "MGItem_l_%d* , i ) ; 

MGItem(l) [i] = XtVaCreateManagedWidget ( 
s , asciiTextWidgetClass , MGBox [ 1 ] , 
XtNeditType, XawtextEdit , NULL) ; 

} 

for ( i=0 ; i<6 ; i++ ) { 

sprintf (s, "MGItem_2_%d' , i) ; 

MGItem[2)(i] = XtVaCreateManagedWidget ( 
s r toggleWidgetClass , MGBox [2 ) , 
XtNradioGroup,MGItem[2 ] [0] f 
XtNradioData, (caddr_t) i+1, 

XtNresize, False, NULL) ; 

XtVaSetValues (MGItem[2j [0] , XtNst ate, True, NULL) ; 

for ( i=0 ; i<2 ; i++ ) { 

sprintf ( s , -MGItem_3_%d* , i ) ? 

MGItem[3] [i] = XtVaCreateManagedWidget ( 
s , labelWidgetClass , MGBox [ 3 ] , 

XtNresize, False, NULL) ; 

} 

for { i = 0 ; i<2 ; i++ ) { , 

sprintf (s, -MGIten\_4_%d- , i) ? 

MGItem[4] ( i ] = XtVaCreateManagedWidget ( 
s , asciiTextWidgetClass , MGBox [4 ] , 
XtNeditType, XawtextEdit , NULL) ; 
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for { 1=0 ; i<4 ; i++) ( 

sprintf (s, •MGItem_5_%d , ‘ , i) ? 

MGItem[5] [i] * XtVaCreateManagedWidget! 

s , commandWidgetClass , MGBox [ 5 ] ,NULL) ; 
XtAddCallback (MGItem [ 5 ] [i] , XtNcallback, 


MGCallback, (XtPointer) 


} 


i) 


/* Create the widgets for the curve fitting dialog V 


CFShell = XtVaCreatePopupShell ( 

•CFShell " , transientShellWidgetClass , topLevel , 
XtNtitle, "Data Curve Fitting", 

NULL); , ■ 

CFForm = XtVaCreateManagedWidget! 

■CFForm" , f ormWidgetClass , CFShell , 

NULL) ; 

for { i=0 ; i<3 ; i++ ) { 

sprintf (s, -CFBox_%d" , i) ; 

CFBox [ i ] = XtVaCreateManagedWidget ( 

s , boxWidgetClass , CFForm, 

XtNhSpace , 1 , 

XtNvSpace, 1, 

XtNtop , XtChainTop , 

NULL) ; 

sprintf {s, •CFTop_%d’' » i) ? 

CFTop [ i 1 = XtVaCreateManagedWidget ( 

s , labelWidgetClass , CFBox [ i ] , 

NULL) ; 


i % 

for (i=l;i<3;i++) x 

XtVaSetValues (CFBox [ i ] , XtNf romHoriz , CFBox [ l-l] # NULL) ; 


for (i=0;i<16;i++) { 

sprintf (s, •CFItemJ)_%d‘ , i) ; 

CFItem[0] [i] = XtVaCreateManagedWidget! 

s , toggleWidgetClass , CFBox [ 0 ] , 

XtNlabel , • " , 

XtNresize , False, 

XtNsensitive, False, NULL) ; 


for (i=0;i<16;i++) { 

sprintf (s, "CFItem_l_%cl" , i) ; 

CFItem[ 1 ) [ i] = XtVaCreateManagedWidget ( 
s , asciiTextWidgetClass , CFBox [ 1 ] , 
XtNeditType, XawtextEdit , NULL) ; 

for (i=0? i<4? i++) ( 

sprintf (s, *CFItem_2_%<3" ,i) ; 

CFItem [ 2 ] [ i ] = XtVaCreateManagedWidget! 

s,commandWidgetClass,CFBox[2] , 

NULL) ; 

) 


} 
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